SARS-CoV-2刺突蛋白Spike RBD N481K突变的功能与结构特征分析

《Archives of Virology》：Functional and structural characterization of the SARS-CoV-2 spike N481K mutation

　　

背景：选择性压力驱动SARS-CoV-2进化，刺突蛋白受体结合域（Receptor-Binding Domain, RBD）的氨基酸取代影响与人血管紧张素转换酶2（Angiotensin-Converting Enzyme 2, ACE2）的结合及抗体逃逸。N481K RBD突变2020年早期在卡塔尔出现，后在全球再现并于2024–2025年在流行毒株中接近固定，但其功能意义此前未被探究。方法：研究人员分析了超1700万条SARS-CoV-2基因组（2020–2025）中N481K的时空分布；采用建模、分子对接及分子动力学（Molecular Dynamics, MD）模拟评估结构效应；利用水泡性口炎病毒（Vesicular Stomatitis Virus, VSV）假病毒中和试验检测感染康复者血清（n＝48）及接种者血清——mRNA疫苗（n＝16）、国药（Sinopharm）灭活疫苗（n＝45）、阿斯利康（AstraZeneca）腺病毒载体疫苗（n＝28）的中和能力。结果：N481K在全球基因组中检出率为4.14％（703,777/17,016,111），2024年起在主流Omicron谱系中接近固定（＞90％），2024年初71.5％升至2025年第36周94.4％。结构分析显示野生型（Wild-Type, WT）RBD与N481K RBD均方根偏差（Root Mean Square Deviation, RMSD）为0.706 ?，N481K与ACE2对接亲和性更高（－313 vs. －308 kcal/mol）；MD确认N481K–RBD–ACE2复合物稳定性增加，结合自由能更低（N481K：－77.07±0.90 vs. WT：－43.09±0.98 kcal/mol）。在WT中和可检样品（n＝89）中，平均中和率由WT的69.76％降至N481K的42.47％（p＜0.0001），阿斯利康、国药及感染组降幅最大，mRNA组保留最高中和率［WT 77.5％ vs. N481K 56.0％］。结论：N481K增强ACE2结合同时部分逃逸抗体中和，提示可提高病毒适应度与持续传播能力。

SARS-CoV-2刺突蛋白Spike RBD N481K突变的功能与结构特征研究解读——《Archives of Virology》论文浓缩总结

研究背景与立项依据

严重急性呼吸综合征冠状病毒2型（SARS-CoV-2）持续通过刺突（Spike, S）蛋白氨基酸突变进化，其中受体结合域（Receptor-Binding Domain, RBD；残基306–534）的取代直接调控与人血管紧张素转换酶2（Angiotensin-Converting Enzyme 2, ACE2）的相互作用及中和抗体识别。深突变扫描表明多数RBD逃逸突变需在抗体逃逸与受体亲和力间权衡，而少数组合可同时恢复或增强亲和力。经典RBD突变（如K417N/T、E484K、L452R、T478K、N501Y）已被广泛研究，但位于受体结合基序（Receptor-Binding Motif, RBM）内的N481K（Asn₄₈₁→Lys，引入正电荷可能影响RBM–ACE2界面静电环境及抗体可及性）虽自2020年起零星检出、2023年后随Omicron后代谱系全球复现并接近固定，其结构与功能后果长期未获充分研究。为此，研究人员整合基因组监测、计算生物学模拟及假病毒中和实验，系统阐明N481K突变的流行特征、结构影响及免疫逃逸表型。

主要关键技术方法概述

研究人员利用CoV-Spectrum平台分析2020–2025年全球＞1700万条SARS-CoV-2基因组中S:N481K的时空与谱系分布；通过Chimera软件模建N481K RBD三维结构并与野生型（Wild-Type, WT）比对均方根偏差（Root Mean Square Deviation, RMSD），使用HDOCK服务器进行RBD–ACE2分子对接并由PDBsum分析氢键/盐桥/非键接触；采用AMBER20（FF19SB力场）运行200 ns全原子分子动力学（Molecular Dynamics, MD）模拟，经CPPTRAJ计算RMSD、回转半径（Radius of Gyration, Rg）、残基均方涨落（Root Mean Square Fluctuation, RMSF）及MM/PBSA（Molecular Mechanics/Poisson–Boltzmann Surface Area）结合自由能；收集感染康复者（n＝48）及三种疫苗免疫者［mRNA疫苗（Pfizer BNT162b2＋Moderna mRNA-1273，n＝16）、国药（Sinopharm BBIBP-CorV，n＝45）、阿斯利康（AstraZeneca ChAdOx1 nCoV-19，n＝28）］血清，通过位点定向诱变法构建N481K Spike表达质粒，包装VSV（Vesicular Stomatitis Virus）ΔG-荧光素酶假病毒，开展配对血清假病毒中和试验并以配对统计检验比较WT与N481K中和差异。

