通过计算方法对与流行性病毒感染相关的MMP8、GZMK和OASL基因中高风险错义变体的结构与功能表征

为何微小的遗传变化会影响大规模疫情

当一种新病毒在群体中传播时，并非所有人都会出现相同程度的疾病。有些人能轻松抵御感染，而另一些人却会发展为危及生命的疾病。本研究提出了一个看似简单的问题：我们免疫系统基因中的小的遗传差异，是否能够解释这些差异？研究人员并非在实验室里直接检测，而是使用强大的计算模拟，探讨特定基因变体如何细微地重塑关键免疫蛋白，进而改变我们对诸如流感、埃博拉或冠状病毒等流行性病毒的应答。

三位免疫“助手”被放大检视

研究团队聚焦于三个人类基因，这些基因在早期的大规模分析中被标为对多种呼吸道病毒应答的核心参与者：MMP8、GZMK 和 OASL。每个基因都编码一种有助于控制我们如何抵抗感染的蛋白质。MMP8 有助于重塑受损组织并调节肺部及其他部位的炎症。GZMK 编码颗粒酶 K，这是一种由杀伤性免疫细胞释放的酶，用于攻击受病毒感染的细胞并影响炎症反应。OASL 产生一种能增强机体抗病毒警报系统并可直接干扰病毒复制的蛋白。在这些基因中，作者挑出了五个罕见但被评为高风险的“错义”变体——即单个碱基变化导致蛋白质中一个氨基酸被替换——因为多个预测工具一致认为它们可能有害。

模拟有害变体的影响

研究者没有在细胞或动物中测试这些变体，而是构建了正常与突变蛋白的三维数字模型，并对它们进行了一系列的计算测试。他们采用成熟的算法评估每一处替换是否被预测为会损害蛋白、会在多大程度上削弱整体稳定性，以及是否触及在物种间高度保守的区域——这是功能重要性的一个标志。随后，他们进行了长时间的分子动力学模拟——本质上是在原子尺度上的物理学“电影”——观察每个突变蛋白在数百纳秒内如何弯曲、压缩或松散，以及这些变化如何影响其与代表性小分子配体的相互作用。

不同基因，不同的结构命运

结果显示，并非所有高风险变体都会以相同方式出问题。在 MMP8 蛋白中，一个变体 D253N 导致结构收缩并减少对周围水的暴露，提示一种更为紧凑但不稳定、并能采样更广泛构象的状态。另一个 MMP8 变体 Y261S 则保持了总体折叠与尺寸接近正常，甚至表现出稍微更为刚性的特征。对于 GZMK，A42P 的影响较轻微，但 L122P 明显增加了局部柔性并改变了蛋白的紧密堆积。在 OASL 中，W216C 变体使结构松散、表面暴露增加并扰乱了远程的内部运动，这与一种放松、较不凝聚的折叠一致。主成分分析（将复杂运动压缩为少数主导模式）证实 D253N、L122P 和 W216C 尤其扩大了所探索构象的范围。

结合看似正常但行为已改变

为评估这些结构变化是否可能影响功能，研究组将每个正常与突变蛋白与选定小分子进行对接，然后通过进一步的模拟和能量计算对这些复合体进行精炼。所有变体仍然能够与其配体结合，有时结合强度的预测变化仅为适度。然而，结合的方式常常不同：蛋白与配体之间的接触模式发生了转移，在某些情况下——例如 MMP8 的 D253N——总结合自由能变得不那么有利，暗示相互作用可能更弱或不那么可靠。有趣的是，尽管 OASL 的 W216C 整体上使蛋白不稳定，但在计算中却显示对测试分子的结合更强，这表明增加的柔性有时能在与配体“抓握”方面带来改善，尽管它会削弱蛋白行为的其它方面。

对未来疫情的意义

对非专业读者来说，关键结论是：我们 DNA 中微小的遗传变化可以在不明显破坏蛋白功能的情况下，细微地重塑重要的免疫蛋白。这五个位于 MMP8、GZMK 和 OASL 的变体似乎都以扰动稳定性、运动与分子接触之间的平衡方式出现，这些改变可能会影响个体对病毒感染的反应——通过改变炎症、细胞杀伤或抗病毒信号传导途径。该研究并不能证明携带这些变体的人在真实疫情中会更好或更差；这需要细致的实验室和临床研究来验证。但通过指出哪些变体最显著地扭曲蛋白动力学，这项计算研究为实验优先级提供了清单，并强调了现代模拟工具在揭示下一次病毒暴发时为何有些人可能更脆弱方面的价值。

引用: Et-tanjaouy, M., Saih, A., Machich, O. et al. In silico structural and functional characterization of high-risk missense variants in MMP8, GZMK, and OASL genes associated with epidemic viral infections. Sci Rep 16, 12973 (2026). https://doi.org/10.1038/s41598-026-40467-w

关键词: 病毒感染易感性, 免疫基因变体, 蛋白质结构动力学, 计算免疫学, 错义突变