植物半胱氨酸蛋白酶抑制剂支架作为结构模块的设计与应用研究

《Proteins: Structure, Function, and Bioinformatics》：The Supporting Scaffold as a Structural Module for Plant Cystatin Design

　　

摘要：蛋白质工程方法已被提出用于改善半胱氨酸蛋白酶抑制剂（cystatins）对植物害虫消化半胱氨酸蛋白酶（Cys proteases）的抑制效力，通常涉及已知与靶蛋白酶相互作用的N末端主干（N-terminal trunk）或抑制性环（inhibitory loops）中的点突变。本研究评估了环支撑框架（即“支架”（scaffold））是否能成为半胱氨酸蛋白酶抑制剂改进中的一个有价值的结构模块。研究人员通过计算机模拟设计了杂合半胱氨酸蛋白酶抑制剂（hybrid cystatins），将供体半胱氨酸蛋白酶抑制剂的N末端主干和抑制性环移植到替代受体半胱氨酸蛋白酶抑制剂的支架上。根据AlphaFold模型推断，支架替换对半胱氨酸蛋白酶抑制剂的折叠和抑制性环的空间取向影响很小，但在若干情况下改变了N末端主干的空间定位以及蛋白酶-半胱氨酸蛋白酶抑制剂复合物的预测结合强度。在大肠杆菌（E. coli）中产生的合成版本杂合体被评估其对木瓜蛋白酶（papain）以及一种草食性昆虫模型消化蛋白酶的抑制效力。在多个案例中，支架替换对所产生半胱氨酸蛋白酶抑制剂的抑制效力产生了影响。例如，携带模式番茄半胱氨酸蛋白酶抑制剂SlCYS8的N末端主干和抑制性环、并移植到其他半胱氨酸蛋白酶抑制剂支架上的半胱氨酸蛋白酶抑制剂，其对木瓜蛋白酶的Ki值变化高达3倍，而对昆虫蛋白酶的抑制率提高多达8倍。研究数据表明支架对植物半胱氨酸蛋白酶抑制剂的抑制活性具有显著影响。这些结果也提示了该结构元件作为植物半胱氨酸蛋白酶抑制剂设计模块的潜力，可用于产生针对草食性害虫消化蛋白酶的功能多样性。

论文解读：植物半胱氨酸蛋白酶抑制剂支架作为结构模块的设计与应用

研究背景与目的

植物半胱氨酸蛋白酶抑制剂（cystatins）是一类能够抑制半胱氨酸蛋白酶（Cys proteases）活性的小分子蛋白，在植物防御草食性节肢动物中发挥重要作用。当害虫取食植物组织时，这些抑制剂进入昆虫消化道，抑制其消化蛋白酶，导致氨基酸缺乏，进而引起生长迟缓、发育缺陷甚至死亡。然而，害虫已经进化出多种应对策略，包括产生对抑制剂不敏感的蛋白酶亚型以及过量表达敏感蛋白酶以补偿消化功能的损失。因此，如何提升植物半胱氨酸蛋白酶抑制剂的抑制效力和拓宽其抑制范围，成为植物保护领域的关键挑战。

传统的蛋白质工程策略主要集中于对功能相关区域——即N末端主干（N-terminal trunk）和两个抑制性环（inhibitory loops）——进行定点突变或随机突变。这些区域直接参与与靶蛋白酶活性位点的相互作用。尽管这些方法取得了一定成效，但关于非功能性的支撑结构（即“支架”scaffold）是否也能作为一个独立的结构模块来调控抑制活性，尚缺乏系统性研究。本研究旨在评估支架替换策略是否能够在不破坏蛋白整体折叠的前提下，有效改变抑制剂的效力，从而为设计更高效的抗虫抑制剂提供新的思路。

研究方法概述

研究人员采用计算机模拟设计（in silico design）构建了20种杂合半胱氨酸蛋白酶抑制剂（hybrid cystatins）。具体而言，将来自四种供体半胱氨酸蛋白酶抑制剂（SlCYS8、PpCYS、StCYS5、CsCYS）的N末端主干和两个抑制性环分别移植到12种不同植物来源的受体支架上。所有杂合体通过大肠杆菌（E. coli）BL21菌株异源表达，利用谷胱甘肽S-转移酶（GST）基因融合系统进行亲和纯化。结构建模采用AlphaFold 3预测工具，结合均方根偏差（RMSD）、质心位移（centroid shift）等距离分析评估结构变化。结合自由能通过分子力学/广义玻恩表面积（MM/GBSA）方法估算。抑制活性测定针对木瓜蛋白酶（papain）和从科罗拉多马铃薯甲虫（Leptinotarsa decemlineata）四龄幼虫中肠提取的消化半胱氨酸蛋白酶（包括组织蛋白酶L样和组织蛋白酶B样酶）进行。

研究结果

3.1 半胱氨酸蛋白酶抑制剂支架替换的一般方案