研究结果

Temporal and geographic circulation of S: N481K (2020–2025)

在17016111条序列中检出N481K阳性703777条（4.14％）。2020年初美国首检，卡塔尔2020年4月出现并在B.1谱系中局部流行（至2021年1月＞50％），全球其余地区此前罕见（＜1％）；2023年第38周起占比突破1％并骤升，2023年末＞60％，2024年初71.5％升至2025年第36周94.4％（第37周略降至～87.5％）。早期检出集中于中东（卡塔尔、巴勒斯坦），2023年后扩散至美、欧、亚。2025年1–9月该突变出现于93.23％测序样本，在JN.1（11.72％）、KP.3.1.1（8.34％）中富集，XEC、LP.8.1.1、NB.1.8.1等谱系携带率最高。卡塔尔2020年4月至2021年1月曾流行后消失，2023年末再现已于2024年起100％检出。

Molecular docking and structural modeling of the S: N481K mutation

N481K RBD相对WT RBD结构叠合RMSD为0.706 ?，提示局部构象微调可能影响ACE2互作。分子对接显示WT RBD–ACE2得分－308 kcal/mol（12条氢键、1个盐桥Lys₄₁₇–Glu₃₀、170个非键接触）；N481K RBD–ACE2得分－313 kcal/mol（13条氢键、同盐桥、155个非键接触），提示N481K提高RBD–ACE2结合亲和性。

Dynamic stability and compactness analysis of wild-type and mutant complexes

200 ns MD显示WT复合物约5 ns平衡后波动较大（平均RMSD≈4 ?），N481K复合物2 ns即平衡且全程收敛更稳（平均RMSD≈3 ?）；回转半径Rg（WT≈31.2 ?，N481K≈31.0 ?）表明N481K复合物略紧凑。较低RMSD与Rg支持N481K–RBD–ACE2复合物动态稳定性增强。

Residual fluctuation and hydrogen bonding analysis of wild-type and mutant complexes

平均氢键数相近（～380条），N481K复合物整体残基均方涨落（Root Mean Square Fluctuation, RMSF）低于WT，提示RBM区柔性降低、界面更刚硬，可能不利于抗体结合。

Binding free energy calculation

MM/PBSA计算总结合自由能：WT RBD–ACE2为－43.09±0.98 kcal/mol，N481K RBD–ACE2为－77.07±0.90 kcal/mol；N481K的范德华（－91.17±0.78 vs. －76.85±0.67 kcal/mol）与静电（－804.61±4.64 vs. －577.02±5.29 kcal/mol）贡献均更强，证实突变显著降低结合自由能、增强RBD–ACE2亲和力。

Functional assessment of the neutralization activity against the N481K mutant PV

全体137份血清WT假病毒平均中和率48.5％，N481K假病毒降至37.9％（p＝0.0002）。在WT中和阳性（抑制≥30％，n＝89）亚组中，WT平均中和69.76％（95％ CI 65.86–73.67），N481K平均42.47％（95％ CI 38.58–46.37）（p＜0.0001）。分组：感染康复组WT 67.1％→N481K 41.2％；mRNA组WT 77.5％（95％ CI 69.7–85.3）→N481K 56.0％（95％ CI 50.6–61.5），中和保留最好（≥50％ cutoff下62.5％仍达标）；国药组WT 59.5％→N481K 30.7％；阿斯利康组WT 73.5％→N481K 39.9％。N481K致所有组别中和显著下降，非mRNA疫苗诱导血清降幅更大（阿斯利康平均降48.0％，感染康复平均降40.7％，国药平均降43.8％，mRNA平均降28.7％）。

讨论与结论

研究人员指出N481K早期局部流行后随Omicron谱系全球再扩张并近固定，暗示强正向选择带来的适应度优势。结构—动力学证据共同表明N481K引起RBD轻微构象偏移、RBM柔性降低、RBD–ACE2复合物更稳定且结合自由能更低，即增强ACE2亲和力；假病毒中和试验证实该突变使各免疫背景血清中和效价显著下降，尤以灭活及腺病毒载体疫苗和天然感染血清明显，mRNA疫苗诱导广谱中和抗体部分抵抗此逃逸。N481K兼具增强受体结合与部分抗体逃逸双重表型，符合SARS-CoV-2适应性进化特征，可能与其它高亲和力RBD突变呈协同而非单纯补偿作用。

结论翻译：综上所述，本研究表明Spike RBD N481K突变与预测的Spike–ACE2结合亲和力升高及抗体介导中和能力下降相关，提示其对SARS-CoV-2病毒适应度与持续存在具潜在协同增益作用。流行病学、结构及功能数据支持N481K对SARS-CoV-2持续进化的适应性贡献。受限于计算机模拟及中和样本量局限，后续需结合高分辨率结构分析与RBD内外共突变模型，进一步明确N481K与其它RBD及非RBD突变组合的叠加影响。