《Frontiers in Bioinformatics》：An integrated subtractive genomics and immunoinformatics approach for designing a universal multi-epitope vaccine against Brucella spp.
编辑推荐：
摘要
引言：布鲁氏菌属（Brucella spp.）是革兰阴性菌，可导致免疫功能低下个体和家畜的布鲁氏菌病。由于慢生长潜伏表型，当前抗生素不足以治疗感染。缺乏批准用于人类针对该病原体的疫苗构成重大公共卫生关切，并表明迫切需要新型预防性干预措施。
方法：
摘要
引言：布鲁氏菌属（Brucella spp.）是革兰阴性菌，可导致免疫功能低下个体和家畜的布鲁氏菌病。由于慢生长潜伏表型，当前抗生素不足以治疗感染。缺乏批准用于人类针对该病原体的疫苗构成重大公共卫生关切，并表明迫切需要新型预防性干预措施。
方法：在本研究中，研究人员将反向疫苗学（reverse vaccinology）方法与泛基因组分析（pan-genome analysis）相结合，以识别潜在疫苗靶标。研究人员筛选了蛋白质的抗原性（antigenicity）、可溶性（solubility）、免疫原性（immunogenicity）和亚细胞定位（subcellular localization）。鉴定了具有高免疫原性和可溶性的B细胞和T细胞表位。研究人员评估了多表位疫苗构建体，并进一步根据其理化性质进行了分析。进行了分子对接（molecular docking）、构象动力学（conformational dynamics）、计算机克隆（in silico cloning）和免疫模拟（immune simulations），以确定最佳疫苗候选物。
结果：四个蛋白——触发因子（trigger factor）、外膜蛋白组装因子BamA、脲酶亚基β（UreB）和脲酶亚基α（UreC1）——被认为是潜在的疫苗靶标。共筛选出26个B细胞和97个T细胞表位，具有显著的免疫原性和可溶性。生成了十二个多表位疫苗构建体，其中Vc7基于结构和理化性质被选为最佳。分子对接分析揭示了与2FSE和2Z65的良好相关性，进一步分析显示Vc7对2FSE表现出更强的结合亲和力（-135.24 kcal/mol），由疏水接触、盐桥和分子间氢键介导，使其成为理想的疫苗复合物，并通过150 ns分子动力学模拟（molecular dynamics simulation）进行了验证。计算机克隆建立了构建体的兼容性，免疫模拟确认了Vc7引发T细胞、B细胞、抗体和细胞因子介导反应的潜力。
结论：Vc7被鉴定为结构稳定且高度免疫原性的构建体，提示其作为预防布鲁氏菌病的通用多表位疫苗候选物的潜力。
论文解读：整合差异基因组学与免疫信息学设计布鲁氏菌通用多表位疫苗
**研究背景与问题**
布鲁氏菌属（*Brucella* spp.）是革兰阴性、兼性胞内寄生菌，导致人兽共患布鲁氏菌病（brucellosis），表现为发热、关节痛、疲劳等症状，严重者可转为慢性。该病在中亚、中东、地中海地区、非洲和印度等地流行，世界卫生组织将其列为最被忽视的人兽共患病之一，每年约160–210万新发病例。目前人类尚无获批疫苗，标准抗生素疗法（如多西环素联合利福平或链霉素）存在疗效有限、复发率高（可达15%）、链霉素需肌注不便等缺陷。布鲁氏菌的慢生长、潜伏态和胞内生存能力导致治疗失败和慢性化。因此，亟需开发安全有效的预防性疫苗。
**研究内容与结论**
研究人员提出将泛基因组分析（pan-genome analysis）与反向疫苗学（reverse vaccinology）以及免疫信息学（immunoinformatics）整合的策略，以设计针对布鲁氏菌属的通用多表位疫苗。通过分析185株完整布鲁氏菌基因组，鉴定核心蛋白，经抗原性、可溶性、非人源同源性筛选，最终选出四个保守蛋白（Tig、BamA、UreB、UreC1）作为疫苗靶标。预测并筛选B细胞和T细胞表位，构建12个多表位疫苗构建体，最终选择Vc7为最优候选。分子对接显示Vc7与HLA等位基因（2FSE和2Z65）具有强结合亲和力，150 ns分子动力学模拟、主成分分析和自由能景观分析证实复合物稳定性。免疫模拟预测Vc7能诱导IgM、IgG1、IgG2抗体、细胞因子（IFN-γ、IL-2）及T/B细胞记忆反应。计算机克隆显示Vc7可在*E. coli* K-12中高效表达。该研究发表在《Frontiers in Bioinformatics》。
**关键技术与方法**
研究人员采用了以下主要技术方法：① 泛基因组分析：利用Roary v3.13.0对185个完整布鲁氏菌基因组进行基因簇分析，鉴定核心基因；② 反向疫苗学筛选：通过BLASTp排除人类同源蛋白，用Gneg-Ploc、CELLO、PSORTb预测亚细胞定位，用VaxiJen v2.0、ANTIGENpro评估抗原性，用SOLpro评估可溶性；③ 表位预测：采用ProPred-I和IEDB工具预测MHC-I和MHC-II结合表位，用BepiPred v3.0和ABCPred预测线性B细胞表位，用ToxinPred和VaxiJen筛选无毒、高抗原性表位；④ 疫苗构建与优化：使用L7/L12核糖体蛋白、HBHA等佐剂和PADRE序列及不同接头（EAAAK、GPGPG、RVRR）组装多表位构建体，用trRosetta进行3D建模，MolProbity等验证结构；⑤ 分子对接与动力学：用ClusPro v2.0对接至HLA等位基因和TLR4，用GROMACS进行150 ns分子动力学模拟，iMODS进行正态模分析；⑥ 免疫模拟与克隆：用C-IMMSim预测免疫反应，用JCat进行密码子优化，虚拟克隆至pET-28a(+)载体。所有基因组数据来源于NCBI数据库。
**研究结果**
3.1 泛基因组分析与潜在疫苗候选鉴定
通过分析185株布鲁氏菌基因组，鉴定出1963个核心基因（共7759个基因簇）。经抗原性、可溶性、非人类同源筛选，最终选出四个蛋白：Tig（触发因子）、BamA（外膜蛋白组装因子）、UreB（脲酶亚基β）和UreC1（脲酶亚基α）。这些蛋白与细菌生存和宿主适应相关。
3.2 免疫原性表位分析
3.2.1 MHC-I类T细胞表位：预测获46个高免疫原性、无毒、高抗原性表位。
3.2.2 MHC-II类T细胞表位：预测获51个非毒性、高抗原性表位。
3.2.3 HLA等位基因功能聚类分析：热图展示了表位与等位基因的相互作用。
3.2.4 全球HLA群体覆盖分析：综合MHC-I和MHC-II表位覆盖率达约96.27%，覆盖多个人群。
3.2.5 线性B细胞表位：通过BepiPred和ABCPred共识预测获得26个重叠表位。
3.2.6 表位保守性分析：所选表位在多种布鲁氏菌种中高度保守（80%序列同一性阈值）。
3.3 疫苗构建体设计
最终构建体Vc7包含上述表位、L7/L12核糖体蛋白佐剂、PADRE序列、EAAAK/GPGPG/RVRR接头和C端His标签，共176个氨基酸。
3.4 疫苗构建体的结构化学属性评估
12个构建体中，Vc7表现出非过敏原、高抗原性、可溶性、负GRAVY值、稳定性指数27.36、脂肪族指数57.84等优良理化特性。
3.5 计算建模与结构验证
Vc7经trRosetta建模和GalaxyRefine优化后，Ramachandran图显示95.98%残基在允许区，Verify 3D通过，结构质量高。
3.6 构象B细胞表位
ElliPro预测Vc7表面有5个不连续B细胞表位，预测评分0.506–0.709。
3.7 二硫键工程
鉴定一对Ala-Gly残基可突变为半胱氨酸形成二硫键以增强稳定性。
3.8 计算机对接分析
Vc7与2FSE（HLA-DR1相关等位基因）结合能最高（-135.24 kcal/mol），通过盐桥和氢键互作；与2Z65结合能为-96.18 kcal/mol。表位残基积极参与受体结合。
3.9 正态模分析
iMODS分析显示Vc7-2FSE复合物特征值较大、相关运动良好、变形性低，表明动态稳定性优于Vc7-2Z65。
3.10 分子动力学模拟
150 ns模拟中，Vc7-2FSE和Vc7-2Z65复合物RMSD、RMSF、Rg、SASA、氢键数均显示结构稳定性。
3.11 主成分分析与自由能景观
PCA和FEL显示复合物达到低能稳定构象，Vc7-2FSE在相空间中聚集更紧密。
3.12 动态交叉相关分析
DCCM显示正相关和负相关运动，Vc7-2FSE中相关运动更协调。
3.13 密码子优化与计算机克隆
Vc7优化后CAI=1，GC含量55.30%，成功虚拟克隆至pET-28a(+)载体。
3.14 免疫模拟
三次免疫后，IgM、IgG1、IgG2滴度升高，IFN-γ和IL-2水平上升，B细胞和T细胞（包括记忆细胞）群体增加。
**讨论与结论**
研究人员在讨论中强调，本研究克服了以往研究局限于单菌株抗原选择或仅进行抗原识别的不足，通过整合泛基因组与反向疫苗学实现了跨种保守疫苗靶标的系统性筛选，并完成了从表位预测到疫苗构建体结构验证的全流程。Vc7候选疫苗在计算机水平上表现出良好的稳定性、受体结合能力和免疫潜力。最终结论指出：Vc7是有前景的布鲁氏菌属通用多表位疫苗候选物，但需要进一步实验验证其治疗潜力。
