蝎(Mesobuthus crucittii, Scorpiones: Buthidae)毒腺转录组及生物活性分析提示毒液中透明质酸酶具双功能(bifunctional)活性

《Frontiers in Molecular Biosciences》：Venom gland transcriptomics and bioactivity profiling suggest bifunctional hyaluronidase activity in the venom of Mesobuthus crucittii (Scorpiones: Buthidae)

摘要：蝎(Scorpion)毒液是具前景生物医学应用的生物活性分子的丰富来源。透明质酸酶(hyaluronidase)是毒液相关酶，可通过降解细胞外基质糖胺聚糖(glycosaminoglycan, GAG)促进毒素扩散，但其于蝎毒液中的结构多样性与底物特异性尚缺乏充分研究。本研究旨在鉴定并表征伊朗特有蝎种Mesobuthus crucittii毒腺中新型透明质酸酶(hyaluronidase)，并评估其生化与结构特性。方法：进行毒腺转录组分析；通过系统发育重建、基序(motif)预测、理化性质计算及结构建模对透明质酸酶(hyaluronidase)序列进行计算机(in silico)分析；开展分子对接(molecular docking)以探究其与透明质酸(hyaluronic acid, HA)及硫酸软骨素(chondroitin sulfate, CS)底物的相互作用；采用浊度法(turbidimetric assay)实验评估酶活性与热稳定性。结果：转录组分析揭示多样化毒素谱，以离子通道调节剂及酶组分为主，包括编码推定透明质酸酶(hyaluronidase)的新型前体，具独特半胱氨酸框架——六个二硫键及三个保守诊断基序(diagnostic motifs)(GDWW、FPDC和GWGS)；结构建模显示其具与糖基水解酶56家族(glycosyl hydrolase family 56, GH56)相符的催化域与结合域；分子对接支持其对透明质酸四糖(hyaluronic acid tetrasaccharide)而非更长或高度硫酸化GAGs具优先结合亲和力；粗毒功能测定确认强透明质酸(hyaluronan)降解能力及较慢硫酸软骨素(chondroitin sulfate)水解能力，提示具双底物活性。讨论：本研究支持存在结构独特的蝎毒透明质酸酶(hyaluronidase)，具推定双功能底物活性，拓展了对蝎毒酶进化的认知并凸显该酶的转化潜力。

《Frontiers in Molecular Biosciences》刊载论文解读：Mesobuthus crucittii毒腺转录组及透明质酸酶双功能活性研究

一、研究背景与立题依据

蝎(Scorpiones)隶属于节肢动物门螯肢亚门，全球已描述近2900种，中东及伊朗为物种多样性热点区域。蝎毒液主要由具多种生物学功能的蛋白与多肽组成，除作用于离子通道的神经毒素(neurotoxins)外，尚含金属蛋白酶(metalloproteases)、丝氨酸蛋白酶(serine proteases)、磷脂酶A₂(phospholipase A₂, PLA₂)、脂解激活肽(lipolysis activating peptides, LVPs)、抗菌肽(antimicrobial peptides, AMPs)及透明质酸酶(hyaluronidase, Hyal)等酶组分，是"毒液向药物(venoms-to-drugs)"研发的重要候选资源。

透明质酸酶(hyaluronidase, 亦称hyaluronoglucosaminidase)，属糖基水解酶家族(glycosyl hydrolase family, GH)，常见于动物毒液，可水解细胞外基质(extracellular matrix, ECM)中透明质酸(hyaluronic acid / hyaluronan, HA)及硫酸软骨素(chondroitin sulfate, CS)为寡糖，增加组织通透性助毒素扩散(spreading factor)，并可诱发炎症及IgE介导过敏反应。尽管蝎毒透明质酸酶(hyaluronidase)在蛰伤全身反应中的作用已有记载，但其生化全面表征及针对硫酸软骨素(chondroitin sulfate, CS)之活性仍缺乏深入探索。目前商业用透明质酸酶(hyaluronidase)源自牛/羊睾丸或重组人源，蝎毒代表潜在新型酶来源。为填补蝎毒透明质酸酶(hyaluronidase)结构多样性与双底物活性认知缺口，研究人员对伊朗特有种Mesobuthus crucittii（Buthidae科）开展毒腺转录组测序及透明质酸酶(hyaluronidase)的in silico与功能表征。

