金黄色葡萄球菌中假设蛋白的功能与结构注释及分子对接研究的生物信息学方法

《BioMed Research International》：Bioinformatics Approaches for Functional and Structural Annotation and Molecular Docking Study of a Hypothetical Protein From Staphylococcus aureus

　　

金黄色葡萄球菌是一种常见的革兰氏阳性菌，定植于人类的皮肤和黏膜。它是医院感染的主要病原体，也是食品和医院中持续的污染源。金黄色葡萄球菌基因组的完整测序已经完成，然而其中编码的许多假设蛋白（hypothetical protein, HP）在功能上仍未表征。在本研究中，研究人员利用生物信息学方法对来自金黄色葡萄球菌HMSC77A05菌株的假设蛋白（登录号OHS92176.1）进行了全面的功能和结构分析。理化性质研究表明该蛋白稳定、可溶且定位于细胞质。利用NCBI保守结构域搜索（NCBI Conserved Domain Search）、InterProScan和Pfam数据库进行的功能注释表明，该蛋白是卤酸脱卤酶（haloacid dehalogenase, HAD）超家族中YqeG家族HAD IIIA型磷酸酶的成员。多重序列比对（multiple sequence alignment, MSA）显示许多金黄色葡萄球菌菌株在进化上非常相似，表明它们可能对其功能重要。PSIPRED和SOPMA均一致预测α-螺旋元素是蛋白质二级结构中最常见的部分。使用SWISS-MODEL和模板A0A2C6WKT6.1.A构建了三维模型，与靶蛋白序列同一性为71.43%。使用ERRAT、Verify3D、QMEAN和PROCHECK分析对模型质量进行了验证，所有分析均证实了其结构可靠性和立体化学准确性。STRING的蛋白质-蛋白质相互作用（protein–protein interaction, PPI）研究发现与核糖体、氨基酸生物合成和核苷酸分解代谢中涉及的酶存在关联。还发现活性位点有四个假定的催化口袋。随后使用六种植物化合物配体——槲皮素（quercetin）、山奈酚（kaempferol）、柚皮素（naringenin）、芹菜素（apigenin）、儿茶素（catechin）和姜黄素（curcumin）——进行了分子对接（molecular docking, MD）模拟，所有配体均表现出良好的结合亲和力。值得注意的是，山奈酚显示出最高的结合亲和力（ΔG = ?8.3 kcal/mol）。这些发现表明，这个以前未表征的假设蛋白可能作为针对金黄色葡萄球菌的抗菌剂的潜在治疗靶点。全面的计算分析为后续的实验验证和基于结构的药物开发提供了坚实的基础。

### 论文解读：金黄色葡萄球菌中假设蛋白的功能与结构注释及分子对接研究

#### 研究背景、现存问题与研究目的

金黄色葡萄球菌（*Staphylococcus aureus*）是一种常见的革兰氏阳性病原菌，可定植于人类皮肤和黏膜，引发从轻度皮肤感染到危及生命的全身性疾病（如心内膜炎、肺炎、败血症）等多种感染。该菌具有多种毒力因子，并能快速产生耐药性，尤其是耐甲氧西林金黄色葡萄球菌（methicillin-resistant *S. aureus*, MRSA）的出现，使得现有治疗方案失效，亟需开发新型抗菌药物。尽管金黄色葡萄球菌的基因组已完全测序，但其中约11%的编码蛋白被归类为“假设蛋白”（hypothetical proteins, HPs），即基于开放阅读框预测但缺乏实验表达或功能验证的蛋白。这些HP是未被探索的潜在生物分子储备，可能在细菌代谢、致病性及耐药性中发挥重要功能。其中，卤酸脱卤酶（haloacid dehalogenase, HAD）超家族包含多种参与磷酸代谢和调控的磷酸酶，是潜在的药物靶点，但金黄色葡萄球菌中许多HAD家族蛋白的功能仍不清楚。因此，研究人员旨在通过生物信息学方法，对来自金黄色葡萄球菌HMSC77A05菌株的一个特定HP（登录号OHS92176.1）进行功能与结构注释，并评估其作为抗菌药物靶点的潜力。该研究发表于《BioMed Research International》。

#### 主要关键技术方法

研究人员采用了多层次的计算策略，主要关键技术方法包括：（1）理化性质分析：使用ProtParam工具预测分子量、等电点、稳定性等。（2）功能域注释：利用NCBI保守结构域数据库（Conserved Domain Database, CDD）、Pfam和InterProScan识别蛋白家族与超家族。（3）同源建模与质量评估：通过SWISS-MODEL基于AlphaFold预测的模板A0A2C6WKT6.1.A构建三维结构，并用PROCHECK（拉氏图）、Verify3D、ERRAT、QMEAN和ProSA-Web进行验证。（4）系统发育分析：通过BLASTp搜索同源序列，使用Clustal Omega进行多重序列比对，并采用邻接法（neighbor-joining, NJ）构建系统发育树。（5）蛋白质-蛋白质相互作用网络分析：利用STRING数据库（版本12）预测互作蛋白。（6）活性位点预测与分子对接：使用PrankWeb识别结合口袋，并用AutoDock Vina（PyRx平台）对六个植物化合物（槲皮素、山奈酚、柚皮素、芹菜素、儿茶素、姜黄素）进行对接。样本来源为NCBI蛋白质数据库中的金黄色葡萄球菌HMSC77A05菌株序列。

#### 研究结果

