该研究由monica l. fernández-quintero、clarissa a. seidler、patrick k. quoika和klaus r. liedl等人完成，主要作者来自奥地利因斯布鲁克大学普通、无机与理论化学系及分子生物科学中心（CMBI）。论文发表于Frontiers in Molecular Biosciences，于2021年4月20日在线发表。

本研究属于结构生物学与计算生物物理学领域，聚焦鲨鱼免疫球蛋白新抗原受体（IgNAR）的可变域（vNAR）在亲和力成熟过程中的构象动力学变化及其抗原结合机制。鲨鱼作为现存最古老的有颌类脊椎动物，其适应性免疫系统具有独特的重链-only抗体形式，即IgNAR。IgNAR由两条重链通过二硫键形成同源二聚体，其可变域vNAR可独立折叠并识别抗原。与常规哺乳动物IgG抗体相比，vNAR仅含两个互补决定区（CDR）——CDR1和CDR3，其中CDR3特别长且结构复杂，可深入抗原的隐蔽表位。这些特性使vNAR在稳定性、溶解性以及穿越血脑屏障等方面表现出优异的生物技术应用潜力。然而，关于亲和力成熟如何影响vNAR的构象柔性、其与抗原的结合机制属于构象选择还是诱导契合，以及结合过程中的动力学与热力学驱动力，尚缺乏系统的原子级描述。因此，本研究旨在利用增强采样分子动力学模拟结合马尔可夫状态模型，对原始germline型与体细胞成熟型pbla8 vNAR在有无抗原状态下的构象系综进行系统表征，并阐明其抗原识别与结合机制。

本研究的工作流程可细分为四个主要步骤。第一步是体系准备与常规分子动力学模拟。研究以四种已解析的晶体结构为起始点，分别是naive型vNAR与抗原溶菌酶复合物（PDB编号2i26）、naive型vNAR游离态（2i27）、成熟型vNAR与溶菌酶复合物（2i25）以及成熟型vNAR游离态（2i24）。成熟型pbla8 vNAR较其祖先序列共携带13个点突变，其中4个位于CDR1、2个位于HV2环、1个位于HV4环、1个位于CDR3环。所有体系在水盒子中溶剂化，采用AMBER ff14SB力场，并经过多步平衡。第二步是对CDR1与CDR3环的增强采样。由于CDR环的构象转变通常发生在微秒至秒时间尺度，远超常规分子动力学模拟的覆盖范围，作者采用well-tempered metadynamics方法，以CDR1和CDR3环的ψ二面角为集合变量，对四种体系各进行了1微秒的模拟。随后对所得轨迹进行层次聚类，以1.5 Å RMSD为截断得到大量簇代表构象。第三步是基于MSM的构象系综分析。从聚类所得的代表构象出发，每个构象进行100 ns常规分子动力学模拟，累积模拟时间从1.0 μs到27.8 μs不等。利用时间滞后独立成分分析（TICA）降维，并结合马尔可夫状态模型（MSM）重构自由能面与动力学网络。第四步是针对抗体-抗原结合路径的增强采样。作者以vNAR与溶菌酶质心距离为集合变量，对naive型与成熟型复合物各进行3条独立的非well-tempered metadynamics模拟，每条10 ns。随后同样进行聚类和100 ns短轨迹回填，并以天然接触的逆距离为输入特征构建TICA和MSM，以解析结合/解离过程中的亚稳态与自由能垒。此外，还通过cpptraj计算了结合界面的静电与范德华相互作用能量，以定量比较不同状态下的驱动力构成。

研究获得了以下几方面的重要结果。第一是关于亲和力成熟导致的刚性化。无论在有无抗原条件下，成熟型vNAR的CDR1/ CDR3构象簇数量均显著低于naive型，表明构象空间大幅收窄。自由能面显示naive型呈现多个宽浅极小值，而成熟型则呈现单一深窄极小值，且该极小值对应于结合活性构象。在分子机制上，成熟型中由CDR1的Arg28与CDR3的Asp93之间形成的盐桥及氢键网络在95%以上的模拟帧中存在，而在naive型中对应的Asn28与Asp93仅形成2%的氢键，从而解释了成熟型CDR3环稳定性提高的原因。第二是关于结合机制模式的界定。作者发现，在没有抗原存在的情况下，naive型与成熟型vNAR的结合活性构象均已预先存在于溶液中，其中成熟型结合活性构象占据主导地位，而naive型则以低概率存在。这一观察支持构象选择机制作为vNAR抗原识别的主导范式。但由于成熟型vNAR在结合前后结构重排极小，其抗原识别过程可被视为“锁钥结合”；而naive型则表现出更明显的构象选择特征。第三是关于结合路径的动力学与热力学描绘。MSM在naive型与成熟型的结合路径上均识别出三个亚稳态：完全未结合态、encounter complex（遭遇复合物）和天然复合态。在naive型中，遭遇复合物为最高占据态，其结合界面静电相互作用显著强于天然复合态；而在成熟型中，天然复合态占据主导，种群从naive型的38%上升到成熟型的75%。结合能垒在成熟型中大幅增加，表明解离速率下降。第四是关于结合界面的关键残基作用。成熟型vNAR中Arg61、Asp51以及Asp101被鉴定为分子识别的关键残基。其中Asp51在成熟前为Lys51，成熟后通过离子键与溶菌酶Arg21（R136）及Lys96（K211）相互作用，并形成与Tyr20（Y135）的氢键。静电相互作用尤其在遭遇复合物形成阶段起主导作用，负责预定向与远程识别，而范德华相互作用则在遭遇复合物向天然复合态过渡时发挥更显著的作用。

本研究的结论是：亲和力成熟通过13个点突变的协同效应，显著降低了vNAR CDR环的构象多样性，并诱导结合活性构象的种群转移，使天然复合态成为最优势态。这一刚性化过程伴随着分子内稳定网络的强化与结合界面接触数量的增加。在结合机制上，vNAR对溶菌酶的识别遵循两步骤机制：首先由长程静电作用驱动的遭遇复合物形成，随后通过侧链重排和去溶剂化过渡到天然复合物。naive型抗体的结合路径偏向遭遇复合物状态，而成熟型则显著偏向天然复合态，表明亲和力成熟实质上是结合自由能景观的再塑造过程。该研究的科学价值在于系统地阐明了鲨鱼vNAR亲和力成熟的原子级机制，并定量描述了抗体-抗原结合的中间态与能垒，这对于理解适应性免疫系统的分子进化具有重要意义。在应用层面，该工作为vNAR作为新型治疗性支架的理性设计提供了结构动力学基础，尤其有助于通过工程化手段稳定vNAR结合位点、优化其亲和力与特异性，并可能推动其在穿越血脑屏障药物递送系统中的应用。本研究的亮点包括：其一，在同一序列框架内直接比较naive型与成熟型vNAR的构象系综，明确揭示了亲和力成熟的柔性-刚性转变；其二，结合metadynamics增强采样与MSM分析，在微秒至毫秒时间尺度上重建了结合路径的完整自由能景观；其三，鉴定了亲和力成熟过程中关键残基突变（如Lys51Asp）在分子识别中的具体作用，为后续定向突变提供了可操作靶点。此外，该研究还提出了一个结合机制的示意模型，直观概括了naive型与成熟型抗体在遭遇复合物与天然复合态之间的种群分布差异，对后续蛋白质-蛋白质对接算法改进和抗体工程研究具有参考价值。