二、主要关键技术方法概述

研究人员采集伊朗胡齐斯坦省(Khuzestan)荒漠中Mesobuthus crucittii成虫，电刺激采集毒液诱导毒腺表达后解剖毒器(telson)提取总RNA，建库并采用Illumina HiSeq 2000平台进行RNA-seq(de novo转录组建集与CD-HIT聚类去冗余)，以Trinotate及tapai流程注释毒素组分比对动物毒素数据库(UniProtKB Animal Toxin Annotation Project)。从转录本中筛选透明质酸酶(hyaluronidase)序列进行MAFFT多序列比对(multiple sequence alignment, MSA)、IQ-TREE最大似然法(maximum likelihood, ML)构建系统发育树；利用InterProScan、SMART、Disulfide by Design 2.0预测保守基序(conserved motifs)、半胱氨酸及二硫键；Expasy ProtParam计算分子量、等电点(pI)、不稳定指数(aliphatic index, GRAVY)；SignalP 6.0预测信号肽；PTMGPT2预测N-糖基化位点(N-glycosylation sites)；以AlphaFold建立三维模型并经SAVES(ERRAT/Verify3D/PROCHECK)及MolProbity(Ramachandran plot)验证；UCSF Chimera进行结构叠合比对；AutoDock 4将建模酶与透明质酸(hyaluronic acid, HA)四糖/六糖及硫酸软骨素A(chondroitin sulfate A, CSA)四糖/六糖进行分子对接(molecular docking)；粗毒以Bradford法定蛋白浓度，浊度法(turbidimetric assay)于400 nm检测HA及CS降解吸光度变化评估酶活及37℃、45℃、60℃、72℃下热稳定性，以羊睾丸及意蜂(Apis mellifera)毒液透明质酸酶(hyaluronidase)为阳性对照。

三、研究结果

3.1 De novo transcriptome assembly and annotation（从头转录组组装与注释）

经质控Trimomatic过滤获44,370,955条双端150 bp清洁读段(clean reads)，Trinity de novo组装得206,153条contigs(>200 bp)，97,421条转录本(N50=1968 bp)，CD-HIT-EST(阈值0.95)聚类为93,303条unigenes。BUSCO评估Arachnida基因集完整度：原始组装92.4%完整BUSCOs(单拷贝8.6%+重复83.8%)，最终组装单拷贝74%+重复16.8%，缺失仅5.5%，CroCo检测无污染。注释比对动物毒素数据库鉴定26个毒成分类别，tapai神经网络验证177条选入转录本：126条离子通道毒素(ion channel toxins)——77条钠通道毒素(sodium channel toxins, NaTx)、46条钾通道毒素(potassium channel toxins, KTx)含与Mesobuthus eupeus λ-MeuTx 92%同源者、3条氯通道毒素(chloride channel toxins, ClTx)；27条酶类；24条其它(CRISPs等)。证实M. crucittii毒腺含典型Buthidae丰富离子通道调节肽及酶组分。

3.2 Phylogenetic analysis and in silico molecular characterization（系统发育与计算机分子特征分析）

