副溶血性弧菌脂蛋白定位系统转运蛋白结构与功能的生物信息学分析

副溶血性弧菌(Vibrio parahaemolyticus)是一种重要的人畜共患病原菌[1]，广泛分布于水产食品及水产养殖环境中。该菌不仅能够引发腹泻、败血症等严重的食源性疾病[2]，对人类健康及食品安全造成了极大的威胁，而且能够感染南美白对虾的幼虾，引起对虾急性肝胰腺坏死病(Acute hepatopancreatic necrosis disease，AHPND)[3]，造成了严重的经济损失，制约着水产行业的可持续发展[3]。脂蛋白是副溶血性弧菌细胞外膜的重要组成部分，与其致病机制密切相关，在该菌细胞膜合成、黏附与侵袭、耐药性产生等重要生理过程中发挥着至关重要的作用[4]。

脂蛋白定位系统(Localization of lipoprotein system，Lol)由LolA–E五个蛋白构成[5]，主要负责副溶血性弧菌细胞内脂蛋白的转运与定位。现有研究表明，Lol系统通过细胞内膜上的ABC转运体LolCD2E复合物，将脂蛋白转运至周质空间，由周质空间内的伴侣蛋白LolA将其运送至细胞外膜，再通过外膜受体蛋白LolB将其接收并正确定位于细胞外膜内侧[6-8]，从而使脂蛋白发挥正常的生理功能。然而，现阶段Lol系统的三级结构还未完全明晰，且大部分研究常基于大肠杆菌[4, 9]、铜绿假单胞菌[10-11]、沙门氏菌[12]等细菌，对于副溶血性弧菌Lol的系统研究尚未展开，大大限制了人们对该菌致病及耐药机制的理解。

因此，本研究通过生物信息学技术，首次分析了副溶血性弧菌Lol系统转运蛋白的理化性质，并对其信号肽及跨膜区进行了预测，进一步基于STRING数据库探究了Lol系统相关蛋白之间的相互作用关系，并通过同源建模技术预测这些蛋白的三级结构，揭示副溶血性弧菌脂蛋白转运过程中关键结构，从生物信息学的角度展示副溶血性弧菌Lol系统转运脂蛋白的核心步骤。本研究可为后续副溶血性弧菌Lol系统转运蛋白的表达纯化、结构解析提供重要的数据基础，为副溶血性弧菌Lol系统功能的进一步揭示提供理论依据。

在NCBI的Protein数据库中，选取副溶血性弧菌菌株RIMD 2210633为研究对象，该菌株为世界上第一株被测序的副溶血性弧菌菌株，具有完整的基因组信息，广泛应用于揭示副溶血性弧菌的致病机理，被视为副溶血性弧菌相关研究中的模式菌株[13]。下载该菌脂蛋白转运系统Lol相关转运蛋白LolA、LolB、LolC、LolD以及LolE的氨基酸序列，并拼接LolCD2E蛋白复合物的氨基酸序列，用于后续的生物信息学分析。

