《BMC Microbiology》：Conserved protein folds underpin the diversification of secreted proteins in a fungal pathogen
编辑推荐：
在宿主定殖过程中，植物病原真菌分泌效应蛋白样蛋白以改变宿主细胞生理并靶向植物关联微生物。然而，效应蛋白样蛋白的快速进化与低序列保守性阻碍了其研究与表征。真菌Zymoseptoria passerinii侵染Hordeum属植物，包含适应野生与驯化大麦的谱系，迄
在宿主定殖过程中，植物病原真菌分泌效应蛋白样蛋白以改变宿主细胞生理并靶向植物关联微生物。然而，效应蛋白样蛋白的快速进化与低序列保守性阻碍了其研究与表征。真菌Zymoseptoria passerinii侵染Hordeum属植物，包含适应野生与驯化大麦的谱系，迄今该物种效应蛋白样蛋白的进化尚未被探讨。研究人员结合多种基于结构与网络的分析解析Z. passerinii的分泌组。首先比较AlphaFold2(AF2)与ESMFold预测以建立结构分析基线；在分泌组中鉴定出72个结构簇，揭示发散序列间的折叠水平关系。研究发现具预测宿主免疫干扰功能的效应蛋白样蛋白源自有限的一组蛋白折叠，而具预测抗菌性质的蛋白分布于不同折叠组。理化比较表明，推定性抗菌效应蛋白主要在Z. passerinii常见效应富集支架上通过氨基酸替换涌现，重构表面电荷与静电势。研究人员通过比较Z. passerinii蛋白与Zymoseptoria属同源物分析选定效应富集家族的种内与种间变异，描述受约束的核心折叠及环区与表面暴露区的局部变异，这与折叠稳定性一致同时允许蛋白多样化。推定性抗菌效应蛋白同源物尽管序列发散仍广布于该属。结论：Z. passerinii分泌组围绕共同结构折叠组织，支撑包括宿主操控与宿主关联微生物互作在内的多样生物角色；保守支架结合表面与理化变异可能促进该真菌效应蛋白样蛋白的快速适应性进化。
研究背景与立项依据
植物病原真菌在宿主定殖中分泌效应蛋白样蛋白（effector-like proteins），可抑制模式触发免疫、扰动活性氧信号、改变激素平衡并靶向宿主成分，同时部分效应蛋白样蛋白具备抗菌（antimicrobial）活性，可在宿主组织内抑制微生物竞争者或重塑叶际微生物组。由于多数真菌分泌蛋白体积小、进化快、缺乏保守序列域，传统基于序列的方法难以推断远缘分类单元间的同源与功能关系。结构生物学研究显示，即便序列不相似，效应蛋白可共享三维折叠，因此基于结构的方法能揭示序列比对无法发现的同源信号。真菌Zymoseptoria属中Z. tritici研究较多，而其近缘种Zymoseptoria passerinii（引起大麦Septoria speckled leaf blotch）包含感染野生与驯化大麦的遗传分化谱系，先前研究虽表征其推定效应子库并发现独特效应子，但分泌蛋白在结构水平的功能多样化与进化机制未受关注。研究人员以分离自野生大麦Hordeum murinum subsp. glaucum的高质量长读长基因组代表分离株Zpa796为模型，开展分泌组规模的结构基因组学分析，旨在明确折叠水平关系是否支撑效应与抗菌功能的分化。
论文发表于《BMC Microbiology》，其意义在于将结构预测、折叠聚类与理化性质分析整合入野生草病原系统，证明保守折叠支架辅以表面电荷重排是效应蛋白样蛋白快速适应的基础。
主要关键技术方法（≤250字）
研究人员基于Z. passerinii Zpa796（源自野生大麦H. murinum subsp. glaucum）预测743个分泌蛋白，并行使用AlphaFold2(AF2)与ESMFold进行结构预测并以pLDDT、TM-score、RMSD评估；以双向TM-score>0.5构建结构相似网络，得72个结构簇；用EffectorP 3.0与AMAPEC分别预测效应与抗菌属性；以CATH、SCOPe、ECOD、PDB、AlphaFoldDB做折叠注释；对选定簇（G.04、G.12）做多结构比对、二硫键推断、密码子水平选择压力分析（纯化/多样化选择）；用OrthoFinder在Zymoseptoria属多物种中找同源并建最大似然系统树；跨物种选Z. tritici IPO323同源物做结构叠加与理化比较。
研究结果
比较AlphaFold2与ESMFold的结构预测
研究人员对Zpa796共743个分泌蛋白建模，AF2成功670个、ESMFold 742个。AF2平均pLDDT 78.23高于ESMFold的63.03，二者pLDDT正相关（r=0.79）。以PDB中228个真菌晶体结构为基准，AF2平均TM-score 0.95、RMSD 0.93 ?，显著优于ESMFold。短肽（<100 aa）两者差异最大。结论：AF2置信度与精度更高，后续采用670个AF2模型。
基于结构的注释拓展分泌组功能
408个蛋白获序列注释、479个获结构注释；402个（60.5%）AF2模型以TM-score≥0.5匹配PDB/CATH/SCOPe/ECOD，另77个命中AlphaFoldDB。EffectorP预测310个效应子、AMAPEC预测341个抗菌蛋白(AMP)，141个重叠为推定性抗菌效应蛋白；结构注释比序列注释多覆盖约68个效应子与53个AMP。
结构相似网络
以双向TM-score>0.5建网，364节点（94% pLDDT>70）分72个连通组件（结构簇），最大簇36蛋白，26个簇仅2蛋白，306个为单例。簇内成员在CATH topology/SCOPe fold/ECOD X-group一致。
结构簇与单例的功能图景
107个推定效应子入31簇，224个AMP入63簇，79个抗菌效应蛋白跨簇分布。6个效应富集簇（p-adj<0.05）：G.03 Gamma-crystallin、G.04 Alpha-beta plaits/Ferredoxin-like、G.06 β-trefoil、G.12 Killer toxin KP4、G.14 Allergen Alt a1、G.21未表征；无AMP富集簇。AMP比例高的簇多具酶活（G.08 Peptidase S8、G.09 Chloroperoxidase、G.11 Acid protease、G.28 Peptidase S28）。单例218个（71%）无分层库注释，多伴无序残基>40%。
网络内理化性质
按功能分四组（抗菌效应79、仅效应28、仅AMP145、其他112），六属性（长度、净电荷pH7、电荷密度、螺旋疏水矩、平均疏水性、表面疏水性）Kruskal-Wallis均显著。抗菌效应与仅效应组最短、绝对净电荷最高；归一化后抗菌效应电荷密度高于仅效应。仅AMP组表面疏水性低而疏水矩/平均疏水性高。簇间比较（46个≥3成员簇）同属性差异显著，同簇内属性跨度大。
保守结构支架支撑序列多样化
选G.04（α+β，核心jg2244.t1为β-α-β-β-α-β，3对二硫键C10-C55、C22-C40、C41-C68）与G.12（α-β-β-α-β-β，核心jg1.t1，C2-C84、C23-C73、C43-C104）。结构对齐显示核心二级结构保守，环区（G.04的β1-α1、β2-β3；G.12的α1-β1）偏离最大。序列同一性低，G.04仅4半胱守恒，G.12仅1半胱1色氨酸守恒；jg1.t1与jg6868.t1 RMSD 0.46、TM-score 0.986但非保守替换多。密码子选择分析：纯化选择位于二级结构与二硫半胱位点，多样化选择位于表面（G.04 α1一处、G.12环区三处）。结论：二硫稳定核心+表面替换=支架上功能发散。
Z. passerinii与Zymoseptoria属同源进化
364网内蛋白359个在属内有同源；G.04/G.12同源广布属内且种间复制数不同。选jg9741.t1（G.04）对Z. tritici IPO323平行源ZtIPO323_014270.1/ZtIPO323_104450.1，jg3663.t1（G.12）对ZtIPO323_026410.1做结构叠加：折叠高度保留，差异仅在表面环与末端延伸；跨种同源净电荷与pI不同——Zpa796同源物在微酸质外体更偏正电，jg9741.t1螺旋α1两亲性疏水面更强。结论：折叠跨种保守，表面理化随宿主与谱系特异化。
讨论与结论翻译
讨论指出，Z. passerinii分泌组中>45%预测为抗菌蛋白，反映微生物互作在其生活史中被低估的重要性；许多推定效应子兼具抗菌性，构成后续功能验证优选集。AF2对含保守域蛋白置信优于ESMFold，短肽常判为无序。
结论原文浓缩：Z. passerinii分泌组围绕共同结构折叠组织，支撑宿主操控与宿主关联微生物互作等多样角色；保守支架结合表面与理化变异可能促进该真菌效应蛋白样蛋白的快速适应性进化。