研究人员设计了一种通用的支架替换方案，将供体半胱氨酸蛋白酶抑制剂的N末端主干和两个抑制性环（Loop 1、Loop 2）移植到受体支架上。共设计了20种杂合序列，涵盖四个供体和十二个受体。大多数杂合体成功在大肠杆菌中表达并纯化，仅有两个含有StCYS5功能元件的杂合体未能获得可溶性产物。

3.2 支架替换改变SlCYS8 N末端主干的空间定位

AlphaFold 3模型显示，所有12种受体支架均能与供体功能元件兼容，保持核心半胱氨酸蛋白酶抑制剂折叠和抑制性环的正确空间取向。全局RMSD值介于0.17至0.56 ?之间，远低于2.0 ?的结构变异阈值，表明整体结构高度相似。然而，N末端主干的RMSD值高达1.5至4.0 ?，显示出显著的位移或构象重排，而两个抑制性环则保持紧密聚集。

3.3 支架依赖的N末端主干定位受靶酶影响

在木瓜蛋白酶结合的复合物模型中，某些杂合体的N末端主干RMSD较游离状态减小（如EgCYS支架杂合体从3.5 ?降至1.1 ?），另一些则增大（如CsCYS支架杂合体从1.8 ?增至4.4 ?）。这表明支架替换对N末端主干的影响在蛋白酶结合后可能被放大或减弱，呈现出支架特异性。MM/GBSA计算显示，不同支架杂合体与木瓜蛋白酶的结合自由能（ΔG）差异显著，其中SlCYS7、HaCYS、JcCYS支架杂合体表现出比野生型SlCYS8更负的ΔG值，暗示更强的预测结合。

3.4 SlCYS8功能元件对木瓜蛋白酶的抑制效力呈支架依赖性

体外抑制实验证实，支架替换确实影响了SlCYS8功能元件对木瓜蛋白酶的抑制效力。最有效的杂合体（EgCYS和StCYS5支架）的K_i值分别降低至野生型的约1/2.1和1/1.7，而最弱的杂合体（ZmCYS和SiCYS支架）的K_i值增加至野生型的2.8倍和1.7倍。值得注意的是，一些高效支架来源于天然具有强抗木瓜蛋白酶活性的半胱氨酸蛋白酶抑制剂，但并非所有情况都符合这一规律。

3.5 SlCYS8功能元件对昆虫消化组织蛋白酶的抑制效力同样呈支架依赖性

针对科罗拉多马铃薯甲虫三种消化组织蛋白酶（IntD4、IntC14、IntB）的AlphaFold复合物模型显示，尽管整体排列与木瓜蛋白酶复合物相似，但MM/GBSA计算的ΔG值在不同支架杂合体间差异明显。体外抑制实验进一步证实，支架替换显著影响了SlCYS8功能元件对昆虫组织蛋白酶L样酶的抑制效力。最有效的杂合体（CsCYS、JcCYS、LjCYS、HaCYS支架）相对抑制率提高了4至10倍，而最弱的杂合体（PpCYS、SlCYS9、ZmCYS支架）几乎没有改善。有趣的是，高效支架往往来源于天然对昆虫蛋白酶活性较低的半胱氨酸蛋白酶抑制剂，而低效支架则来源于天然高活性抑制剂。对于通常对半胱氨酸蛋白酶抑制剂不敏感的组织蛋白酶B样酶，PpCYS和CsCYS支架意外地使抑制率提高了2至5倍。

3.6 天然支架影响植物半胱氨酸蛋白酶抑制剂的蛋白酶抑制效力

当将PpCYS和StCYS5的功能元件移植到替代支架上时，其对昆虫组织蛋白酶L样酶的抑制效力显著下降，表明这些天然高效抑制剂的支架对其高效性是必要的。相反，弱抑制剂CsCYS的功能元件即使移植到其他支架上也未能显著提升活性。这揭示了功能元件和支架共同决定最终抑制效力的协同关系。分布图分析显示，支架替换策略产生的功能变异系数（CV）为62%，与此前基于功能元件替换的策略（CV=61%）相当，且平均抑制率提高了近三倍。

结论与意义

本研究证明，植物半胱氨酸蛋白酶抑制剂的支架区域是一个有价值的结构模块，可用于设计和优化抑制剂的活性。AlphaFold模型和距离分析表明，支架替换对整体折叠和抑制性环的影响较小，但对N末端主干的定位有显著影响，进而影响蛋白酶-抑制剂复合物的稳定性。体外实验证实，支架替换可使SlCYS8功能元件对木瓜蛋白酶的K_i值变化高达3倍，对昆虫蛋白酶的抑制率提高多达8倍。研究结论指出：“支架区域对植物半胱氨酸蛋白酶抑制剂的抑制活性具有显著影响……这些结果也提示了该结构元件作为植物半胱氨酸蛋白酶抑制剂设计模块的潜力，可用于产生针对草食性害虫消化蛋白酶的功能多样性。”未来的研究需要进一步比较最优支架杂合体与已有功能元件杂合体及单突变体的植物保护效果，并评估其在转基因植物中的潜在脱靶效应。