使用ExPASy在线软件(http://web.expasy.org/ protparam/)分析了副溶血性弧菌Lol系统转运蛋白LolA、LolB、LolC、LolD、LolE和LolCD2E蛋白复合物的分子结构式、相对分子质量、平均亲水性、理论等电点、摩尔吸光系数等一系列理化性质[14]。

使用SignalP 4.0 Server在线软件(http://www. cbs.dtu.dk/services/SignalP-4.0/)对副溶血性弧菌Lol系统转运蛋白的信号肽进行预测分析[15]。

使用TMHMM-2.0 (http://www.cbs.dtu.dk/ services/TMHMM-2.0/)对副溶血性弧菌Lol系统转运蛋白是否含有跨膜区进行预测，并对跨膜区的个数及所在位置进行分析[16]。

基于副溶血性弧菌Lol系统转运蛋白的氨基酸序列，通过STRING数据库(https://string-db.org/)绘制了蛋白质互作网络图(Protein Protein Interaction Network，PPI)，分析了副溶血性弧菌LolA-E蛋白之间的相互作用关系，以及Lol系统相关蛋白之间的互作关系。通过以上蛋白的注释、结构及结构域、基因共表达、邻接、共生等信息，分析副溶血性弧菌Lol系统转运蛋白的功能[17]。

进一步运用SWISS-MODEL在线数据库(https://swissmodel.expasy.org/)对副溶血性弧菌Lol系统转运蛋白进行同源建模分析，预测并构建了LolA-E的三级结构[18]，并着重分析了副溶血性弧菌LolA与LolC之间发挥脂蛋白转运功能的关键的Hook结构[6]。

利用ExPASy软件分析了副溶血性弧菌Lol系统转运蛋白LolA、LolB、LolC、LolD、LolE以及LolCD2E复合物的理化性质，结果如表 1所示。LolA蛋白由208个氨基酸构成，分子式为C1055H1621N275O321S3，相对分子质量为23389.31 Da。LolB蛋白由212个氨基酸构成，分子式为C1084H1709N285O319S6，相对分子质量为24030.57 Da。LolC蛋白由405个氨基酸构成，分子式为C1966H3155N521O555S22，相对分子质量为43676.16 Da。LolD蛋白由235个氨基酸构成，分子式为C1121H1834N322O349S9，相对分子质量为25695.38 Da。LolE蛋白由414个氨基酸构成，分子式为C2042H3276N540O567S13，相对分子质量为44880.53 Da。LolCD2E蛋白由1289个氨基酸构成，分子式为C6250H10093N1705O1817S53，相对分子质量为139893.40 Da。

进一步分析显示：LolA、LolB和LolD蛋白呈现酸性，理论等电点(pI值)分别为5.19、6.19和5.71，LolCD2E蛋白复合物偏中性，pI值为6.90，LolC和LolE蛋白为碱性蛋白，pI值分别为8.89和9.40。LolA–E及LolCD2E的摩尔吸光系数分别为：1.195、1.683、0.777、0.174、1.011和0.633。明晰这些蛋白的等电点和摩尔吸光系数，有助于后续这些蛋白提取与纯化方法的构建与优化[19]。

此外，LolA、LolB和LolD蛋白平均亲水性 < 0，为亲水性蛋白，其余为疏水性蛋白。疏水性蛋白有可能含有跨膜结构[20]，因此，本研究进一步分析了Lol系统转运蛋白的信号肽和跨膜结构。

号肽通常位于分泌蛋白和膜蛋白氨基酸序列的N端，由16到30个氨基酸残基组成[21]。信号肽的主要作用是促进蛋白分泌到细胞外，信号肽的有无直接决定着重组蛋白表达与纯化方式的选择[22]。因此，本研究通过在线数据库SignalP 4.0 Server，根据默认的神经网络模型(Neural networks)和隐马尔可夫模型(Hidden Markov models)对副溶血性弧菌Lol转运系统转运蛋白进行信号肽分析。结果如图 1所示，LolA和LolB蛋白含有信号肽切割位点，信号肽分值(S-score)的峰值分别位于第20位氨基酸和第25位氨基酸处。而LolC、LolD、LolE及LolCD2E均无显著的信号肽切割位点。

使用在线分析软件TMHMM-2.0对Lol系统转运蛋白的跨膜方式进行分析，结果如图 2所示，LolA、LolB和LolD蛋白结构中无明显的跨膜区域。LolC蛋白含有5个跨膜区，其跨膜螺旋区(TMhelix)预测位置分别在29–51、272–294、314–336、343–365、369–391的氨基酸处。LolE蛋白含有4个跨膜区，其跨膜螺旋区(TMhelix)预测位置分别在23–45、267–289、318–340、375–397的氨基酸处；LolCD2E蛋白含有9个跨膜区，其跨膜螺旋区(TMhelix)预测位置分别在29–51、272–294、314–336、343–365、369–391、898–920、1142–1164、1193–1215和1250–1272的氨基酸处。

利用STRING数据库，对副溶血性弧菌Lol系统相关蛋白的相互作用网络进行分析，结果如图 3所示。每个圆形的节点代表一种蛋白，节点之间的连线表示两个蛋白之间存在相互作用关系，不同颜色的连线表示不同的相互作用途径，两种蛋白之间的连线越丰富，表示这两种蛋白间具有更为密切的相互作用关系[23]。对于Lol系统转运蛋白之间的相互作用分析可知，LolA–E五个蛋白的编码基因均共表达，负责脂蛋白的合成与转运。LolC、LolE和LolD位于副溶血性弧菌的细胞内膜上，共同组成LolCD2E蛋白复合物，其中LolD是脂蛋白转运系统ATP结合蛋白，通过消耗ATP的方式，将脂蛋白运输通过细胞内膜。LolA为伴侣蛋白，参与脂蛋白从细胞内膜到细胞外膜的转运。LolB是一种外膜脂蛋白，在脂蛋白被LolA蛋白释放后与外膜结合的过程中起关键作用。此外，STRING蛋白相互作用网络显示，Lol系统转运蛋白LolA-E与BamA、Pal、MacB、CmeC等蛋白均具有密切的互作关系。

利用SWISS-MODEL在线数据库对副溶血性弧菌Lol系统转运蛋白进行同源建模，结果显示，副溶血性弧菌LolA蛋白与其同源蛋白大肠杆菌(Escherichia coli) LolA[6]相似度为43%，LolB与E. coli外膜脂蛋白受体LolB[24]序列相似度为35%。LolC、LolD与LolE均依据放线共生放线杆菌(Aggregatibacter actinomycetemcomitans) ATP结合蛋白MacB结构进行同源建模[25]，相似度分别为29%、40%和29%。选取以上蛋白的结构作为模板，预测了副溶血性弧菌LolA–E三级结构预测模型，并结合上文对Lol系统的功能分析，模拟了副溶血性弧菌Lol系统转运蛋白在细胞膜上的位置分布，及运输脂蛋白的过程(图 4)。

在大肠杆菌、霍乱弧菌、沙门氏菌、鼠疫耶尔森菌、流感嗜血杆菌、铜绿假单胞菌等细菌的Lol系统中，LolC和LolE蛋白中存在一段特殊的Hook结构，是LolCD2E复合物与LolA结合并转运脂蛋白的关键区域[6]，但尚未有研究揭示副溶血性弧菌是否存在类似的结构。本研究以其他革兰氏阴性菌Hook结构为基础，比对了副溶血性弧菌与其他革兰氏阴性菌LolC和LolE蛋白的氨基酸序列，结果如图 5-A所示，副溶血性弧菌LolC和LolE蛋白中均存在一段保守的Hook结构，位于2个β折叠的序列之间，LolC的Hook结构在第172位与第182位氨基酸之间，LolE的Hook结构在第172位与第181位氨基酸之间。并进一步结合副溶血性弧菌LolC与LolE蛋白的同源建模结果，绘制了Hook区域的三维结构，结果如图 5-B所示，该Hook区域呈现特有的β-发夹结构，使LolCD2E具有结合并转运脂蛋白的特殊功能。

本研究首次针对副溶血性弧菌Lol系统转运蛋白相关的氨基酸序列进行了系统的生物信息学分析。首先分析了副溶血性弧菌Lol系统转运蛋白的基本性质，LolA、LolB和LolD为酸性亲水蛋白，无跨膜区域，研究显示LolB是位于外膜上的脂蛋白受体蛋白，因此LolB可能是通过特定锚定方式与外膜结合[26-27]。LolC和LolE为碱性疏水蛋白，包含4–5个跨膜区，LolCD2E为中性疏水蛋白复合物，包含9个跨膜区。LolA和LolB蛋白含有信号肽切割位点，LolC、LolD、LolE及LolCD2E均无显著的信号肽切割位点。因此，在后续该菌LolA–E、LolCD2E的重组表达载体设计中，应根据实验的需要添加或删除信号肽区域，以保证蛋白的正常表达。并针对LolC、LolE、LolCD2E这3个具有跨膜区域的膜蛋白，结合去垢剂等物质进行处理，以确保此类蛋白正确地提取与纯化。

进一步通过蛋白质-蛋白质互作网络分析阐述了副溶血性弧菌Lol系统相关蛋白互作关系，为后续副溶血性弧菌Lol系统转运蛋白功能的研究提供理论参考。通过分析STRING数据库提供注释信息可知，LolA–E五个蛋白的编码基因均共表达，负责脂蛋白的合成与转运。LolC、LolE和LolD共同组成LolCD2E蛋白复合物，位于副溶血性弧菌的细胞内膜上，拥有类似于MacB的周质核心区域结构。此类结构的存在，赋予了副溶血性弧菌Lol系统消耗ATP运输脂蛋白的重要功能。LolA为周质空间内的伴侣蛋白，负责将脂蛋白从细胞内膜转运到细胞外膜上，再通过外膜受体蛋白LolB将其接收并正确定位于细胞外膜内侧。

此外，蛋白相互作用网络还显示，副溶血性弧菌Lol系统与BamA、Pal、MacB、CmeC等外膜蛋白均具有密切的互作关系。其中，BamA蛋白为外膜蛋白装配系统的重要组成部分[28]，Pal蛋白参与细菌外膜肽聚糖的定位[29]。MacB蛋白含有特定ABC转运体结构域，能够与Lol系统共同表达[30]。CmeC为重要外膜外排泵蛋白，含有OEP家族特定结构域，能够形成三聚体通道，允许革兰氏阴性菌输出多种底物[31]。由此可见，副溶血性弧菌的Lol系统不仅负责细胞脂蛋白的合成与转运，还在细胞外膜结构的组装、细菌生命活动的调控中发挥着重要的作用。

最后，本研究预测并分析了副溶血性弧菌Lol系统转运蛋白的三级结构，发现其具有与大肠杆菌相似的LolA和LolB结构，推测其可能拥有相似的脂蛋白转运机制[32]。副溶血性弧菌的LolC、LolD和LolE蛋白均含有MacB的同源结构，与Wang等[33]以及Khwaja等[34]的结论相符，也与上述STRING数据库注释相一致，进一步证明了LolCD2E可能和MacB具有类似的通过消耗ATP转运物质的机制[35]。此外，先前研究表明，大肠杆菌、沙门氏菌等革兰氏阴性菌的LolC和LolE蛋白中包含一段特殊的Hook结构，被视为LolCD2E蛋白复合物与LolA蛋白识别并转运脂蛋白过程中的关键结构[6]。本研究Lol系统转运蛋白结构分析显示，副溶血性弧菌LolC和LolE蛋白中也存在这样一段保守的Hook结构，揭示了副溶血性弧菌Lol系统转运脂蛋白的关键区域。若将该Hook结构进行破坏或阻断，可终止副溶血性弧菌脂蛋白的转运，从而抑制副溶血性弧菌的生长与繁殖。后续研究可针对副溶血性弧菌LolC和LolE蛋白中的Hook结构进行开展，进一步探究其脂蛋白转运功能，并以此为靶标研发新型抗菌类药物。