SARS-CoV-2刺突蛋白N端结构域（NTD）抗原重塑：序列变异与中和抗体超位点结构重组的关联研究

《Frontiers in Cellular and Infection Microbiology》：Antigenic remodeling of the SARS-CoV-2 Spike N-terminal domain links sequence variability to structural reorganization of a neutralizing antibody supersite

　　

摘要 引言： 严重急性呼吸综合征冠状病毒2型（SARS-CoV-2）刺突蛋白（Spike）的N端结构域（N-terminal domain，NTD）含有一个中和抗体超位点（neutralizing antibody supersite），且是抗原性变化的频繁发生位点。然而，疫苗接种前后免疫景观中的序列变异性如何重塑抗体识别，目前仍未得到充分表征。 方法： 研究人员将临床基因组监测与结构建模和分子动力学（molecular dynamics，MD）相结合，以识别具有抗原相关性的NTD变化，并探究其对中和抗体界面的结构影响。研究人员对两个反映不同免疫背景时期的SARS-CoV-2谱系进行了比较：以B.1.1.33为主的疫苗接种前时期和以Omicron亚谱系BQ.1.1为主的疫苗接种后时期。一个本地队列（n=126）支持描述性表征。一个整合了本地和GISAID基因组并统一了临床元数据的扩展队列（n=534；401例疫苗接种前，133例疫苗接种后）被用于根据突变与临床严重程度的关联来优先筛选候选位点。 结果： 在调整模型中，年龄在两个时期均与严重程度独立相关，而性别则无此关联。在疫苗接种前队列中，NTD中的Spike S突变与更严重的疾病相关，但未观察到按临床结局的系统发育分离。该关联属于探索性发现，被用于优先筛选一个位于抗原区域的NTD变异进行结构分析。Q23K在圣埃斯皮里图州（Espírito Santo）呈现区域性富集（57.4%），而在巴西其他州（2.7%）和国际基因组（1.4%）中频率较低。在疫苗接种后队列中，一个重叠的缺失事件（ORF9b Δ27-29/N Δ31-33）与非重症疾病相关，但鉴于元数据异质性、队列结构及该位点的性质，应谨慎解读。然而，在该时期，临床结局与测序来源存在部分混杂。一项仅限于GISAID样本的按来源分层的敏感性分析显示该关联仍然存在，但无法排除来源或流程依赖的插入缺失（indel）判定伪影，因此该发现被报告为探索性结果。尽管Q23K本身未持续存在，但相邻的NTD变化（T19I/A27S/Δ24-26）在疫苗接种后变得频繁出现，表明同一抗原邻域发生了重塑，而非单一突变的持续存在。NTD与抗体1-87复合物的结构建模和分子动力学模拟支持所建模的抗体-抗原界面发生重组，表现为超位点接触的重新分布，而非结构相容性的丧失。 讨论： 这些发现将不同免疫背景下的NTD变异性与所建模的抗体-抗原界面的结构重组联系起来，并支持将基因组监测与结构分析相结合，以优先筛选具有抗原相关性的变异进行功能验证。

论文解读：SARS-CoV-2刺突蛋白N端结构域抗原重塑的临床基因组学与结构生物学整合研究

一、研究背景与科学问题

严重急性呼吸综合征冠状病毒2型（SARS-CoV-2）的持续演化对疫情预测和控制构成重大挑战。病毒刺突蛋白（Spike）是病毒入侵宿主细胞的主要决定因素，同时也是中和抗体的主要靶标。尽管受体结合结构域（receptor-binding domain，RBD）受到广泛关注，但N端结构域（N-terminal domain，NTD）同样是重要的抗原区域。NTD含有一个被强效中和抗体靶向的超位点（supersite），并且是反复出现的缺失和替换突变的热点区域，这些突变常与抗体逃逸相关。然而，在疫苗接种前后的不同免疫背景下，NTD序列变异如何重塑抗体识别机制，此前尚未得到充分阐明。

巴西自2021年1月开始大规模疫苗接种，为研究不同免疫选择压力下病毒突变谱的变化提供了独特背景。本研究旨在通过整合临床流行病学数据、基因组监测数据和结构生物学分析，系统探究疫苗接种前后两个不同时期NTD变异与中和抗体识别之间的关系。

二、研究设计与主要方法

