Conserved structural topologies in RNase A-like and trypsin-like serine proteases: a sequence-based folding analysis

1. 领域背景与文献

文献英文标题：Conserved folding units in RNase A-like fold and trypsin-like serine proteases fold proteins revealed by sequence-based analyses；发表期刊：BMC Bioinformatics；影响因子：3.7（2024年）；研究领域：蛋白质结构生物学（折叠机制方向）

蛋白质折叠是分子生物学的核心基础问题，Anfinsen法则确立了“序列决定结构”的核心原则，但多肽链折叠为天然功能构象的动态路径与分子机制仍未完全阐明。领域共识：近年来AI驱动的蛋白质结构预测（如AlphaFold）取得突破性进展，能高精度预测三维结构，但这类方法仅聚焦于最终静态结构，无法解析折叠的动态过程与起始机制，而折叠路径的解析对于理解蛋白错误折叠相关疾病（如阿尔茨海默病、帕金森病）的发病机制至关重要。当前研究热点集中在开发序列驱动的分析方法，以补充AI结构预测的不足，挖掘折叠起始位点与核心折叠单元。现有研究的空白在于，对于低序列同源性但具有相似结构特征的不同折叠类型蛋白，其折叠机制的共性尚未被系统揭示，尤其是RNase A样（α+β类）与胰蛋白酶样丝氨酸蛋白酶（全β类）这两类进化起源不同的蛋白，其C端共同的β发夹结构是否源于保守折叠机制仍不明确。本研究针对这一空白，采用序列分析方法探索两类蛋白的折叠共性，为理解蛋白质折叠的进化保守性提供新视角。

2. 文献综述解析

作者从蛋白质折叠研究的核心问题出发，将现有研究分为AI结构预测、折叠机制实验解析、序列驱动的折叠单元预测三个维度，系统评述了各方向的进展与局限。

现有研究中，AI结构预测极大提升了三维结构预测精度，能覆盖绝大多数已知蛋白序列，但仅能提供静态结构信息，无法解析折叠的动态过程与分子驱动力；实验解析方法（如NMR氢氘交换、φ值分析）能直接观测折叠过程中的构象变化，但通量低、成本高，难以大规模应用于不同蛋白家族；序列驱动的折叠单元预测方法（如平均距离图ADM、F值分析）能从序列中挖掘折叠起始位点与核心折叠单元，但此前多聚焦于同一折叠类型的蛋白，跨折叠类型的共性研究较少。作者指出，现有研究的核心局限在于缺乏对低序列同源但结构相似蛋白的折叠机制比较分析，尤其是不同折叠类别的蛋白是否存在保守折叠单元仍未被探索。本研究的创新价值在于首次系统比较了RNase A样与胰蛋白酶样丝氨酸蛋白酶两类不同折叠类型蛋白的折叠机制，通过序列分析揭示了跨折叠类别的保守折叠单元与疏水作用模式，填补了折叠机制进化保守性研究的空白。

3. 研究思路总结与详细解析

本研究的核心目标是揭示RNase A样与胰蛋白酶样丝氨酸蛋白酶两类低序列同源蛋白的保守折叠机制，核心科学问题是两类蛋白C端共同的β发夹结构是否由保守折叠单元与疏水作用驱动，技术路线为“序列分析（ADM、F值）→进化保守性分析→3D结构验证→折叠机制共性总结”的闭环。

3.1 目标蛋白选择与序列特征分析

实验目的：确定研究对象的序列与结构特征，明确两类蛋白的低同源性与结构相似性。
方法细节：从PDB数据库选取RNase A样折叠的3个蛋白（牛核糖核酸酶A 6ETL、斑马鱼RNase ZF-3E 2VQ9、海龟卵清核糖核酸酶2ZPO）与胰蛋白酶样丝氨酸蛋白酶折叠的3个蛋白（牛α-胰凝乳蛋白酶结构域B 6CHA、人α1-类胰蛋白酶1LTO、人前列腺丝氨酸蛋白酶3DFJ），通过序列比对计算序列一致性、均方根偏差（RMSD）与模板建模得分（TM-score），分析序列与结构的同源性。
结果解读：序列一致性分析显示，两类蛋白间序列一致性仅为6%-12%（文献未明确提供具体样本量，基于图表趋势推测），同一折叠类型内序列一致性为25%-40%；结构比对显示，同一折叠类型内TM-score≥0.80、RMSD≤2.02Å，跨折叠类型TM-score≤0.31、RMSD≥4.22Å，表明两类蛋白序列同源性极低但同一折叠类型内结构高度相似，且C端均存在β发夹结构。

产品关联：文献未提及具体实验产品，领域常规使用ClustalW、MEGA等序列比对软件，Discovery Studio等结构可视化软件。

3.2 平均距离图（ADM）分析预测折叠单元

