《ACS Omega》：Molecular Mimicry Inference, Therapeutic Target Identification, and Natural Product Inhibition Study in Cardiobacterium valvarum through Bioinformatics Approach
编辑推荐：
Cardiobacterium valvarum是感染性心内膜炎（infective endocarditis, IE）的已知病因。自身免疫反应被认为参与了IE中观察到的炎症过程。分子模拟显著影响病原体与宿主免疫系统之间的相互作用，常导致自身免疫反应。本研究采
Cardiobacterium valvarum是感染性心内膜炎（infective endocarditis, IE）的已知病因。自身免疫反应被认为参与了IE中观察到的炎症过程。分子模拟显著影响病原体与宿主免疫系统之间的相互作用，常导致自身免疫反应。本研究采用生物信息学方法，识别C. valvarum中可能引发免疫反应并促进炎症的具有分子模拟特征的候选蛋白。基于同源性和抗原性，共有18种蛋白显示出自身免疫潜力。进一步应用减除蛋白质组学来识别病原体中的潜在药物靶点，并利用purK基因产物（5-(carboxyamino)imidazole ribonucleotide synthase）进行下游分析。采用分子对接评估来自Ocimum tenuiflorum、Leonurus cardiaca和Panax ginseng的天然产物的抑制潜力，这些植物先前已被报道对心内膜炎具有保护作用。根据AutoDock Vina的结合评分，优先选出5种化合物：(1) 来自P. ginseng的一种复杂三萜类化合物，(2) 来自L. cardiaca的一种黄酮类衍生物，(3) 来自O. tenuiflorum的一种甾体化合物和一种黄酮类化合物。其中，LTS0201798和LTS0158828似乎更适合口服给药，而使用纳米载体或与环糊精复合也可提高其他化合物的适用性。LTS0201798和LTS0158828在100 ns分子动力学模拟中显示出良好的结合亲和力和稳定的复合物形成。因此，研究人员提出未来研究应专注于通过实验测定验证这些发现，以确认这些化合物在临床环境中的疗效。
**研究背景**
Cardiobacterium valvarum是一种革兰阴性菌，属于人类鼻咽和口腔正常菌群，主要引起感染性心内膜炎（infective endocarditis, IE）。IE中自身免疫反应被认为参与炎症过程，而分子模拟（即病原体抗原与宿主分子结构相似导致交叉反应）在自身免疫病发生中起关键作用。然而，C. valvarum在IE中的分子模拟机制及针对性治疗靶点尚未充分研究。为此，研究人员开展本研究，旨在通过生物信息学方法评估C. valvarum蛋白与人蛋白的分子模拟潜力，并利用计算机辅助药物设计（computer-aided drug design, CADD）识别针对非同源蛋白的天然产物抑制剂，为IE治疗提供新靶点和候选化合物。论文发表在《ACS Omega》。
**主要关键技术方法**
研究人员从NCBI和UniProt获取C. valvarum（菌株NCTC13294）和人类蛋白质组数据。通过BLAST比对识别同源序列（e值=0，序列一致性≥50%，模拟片段≥10个氨基酸），用VaxiJen评估抗原性，利用TM-align进行结构叠加（RMSD<2 ?）。通过ADEx数据库关联自身免疫病。在药物靶点挖掘中，用CD-HIT去除冗余（阈值90%），与必需基因数据库DEG和CEG比对（e值≤1×10
-5，bit-score≥100），排除人类同源物和肠道菌群同源物，用SPIDER预测药物性（阈值≥0.95），CELLO预测亚细胞定位，进一步根据分子大小（200~500氨基酸）和功能类别筛选，最终选定purK酶（5-(carboxyamino)imidazole ribonucleotide synthase）作为靶点。从LOTUS数据库获取来自Ocimum tenuiflorum、Leonurus cardiaca和Panax ginseng的化合物，用AutoDock Vina进行分子对接。利用pkCSM、SWISS-ADME和ProtoX III评估ADMET性质，用CALCES平台预测环糊精复合自由能。用GROMACS进行100 ns分子动力学模拟，分析RMSD、RMSF、半径回旋半径（radius of gyration, Rg）、氢键以及MM/GBSA和吉布斯自由能。
**研究结果**
- **分子模拟（Molecular Mimicry）**：通过BLAST和VaxiJen分析，29个同源蛋白中有20个至少含一个模拟序列，18个至少含一个抗原性模拟序列。DnaK和ATP合酶β亚基模拟数量最多（各8个），胸苷酸合酶等仅含1个模拟。结构叠加显示抗原性模拟蛋白的RMSD<2 ?。ADEx数据库关联显示，TYMS、GSR、PGD、TALDO1、MDH、GALE1、ATP5F1A在至少一种自身免疫病（系统性红斑狼疮SLE、类风湿关节炎RA、干燥综合征SjS、系统性硬化症SSc、1型糖尿病T1D）中上调，部分基因（如ENO2、HSPD1）下调。
- **药物靶点挖掘（Drug Target Mining）**：初始蛋白质组2595个蛋白经CD-HIT去冗余后得到2509个唯一序列。与DEG和CEG比对分别获得969和888个蛋白，858个共有。排除人类同源物后剩余403个。经SPIDER药物性筛选（≥0.95）保留25个，排除膜相关或周质蛋白后剩余14个，再经大小（200~500氨基酸）和功能类别过滤，最终获得8个胞质蛋白。基于药物性、结构可处理性、功能相关性及新颖性，选定purK酶（属于连接酶家族，参与嘌呤从头合成）作为最终靶点。
- **对接（Docking）**：对三种植物来源的化合物进行分子对接，优先选出5个化合物。LTS0177577（来自P. ginseng的三萜类）结合亲和力最高（-10.7 kcal/mol），其次为LTS0032845（-10.3 kcal/mol），其余三个（LTS0201798、LTS0158828、LTS0015495）均为-10.0 kcal/mol。共同参与相互作用的残基包括Glu257、Ile69、Tyr144、Gly143、Gln10等。
- **ADMET**：化合物分子量576.9~783.0 g/mol，LogP值-1.69~5.85，水溶性低至中等。LTS0032845和LTS0015495具较高Caco-2渗透性和肠吸收。所有化合物均为P-糖蛋白（P-gp）底物。LTS0158828有正分布容积（VDss），组织分布良好；LTS0032845和LTS0015495血脑屏障渗透性较高。LTS0158828为AMES阳性（致突变性）。LTS0201798和LTS0158828符合Lipinski规则，适合口服；其它化合物因分子量大或氢键供体/受体多而有口服生物利用度挑战。Protox III预测显示LTS0201798和LTS0158828毒性类别6（低毒），但存在心脏毒性风险（hERG II抑制）。
- **制剂增强（Formulation Enhancement）**：环糊精（CD）复合计算显示，LTS0201798与γ-CD复合自由能最负（-19.087 kJ/mol），LTS0158828与多种CD均显示有利复合。LTS0177577与DMCD、TMCD、RMCD复合良好。LTS0015495和LTS0032845复合自由能较差。
- **分子动力学模拟（MD Simulation）**：100 ns模拟显示，apo形式purK酶的平均RMSD高于3 ?，而5个化合物结合后RMSD均低于2 ?，其中LTS0158828、LTS0201798和LTS0177577复合物RMSD低于1.5 ?。RMSF分析显示配体结合降低了残基柔韧性，尤其Asp145等波动减小。Rg在20~21 ?之间稳定。MM/GBSA计算表明LTS0158828总结合能最低（-46.03 kcal/mol），LTS0015495最高（-19.79 kcal/mol）。吉布斯自由能景观显示LTS0201798复合物具有最低能量最小值，稳定性最高。
**讨论与结论**
讨论部分指出，C. valvarum蛋白与人蛋白的结构相似性可能导致交叉反应抗体产生，损伤心脏组织。16个模拟相关基因在人胎儿心脏中表达（除ATP5F1A和ATP5F1B），支持分子模拟在IE中的潜在作用。药物设计聚焦于非同源蛋白purK酶，通过CADD筛选出5个候选化合物。LTS0158828在MM/GBSA中结合最稳定，LTS0201798在吉布斯自由能中构象最稳定。但部分化合物存在hERG II抑制、AMES阳性等风险，需通过结构修饰、纳米制剂或环糊精复合优化。研究结论翻译如下：C. valvarum的分子模拟计算推断揭示了可能参与IE自身免疫反应的蛋白和通路。多个含有抗原性模拟序列的人同源蛋白通过ADEx数据库关联到免疫介导疾病（包括干燥综合征SjS、类风湿关节炎RA、系统性红斑狼疮SLE、1型糖尿病T1D）。尽管这些关联不能确立因果关系，但突出了感染期间易受交叉反应影响的共享免疫相关通路。通过精心选择最小化与人蛋白模拟的靶点，可能开发出有效治疗感染而不引发或加剧自身免疫反应的疗法。据此，purK酶因其与人类蛋白无同源性而被确定为有前景的靶点。经文献综述筛选出三种对IE有保护作用的植物，其代谢物与purK酶对接后基于结合亲和力优先选出5个化合物，它们具有进一步优化以改善ADMET性质和降低毒性的潜力。环糊精复合可提高溶解度。为验证生物信息学分析结果，需要进一步开展体外和体内研究，重点阐明观察到的毒性机制并制定降低风险、最大化抗C. valvarum疗效的策略。