转录组鉴定两条与东亚钳蝎(Olivierus martensii, 原Mesobuthus martensii)P86100透明质酸酶(hyaluronidase)96%一致的转录本， consensus序列GenBank登录号PP347720(mRNA)及WWA73958(蛋白)，录入ScorpDB。与71条蝎源Hyal多序列比对(multiple sequence alignment, MSA)，ML树显示Buthidae与Iurida两大单系分支(支持率94%)，M. crucittii Hyal归Buthid支内"Buthus"组进化枝。所有蝎Hyal具三保守基序：GDWW(Gx₃IDWExWRPxWx₃W)、FPDC(RPx₃WxYYxFPDCYxG)、GWGS(Gx₂Gx₃WGSS)；M. crucittii Hyal前体含28 aa信号肽，成熟肽382 aa，共12个Cys形成6对二硫键(C1?C6, C2?C5, C3?C4, C7?C9, C8?C10, C11?C12)，C端具EGF样域(EGF-like domain, EGF_3)含C7?C9/C8?C10/C11?C12三对二硫键。ProtParam预测成熟肽分子量47.5 kDa，理论pI 8.13，GRAVY ?0.424(亲水区)，不稳定指数37.58(<40稳定)，脂肪族指数(aliphatic index)79.46；Bepipred预测18个线性表位(linear epitopes)。

3.3 Molecular modeling（分子建模）

AlphaFold生成五模型取最高评分者(Ramachandran图94.3% favoured/99.7% allowed；ERRAT 95.756；Verify3D 87.27)。与O. martensii(P86100)及西方蜜蜂(Apis mellifera, PDB:1FCQ)晶体结构叠合显示Glyco_hydro_56域及EGF_3域高度保守；蝎Hyal含12 Cys(6二硫键)平分两域，蜜蜂Hyal仅4 Cys(2二硫键)对应其中子集。InterProScan/Simple Modular Architecture Research Tool(SMART)确认分Glyco_hydro_56催化域与EGF_3域；PTMGPT2预测并建模保留4个保守N-糖基化位点(N-glycosylation sites, Asn48/111/256/392)。催化位点选残基120?140区间(含UniProt标注催化Glu130/E130及GDWW基序)。

3.4 Molecular docking（分子对接）

以AutoDock 4将酶与HA四糖、HA六糖、CSA四糖、CSA六糖对接。结合能(binding energy, ΔG)：HA四糖?6.15 kcal/mol(最优，11条经典氢键，残基L85/E86/S87/A131/R133/Y138/W140/R142参与)>HA六糖?4.58 kcal/mol(9氢键)>CSA四糖?3.61 kcal/mol(6氢键，E130/R133/W140/R142)>CSA六糖+0.38 kcal/mol(不利作用，仅2氢键)。R133与R142在所有配体对接中均参与氢键/静电作用，为关键底物识别残基。HA四糖未直接结合E130但靠周边残基及空间排布获最稳定结合，提示四糖为适宜底物长度，CSA因电荷与空间不匹配结合较弱，HA六糖以上过长致结合质量略降。

3.5 Hyaluronidase activity measurement（透明质酸酶活性测定）

粗毒蛋白浓度56.83 μg/mL。浊度法显示M. crucittii粗毒15 min及30 min均可显著降解hyaluronan(HA)，活性与蜜蜂毒及羊睾丸阳性对照相当；延长孵育(0.5?120 h)证实可降解chondroitin sulfate(CS)但速率明显低于HA降解，提示双底物(bifunctional/substrate-promiscuous)活性。热稳定性实验显示37℃、45℃、60℃、72℃预处理后相对酶活无显著下降，表明该Hyal于全毒液中具高热稳定性(thermostable)。

四、讨论与结论翻译

讨论归纳如下：M. crucittii毒腺转录组丰富度与近缘Buthidae相符，首次以RNA-seq解析Mesobuthus属毒组成。Hyal序列系统发育追踪蝎演化史，Buthida与Iurida分开，支持作分子诊断特征(molecular diagnostic character)。in silico预测47.5 kDa成熟肽(实验因N-糖基化N-glycosylation常现55?60 kDa)具GH56家族特征、三保守基序(GDWW/FPDC/GWGS)及六二硫键增强热稳。