研究人员选取巴西圣埃斯皮里图州（Espírito Santo）的两个时间窗口进行分析：疫苗接种前时期（2020年7月9日至2021年3月29日），以B.1.1.33谱系流行为主；疫苗接种后时期（2022年11月14日至2022年12月28日），以Omicron亚谱系BQ.1.1流行为主。研究构建了两个分析数据集：本地队列（n=126）用于描述性临床基因组表征，扩展队列（n=534）整合本地基因组和GISAID公共基因组，用于严重程度相关突变信号的推断性分析。研究采用的主要技术方法包括：全基因组测序（Oxford Nanopore平台）、最大似然系统发育分析（IQ-TREE）、AlphaFold3结构预测、GROMACS分子动力学模拟（CHARMM36力场，每个体系3次独立200纳秒重复）、PLIP界面相互作用预测、PRODIGY结合自由能估算以及FreeSASA埋藏表面积计算。

三、主要研究结果

3.1 本地队列的临床和流行病学特征

本地队列中，疫苗接种前时期包含53个基因组（11例轻症，42例重症），疫苗接种后时期仅包含73例轻症。疫苗接种前重症患者年龄较大，以呼吸系统表现为主，呼吸困难（dyspnea）和去饱和（desaturation）为主要症状，中位住院时间16.5天，病死率54.8%。疫苗接种后队列以女性为主（72.6%），症状以非特异性表现为主。两个时期的主要流行谱系分别为B.1.1.33（60.4%）和BQ.1.1（75.3%）。

3.2 基因组分析与突变谱

扩展队列分析显示，疫苗接种前B.1.1.33组（n=401）存在593个独特的氨基酸变化，呈长尾分布；疫苗接种后BQ.1.1组（n=133）存在388个独特变化，突变在Spike中更为集中。在调整年龄和性别后，年龄在两个时期均与严重程度独立相关。疫苗接种前时期，Q23K突变在重症中富集（11.0% vs 2.0%），调整后优势比（odds ratio，OR）为6.94（95% CI 1.83-26.28，p=0.004），经Benjamini-Hochberg校正后仍显著（q=0.004）。系统发育分析显示重症和轻症Q23K阳性基因组在多个分支上交错分布，未形成结局特异性进化枝。地理分布分析表明Q23K在圣埃斯皮里图州显著富集（57.4%），远高于巴西其他州（2.7%）和国际序列（1.4%）。疫苗接种后时期，ORF9b:Δ27-29/N:Δ31-33重叠缺失事件在轻症中更常见（92.6% vs 49.2%），调整后OR为0.038（95% CI 0.011-0.13），但该发现受测序来源混杂影响，被报告为探索性结果。

3.3 结构建模与分子动力学模拟

研究人员以抗体1-87作为结构探针，对野生型（WT）、Q23K和T19I/A27S/Δ24-26（LPP）三种NTD模型进行了分子动力学模拟。结果显示，与WT相比，突变体复合物表现出更紧凑的构象：LPP的均方根偏差（root-mean-square deviation，RMSD）最低（0.40±0.07 nm），WT最高（0.64±0.11 nm）；突变体的溶剂可及表面积（solvent-accessible surface area，SASA）和回旋半径（radius of gyration，Rg）均低于WT。埋藏表面积（buried surface area，BSA）从WT到Q23K和LPP依次降低（30.5±6.42、26.0±3.3、21.7±2.0 nm2），效应量最大且最可重复（?2=0.563）。接触持续性分析表明，WT界面较为弥散，仅有一个超稳定接触；Q23K增加了接触数量和持续性；LPP虽总接触数较少，但具有最高的中位占有率（8.7%）和最多的超稳定（12个）及高稳定（5个）接触。关键锚定残基包括WT中的ASP253-GLY:B:411、Q23K中的TRP152-PRO:B:414，以及LPP中的LYS147-GLU:B:379和GLY252-ASP:B:423。

四、讨论与结论

本研究通过整合临床基因组监测和结构生物学分析，揭示了SARS-CoV-2 NTD在不同免疫背景下的抗原重塑机制。Q23K作为疫苗接种前时期的区域性富集突变，虽与严重程度存在探索性关联，但其生物学意义主要在于抗原性而非受体结合。疫苗接种后时期，虽然Q23K本身未持续存在，但同一抗原邻域的LPP基序（T19I/A27S/Δ24-26）成为高频变异，表明NTD超位点经历了持续的接触网络重组。结构模拟表明，突变体并未丧失与抗体1-87的结构相容性，而是通过接触重新分布和巩固形成了更紧凑的相互作用模式。这种重组可能通过改变表位几何结构、溶剂可及性或糖基化遮蔽效应影响抗体的中和效力，但这些机制超出了当前单体、单抗体、预形成复合物模型的解释范围。研究结论支持将基因组监测与结构分析相结合，以优先筛选具有抗原相关性的变异进行功能验证，为理解SARS-CoV-2持续演化背景下抗体识别的动态变化提供了重要框架。