《Viruses》：Analysis of HA5 Hemagglutinins from Influenza A Viruses for Mass-Charge Patterns to Understand Antigenic Drift Related to Immune Escape
编辑推荐：
甲型流感病毒因其突变和逃避宿主免疫防御的能力而持续构成全球健康威胁，这一过程主要由血凝素（HA）的变化驱动。HA是流感病毒颗粒表面的一种糖蛋白，通过结合细胞受体在病毒进入宿主细胞过程中发挥关键作用。本研究利用预测建模算法，对2023年至2024年期间在中国大陆
甲型流感病毒因其突变和逃避宿主免疫防御的能力而持续构成全球健康威胁，这一过程主要由血凝素（HA）的变化驱动。HA是流感病毒颗粒表面的一种糖蛋白，通过结合细胞受体在病毒进入宿主细胞过程中发挥关键作用。本研究利用预测建模算法，对2023年至2024年期间在中国大陆从野鸟、海鸥、鸭、鹅、鸡以及人类中采集的多种流感病毒株的23条HA5蛋白序列进行了质量-电荷模式分析，以理解病毒进化的非线性趋势，特别是其抗原漂移和免疫逃逸潜力。研究人员识别并分析了HA5蛋白中的两个关键区域。第一个区域是“可突变区”（由质量-电荷模式分析鉴定的负重力中心），该区域频繁累积多种突变，导致抗原漂移和不稳定性，使病毒能够逃避已有的群体免疫。分析和监测质量-电荷变化可能有助于预测病毒的进化趋势以及在这些病毒采集区域未来变异株的出现，从而为所分析病毒及相关毒株的疫苗设计和开发提供信息。第二个区域是“稳定核心”（正重力中心），这是一个对病毒功能至关重要的高度保守区域。由于该区域在不同毒株间具有保守性，它可能代表开发通用疫苗和广谱抗体的有前景的靶点，这些疫苗和抗体能够同时针对所分析病毒及相关毒株提供针对多种变异株的保护。总之，研究人员通过预测建模算法对HA5蛋白序列进行质量-电荷模式分析，为病毒进化和免疫逃逸潜力评估提供了一种新的度量指标，并为针对所分析病毒及相关毒株设计疫苗提供了策略框架。
**1. 研究背景与问题**
流感是由流感病毒引起的急性呼吸道传染病。根据核蛋白和基质蛋白，流感病毒分为A、B、C、D四型；依据表面血凝素（HA）和神经氨酸酶（NA）糖蛋白，又可进一步分为多种亚型，如H1N1、H5N1等。过去一百年间，全球发生过四次主要流感大流行。在流感病毒中，A型流感病毒变异率最高，其表面抗原的漂移与转变直接决定流感的流行规模。抗原漂移可产生新毒株并引起中小型流行，而抗原转变则可产生全新亚型，易导致大规模甚至全球性大流行。因此，阐明A型流感病毒抗原漂移与免疫逃逸机制，对理解病毒进化和改进疫苗设计至关重要。现有研究虽已通过深度突变扫描、全基因组测序、结构生物学等手段揭示HA蛋白变异与免疫逃逸的部分特征，但常规分析方法常面临人工效率低或过度依赖计算资源等问题。本研究旨在利用一种基于氨基酸质量-电荷关系的半协方差系数分析工具，对2023—2024年中国大陆采集的A型流感病毒H5N1和H5N6亚型HA5蛋白序列进行分析，以揭示其非线性进化趋势、抗原漂移规律及免疫逃逸潜力。
**2. 研究内容与结论**
研究人员从GISAID数据库收集了23条HA5蛋白序列，包括19条禽源和4条人源序列，采样地点覆盖中国大陆不同地理区域及春、初夏、秋、冬等季节。通过Clustal Omega进行序列比对，将蛋白序列转换为电荷序列（+1、?1、0），利用Langlands公式计算质量-电荷值，并结合Pearson相关系数、半协方差系数和重力中心等指标分析序列间的线性与非线性关系。研究发现：禽类传播过程中HA5蛋白序列高度保守，正半协方差系数>0.994，负半协方差系数<0.004；而跨物种传播至人类时，序列出现明显分化，正半协方差系数降至0.977543，负半协方差系数升至0.014605，Pearson相关系数降至0.962938。正漂移重力中心集中在第303~305位氨基酸残基，构成高度保守的功能区域；负漂移重力中心广泛分布于第61~234位残基，形成突变热点区域，可能与抗原漂移和免疫逃逸相关。将漂移位点映射到HA5蛋白结构域后发现，正漂移中心位于茎部F′亚结构域，负漂移中心分布于伪酯酶（VE）和受体结合（RB）亚结构域。按月趋势预测显示，4—6月期间跨物种传播活动增强，正漂移中心移至第304~305位，负漂移中心与月份呈负相关（相关系数?0.43），且该时期VE亚结构域出现独特漂移模式。这些结果表明，半协方差系数分析结合质量-电荷漂移可视化可作为表征病毒蛋白进化与免疫逃逸机制的有效框架，为疫苗设计提供预测性信息。
**3. 主要关键技术方法**
本研究采用的主要方法包括：（1）序列比对：使用Clustal Omega对HA5蛋白序列进行多序列比对；（2）电荷序列转换：将氨基酸残基按电荷性质赋值为+1、?1或0，缺失片段以均值填充；（3）质量-电荷计算：采用Langlands公式整合氨基酸分子质量、疏水性、电荷和等电点信息，并引入反映病毒蛋白与宿主受体间Finsler距离的缩放因子λ（取值1、3、5、7、10），初步验证表明λ值不影响漂移模式，结果示例以λ=1计算；（4）半协方差系数分析：通过ReLU类非线性变换捕捉序列间非线性相关关系，区分正半协方差（第一、三象限）和负半协方差（第二、四象限）；（5）重力中心计算：分别计算正漂移（趋同）和负漂移（趋异）的重力中心，确定突变聚集区域。样本来源于GISAID数据库所收录的2023—2024年中国大陆不同地区的H5N1和H5N6毒株。
**4. 研究结果**
*4.1 质量-电荷漂移可视化与线性/非线性相关性分析*
以EPI4158283（A/Seagull/Hebei/qhd6/2024_HA）为初始序列，计算各HA5蛋白序列相对于初始序列的名义质量-电荷变化并绘制漂移图。结果显示，负漂移次数与宿主种类和纬度差相关：禽源宿主变幅为1~4次，人源宿主为5~10次，且随纬度差增大而增加。半协方差系数表明，禽类传播中HA5序列高度保守，而跨物种传播至人类时序列差异显著增大，Pearson相关系数降低，线性相关性减弱。该分析在HA5蛋白序列中识别出质量-电荷值呈增加或减少趋势的特定区域，提示局部突变可能影响病毒结合特性和抗原性。
*4.2 重力中心定位与功能结构域映射*
正漂移重力中心紧密聚集于第303~305位氨基酸残基，提示该区域及其邻近区域为HA5蛋白高度保守的功能结构域。负漂移重力中心广泛分布于第61~234位残基，形成显著突变热点。这些负漂移中心与人类流感病毒HA序列中经遗传学和质谱研究识别的突变热点重叠，间接验证了分析方法的可靠性。参照Zheng等人（2018）描绘的HA5亚型结构域图谱，所有分析序列的正漂移重力中心均位于F′亚结构域（茎部），而负漂移重力中心分布于VE和RB亚结构域（球状头部）。F′亚结构域结构保守且功能重要，可能是下一代流感疫苗和治疗性抗体的关键靶点；RB亚结构域是适应性突变热点，影响宿主特异性并促进人畜共患传播。
*4.3 按月进化趋势预测*
应用半协方差方法预测HA5蛋白的月度进化轨迹。所有禽源毒株的正漂移重力中心一致位于第303位残基，显示高度进化保守性。4—6月期间，跨物种传播至人类的活动增强，正漂移中心移至第304~305位。负漂移重力中心广泛分布于第61~234位，多数位于RB亚结构域，且与月份呈负相关（相关系数?0.43）。值得注意的是，4—6月期间负漂移在相对保守的VE亚结构域出现独特模式，提示保守区域的进化动态可能受季节性调节机制影响。
**5. 讨论与结论总结**
讨论部分指出，基于半协方差方法计算正负漂移重力中心，可识别与待预防病毒相似且已有获批疫苗的多种流感病毒，将这些现有疫苗组合使用，可能弥补新疫苗临床试验相对病毒流行的持续滞后。该方法避免了手工系统发育推断的低效和深度突变扫描对大量数据及计算资源的依赖。研究局限性包括：样本量有限（23条序列）、仅涵盖中国大陆区域、仅考虑一维序列信息而未纳入三维结构特征、未对神经氨酸酶（NA）共进化进行分析、缺乏体外功能实验验证。未来研究方向将扩大地理范围和病毒亚型，纳入NA共进化分析，结合HA和NA三维结构特征，并采用人工智能方法优化建模。结论部分总结：HA5蛋白序列的平均值和重力中心位置可预测与免疫逃逸相关的突变特征；负重力中心跨越第61~234位氨基酸，该区域频繁突变导致抗原不稳定或抗原漂移，使病毒逃避已有群体免疫；正重力中心位于第303~305位及其邻近区域，进化中高度保守，是疫苗和抗体开发的理想靶点。半协方差系数分析不仅可用于预测A型流感病毒HA5蛋白的免疫逃逸潜力，还可为相关病毒及毒株的疫苗和抗体开发提供依据。