**3.1 蛋白的理化特征**：该HP由175个氨基酸组成，分子量为20,171.45 Da，等电点（pI）为9.73（碱性），亲水性均值（GRAVY）为?0.298（亲水），不稳定指数（II）为29.61（稳定），表明该蛋白为稳定的、碱性的、亲水的胞质蛋白。

**3.2 胞内分布与可溶性**：CELLO v2.5预测该蛋白主要定位于细胞质（得分3.279），Protein-Sol分析显示其可溶性分数（QuerySol）约为0.7，显著高于群体平均可溶性（PopAvrSol）0.45，提示其在生理条件下可溶。

**3.3 功能域与家族/超家族分析**：NCBI-CDD、Pfam和InterProScan一致鉴定该蛋白属于HAD超家族中的YqeG家族HAD IIIA型磷酸酶，保守域覆盖残基6-168（E值1.48e-86），Pfam确认HAD样域（7-170，E值~1e-70），表明其可能水解5′-核苷酸磷酸基团。

**3.4 多重序列比对与系统发育树**：BLASTp结果显示该蛋白与多种金黄色葡萄球菌菌株的YqeG同源物有极高序列同一性（99.43%-100%，E值1e-123至3e-123）。MSA和系统发育树显示所有序列紧密聚类，反映进化保守性，暗示这些蛋白在细菌代谢中具有重要功能。

**3.5 二维结构分析**：SOPMA和PSIPRED预测二级结构以α-螺旋为主（SOPMA: 48%; PSIPRED: 40%），其次为随机卷曲、延伸链和β-转角，符合HAD超家族典型的α/β折叠特征。

**3.6 三维结构预测与质量评估**：基于模板A0A2C6WKT6.1.A（序列同一性71.43%，覆盖率93%，GMQE=0.93）构建三维模型。拉氏图显示94.9%残基位于最适区，5.1%在额外允许区，无残基在禁止区；ERRAT质量因子98.788；QMEAN4得分为0.72（|z-score|<1）；ProSA-Web计算模型z-score为?7.74（模板?7.44），均表明模型质量可靠，结构稳定性好。

**3.7 蛋白质-蛋白质相互作用分析**：STRING分析揭示该HP（对应登录号ABD30775.1）与多个重要蛋白高度互作（最高得分0.910），包括参与核糖体沉默（rsfS）、芳香族氨基酸生物合成（aroE）、NAD生物合成（nadD）、甲硫氨酸补救途径（mtnN）的蛋白以及多个保守HP，提示其参与关键的代谢和调控通路。

**3.8 活性位点检测**：PrankWeb识别出四个配体结合口袋，其中排名第一的口袋（Rank 1）得分为23.75，概率0.875，含19个残基，是主要活性位点。

**3.9 分子对接分析**：对六种植物化合物进行分子对接，所有配体均与蛋白活性位点结合，结合能范围为?7.9至?8.3 kcal/mol。山奈酚结合能最低（?8.3 kcal/mol），其次为芹菜素和柚皮素（?8.2 kcal/mol），槲皮素和儿茶素（?8.1 kcal/mol），姜黄素（?7.9 kcal/mol）。山奈酚通过三个氢键与Asp36、Ser71和Asn73相互作用，表明其作为先导化合物的潜力。

#### 讨论与结论

讨论部分指出，该研究通过集成理化分析、保守域鉴定、结构建模、进化分析、PPI作图及分子对接，为HP OHS92176.1构建了连贯的功能与结构框架。理化性质与YqeG家族成员相符；保守域分析强有力地将其归为HAD IIIA型磷酸酶；进化分析显示高度序列保守性（~99%同一性），表明其功能重要性；三维模型经多种方法验证可靠；PPI网络揭示其与NAD、氨基酸生物合成及核糖体调控等关键通路连通，进一步支持其生物学意义。分子对接结果中，山奈酚等植物化合物表现出优异结合力，提示该HP可作为抗菌药物开发的潜在靶点。但讨论强调计算研究存在局限性，需通过生化实验、定点突变及分子动力学模拟进行验证。结论部分翻译如下：在本研究中，研究人员使用多种生物信息学工具系统地表征了金黄色葡萄球菌编码的一个假设蛋白，提供了关于其结构身份、功能注释和治疗潜力的强计算依据。基于结构域的顺序分析表明，该蛋白属于HAD超家族IIIA型磷酸酶的YqeG家族，在磷酸转移和代谢调节反应中起重要作用。这一分类得到多种金黄色葡萄球菌菌株中YqeG同源物约99%序列同一性的进一步支持，表明显著的功能保守性。理化性质证实该蛋白具有稳定、可溶性胞内酶的特性。基于同源性的结构建模经多种立体化学质量评估方法验证，确认其可靠性和适用于下游计算机研究。PPI网络分析表明该蛋白与金黄色葡萄球菌生存和稳态所必需的关键细胞磷酸化和代谢调节过程功能相连。涉及多种抗菌剂的分子对接研究表明该蛋白可作为开发新型抗金黄色葡萄球菌治疗药物的可成药靶点。然而，预测的胞质定位限制了其作为表面暴露疫苗候选的应用。该研究结果为未来基于结构的药物发现工作提供了强大的计算框架，以开发金黄色葡萄球菌YqeG家族HAD磷酸酶作为新型治疗药物。本研究完全基于计算预测，因此建议通过体外和体内实验对蛋白的抑制敏感性和生物活性进行实验验证。这些发现促进了对细菌中功能尚未充分表征的蛋白质的后续计算和实验研究。