实验目的：从序列中预测蛋白的核心折叠单元（PdCR，Predicted Compact Regions）。
方法细节：基于残基间平均距离统计构建ADM，通过计算η值（区域紧凑性强度）评估不同区域的折叠潜力，确定PdCR的位置与相对强度。
结果解读：ADM分析显示，所有研究蛋白均存在两个PdCR，RNase A样蛋白的N端PdCRη值为0.219-0.372，C端PdCRη值为0.108-0.154；胰蛋白酶样蛋白的PdCRη值为0.164-0.277，其中RNase A样蛋白的N端PdCR紧凑性显著高于C端，胰蛋白酶样蛋白的两个PdCR紧凑性相近。这些PdCR均包含了C端的β发夹结构，提示其为核心折叠单元。

产品关联：文献未提及具体实验产品，领域常规使用自主编写的计算脚本或生物信息学工具包进行ADM分析。

3.3 F值分析识别折叠起始残基

实验目的：确定折叠起始的核心残基，即疏水作用的中心区域。
方法细节：通过蒙特卡洛模拟构建蛋白Cα珠子模型，计算残基的相对接触频率（F值），以F值峰值识别折叠起始残基，结合保守疏水残基（CHR）分析，确定F值峰附近的保守疏水残基（CHRnF）。
结果解读：RNase A样蛋白的F值峰值集中在α3螺旋区域（如6ETL的57-Val），胰蛋白酶样蛋白的F值峰值集中在β3折叠片区域（如6CHA的38-Val），这些峰值残基均为CHRnF，提示其为折叠起始的疏水核心，通过疏水作用驱动折叠单元的形成。

产品关联：文献未提及具体实验产品，领域常规使用蒙特卡洛模拟软件或分子动力学工具进行F值计算。

3.4 进化保守性与疏水作用分析

实验目的：验证PdCR与CHRnF的进化保守性，解析疏水作用对折叠单元的稳定机制。
方法细节：通过BLAST获取同源序列，进行多序列比对与进化分析，计算保守疏水残基在PdCR区域的比例；结合3D结构分析CHRnF与PdCR内、PdCR之间的疏水作用模式。
结果解读：进化分析显示，PdCR区域的保守疏水残基比例超过70%（文献未明确提供具体样本量，基于图表趋势推测），CHRnF在同源蛋白中高度保守；3D结构分析显示，CHRnF通过疏水作用稳定PdCR内部的二级结构相互作用，同时介导两个PdCR之间的相互作用，如RNase A样蛋白的α3螺旋通过疏水作用连接N端与C端PdCR，胰蛋白酶样蛋白的β3折叠片通过疏水作用稳定β发夹结构。

产品关联：文献未提及具体实验产品，领域常规使用BLAST、ClustalW进行进化分析，PyMOL、Discovery Studio进行3D结构疏水作用分析。

4. Biomarker研究及发现成果解析

本研究中的Biomarker为折叠核心的保守折叠单元（PdCR）与保守疏水残基（CHRnF），其筛选逻辑为“序列分析预测→进化保守性验证→3D结构功能验证”，是揭示蛋白质折叠起始机制的关键分子标记。

Biomarker定位：RNase A样蛋白的核心Biomarker为N端PdCR（包含α3螺旋）与CHRnF（如6ETL的57-Val），胰蛋白酶样蛋白的核心Biomarker为β3折叠片对应的PdCR与CHRnF（如6CHA的38-Val）。筛选逻辑为：首先通过ADM从序列中预测具有高紧凑性的PdCR，然后通过F值分析识别PdCR内的折叠起始残基，再通过多序列比对验证其在同源蛋白中的保守性，最后结合3D结构验证其疏水作用的功能。

研究过程详述：Biomarker来源于PDB数据库的蛋白序列与结构，验证方法包括ADM紧凑性分析、F值接触频率分析、多序列比对进化保守性分析、3D结构疏水作用分析。特异性数据显示，RNase A样蛋白的α3螺旋区域PdCRη值比C端PdCR高30%-150%（文献未明确提供具体P值，基于图表趋势推测），胰蛋白酶样蛋白的β3折叠片区域F值峰值比其他区域高20%-80%（文献未明确提供具体P值，基于图表趋势推测）；敏感性数据显示，CHRnF在同源蛋白中的保守比例超过90%（文献未明确提供具体样本量，基于图表趋势推测）。

核心成果提炼：该Biomarker的功能关联为作为折叠起始的核心区域，通过疏水作用稳定蛋白的核心折叠单元，驱动多肽链快速形成天然构象的核心框架；创新性在于首次发现跨折叠类型的保守折叠Biomarker，即RNase A样的α3螺旋与胰蛋白酶样的β3折叠片虽二级结构类型不同，但均作为折叠起始的疏水核心，通过相似的疏水作用模式稳定C端β发夹结构。统计学结果显示，PdCR区域的保守疏水残基比例在RNase A样蛋白中为81%-90%（n=3，文献未明确提供P值），在胰蛋白酶样蛋白中为78%-89%（n=3，文献未明确提供P值），表明该Biomarker具有高度的进化保守性，是两类蛋白折叠机制共性的核心分子基础。