《International Journal of Molecular Sciences》：Structural and Biophysical Profiling of HLA Class II Mismatches Using AlphaFold 3: A Computational Study
编辑推荐：
传统线性表位匹配算法（如HLAMatchmaker，HLAMM）在评估人类白细胞抗原（human leukocyte antigen，HLA）位点的氨基酸序列多态性时，未整合三维（three-dimensional，3D）空间可及性和局部物理化学约束。本研究开
传统线性表位匹配算法（如HLAMatchmaker，HLAMM）在评估人类白细胞抗原（human leukocyte antigen，HLA）位点的氨基酸序列多态性时，未整合三维（three-dimensional，3D）空间可及性和局部物理化学约束。本研究开发了一套探索性计算抗原性分析流程，利用AlphaFold 3评估10,000对模拟欧洲供者-受者配对，覆盖预先计算的、全长HLA I类（HLA-A、-B、-C）和异二聚体II类（HLA-DR、-DQ、-DP）复合物，这些复合物以标准化的无肽构象进行建模。候选表位依次通过空间溶剂可及表面积阈值（SASA ≥ 10.0 ?2）和经验性生物物理排除门控（绝对净电荷|Qnet|的第5百分位数 ≥ 0.3326 e且总亲水性Htotal ≥ ?7.265，基于N = 170个抗体验证表位校准）进行筛选。在模拟队列（n = 10,000）中，过滤将中位数[Q1；Q3]表位计数从19 [12；26]降至17 [10；23]（HLA-A），22 [15；29]降至18 [12；24]（HLA-B），12 [7；17]降至10 [7；15]（HLA-C），15 [10；19]降至11 [7；14]（HLA-DR），12 [6；17]降至9 [4；14]（HLA-DQ），以及7 [3；12]降至6 [2；9]（HLA-DP）。组分消融分析表明，3D空间可及性（SASA ≥ 10.0 ?2）贡献了绝大多数表位排除——在HLA-DR和HLA-DP中占87.8%，在HLA-DQ和I类位点中占58.2%至69.1%——而注册表校准的生物物理门控作为次级精筛发挥作用。在最初超过任意≥10表位线性阈值的模拟配对子集中，过滤后降至10个表位以下的比例在I类中较低（HLA-A：4.4%，HLA-B：8.0%，HLA-C：13.2%），但在II类高负荷子集中显著更高：HLA-DQ（21.8%）、HLA-DR（26.0%）和HLA-DP（3693对基线高负荷配对中的33.7%）。此外，II类位点，特别是HLA-DP，在1:1计数匹配子采样（ρmatched = ?0.268）和偏相关调整（ρpartial = ?0.385）后，启发式电荷与亲水性指数之间的协方差仍然减弱。该双过滤流程定义了一个受限的数学子集，包含空间可及且生物物理超阈值的候选表位。由于这是一个假设生成性的计算机模拟框架，在临床转化之前，严格需要多中心回顾性注册研究以及评估新生供者特异性抗体形成、活检证实的排斥反应和移植物存活的前瞻性临床试验。
论文解读
一、研究背景与意义
在肾移植同种免疫风险分层领域，传统的血清学抗原分型已逐步演进为基于表位（eplet）水平的分子相容性精细评估。表位是指HLA分子表面能够引发抗体识别的离散多态性斑块。以HLAMatchmaker（HLAMM）为代表的线性序列比对算法已成为评估供者特异性抗体（donor-specific antibody，DSA）形成风险和预测移植物结局的标准工具。尽管PIRCHE-II等补充框架增强了对间接T细胞同种反应性的预测，但准确定量B细胞相关的表面可及性仍是持续性挑战。
然而，这些广泛采用的方法从根本上依赖于HLA一级氨基酸序列的一维（one-dimensional，1D）分析。这一简化忽略了HLA分子的固有三维（three-dimensional，3D）结构，可能将埋藏、结构受限或空间不可及的残基错误分类为具有免疫原性，导致风险评估膨胀，并可能排除潜在相容的供者。近年来，Snowflake算法、HLA3D工具包和HLA-EpiCheck等计算机模拟方法尝试整合3D结构数据，但这些工具往往依赖异质的现有X射线晶体数据集、主要作为可视化数据库运作而不具备定量临床风险重分类能力，或需要当前临床交叉配型实践中不切实际的大量计算资源。
在此背景下，研究人员严格区分了抗原性（antigenicity）与免疫原性（immunogenicity）两个操作概念。抗原性被定义为多态性表面斑块在理论上结合抗体的结构可及性和物理化学能力；免疫原性则涵盖导致可测量临床结局（如DSA形成或移植物排斥）的实际体内生物学级联反应。由于静态计算机模拟模型无法模拟完整的生物学反应，该计算框架严格设计用于评估理论抗原性。AlphaFold 3的引入标志着结构免疫生物学的重要进展，其利用生成式扩散架构实现近原子精度的生物分子复合物预测，在重建II类异二聚体（HLA-DR、-DQ、-DP）关键的蛋白-蛋白界面方面表现出显著改善，从而实现了从理论多态性替换枚举向分子表面严格物理化学表征的转变。
二、研究设计与主要技术方法
研究人员开发了一套自动化探索性组织相容性评估流程，核心基于AlphaFold 3生成的全长结构模型。研究使用官方AlphaFold 3服务器构建了包含80个模型的HLA结构库（15个HLA-A、15个HLA-B、14个HLA-C、11个HLA-DR、10个HLA-DQ和15个HLA-DP复合物），其中76个高频欧洲等位基因和α/β单倍型（群体频率≥1.0%）用于10,000对供者-受者配对的蒙特卡洛模拟。所有I类重链-β2-微球蛋白复合物和II类α/β异二聚体均以标准化无肽构象共同建模，以避免序列特异性肽段偏差并提供溶剂可及表面积（solvent-accessible surface area，SASA）的保守上界估计。结构验证通过对13个模型与Protein Data Bank（PDB）高分辨率晶体结构的Cα原子均方根偏差（root mean square deviation，RMSD）叠加比对完成。表位负荷计算采用三种方法：标准HLAMM方法、仅验证表位方法，以及研究开发的双过滤结构-生物物理方法。后者依次应用空间可及性约束（SASA ≥ 10.0 ?2，采用Shrake-Rupley算法）和经验性生物物理排除门控（基于170个抗体验证表位校准的第5百分位数绝对净电荷和总亲水性阈值）。物理化学特征分析采用Henderson-Hasselbalch方程计算pH 7.4下的净电荷，以及Kyte-Doolittle亲水性标度计算总疏水性潜力。
三、研究结果
3.1 结构验证与模型质量
内部质量评估显示，整个结构库具有一致的高预测置信度，总体平均pLDDT（predicted Local Distance Difference Test）评分为86.50 ± 0.82，折叠的细胞外抗原呈递结构域的中位数残基水平pLDDT为95.1-98.6。总体中位数pTM（predicted Template Modeling）评分为0.81，总体中位数ipTM（interface pTM）评分为0.94（I类）和0.78-0.80（II类）。与实验结构比对显示，13个测试等位基因的平均Cα RMSD为0.56 ± 0.20 ?，所有参考结构的Cα偏差均低于0.90 ?，确认了适用于溶剂可及性分析的近天然三级几何结构。
3.2 表位错配负荷的比较分析
双过滤方法在所有六个位点均系统性地降低了计算错配计数。I类中，HLA-B位点的降幅最显著（中位数从22降至18）；II类中，HLA-DR从15降至11，HLA-DQ从12降至9，HLA-DP从7降至6。双过滤模型始终占据标准HLAMM库与仅验证子集之间的中间数学位置。五臂组分消融研究证实，仅SASA结构过滤在HLA-DR和HLA-DP中贡献了总表位排除的87.8%，在HLA-DQ和I类位点中贡献58.2%至69.1%，而单独应用生物物理过滤仅产生最小减少，确认II类错配负荷下降主要由3D空间可及性驱动。
3.3 代表性II类错配的结构解卷积
以供者DRB1 * 04:01/受者DRB1 * 07:01为代表配对为例，标准算法识别出19个氨基酸错配，均匀分布于整个分子，包括肽结合槽底部。应用SASA > 10 ?2过滤后，位于深槽袋内的残基（特别是Tyr37）因被α螺旋壁空间屏蔽而低于阈值，被重新分类为结构惰性，仅保留表面暴露的α螺旋区域作为有效靶标。
3.4 生物物理与静电分析
标准HLAMM算法在所有位点均产生比双过滤模型更高的累积绝对电荷分数。被排除的可电离残基主要位于β折叠底部和肽锚定袋（P1、P4、P6、P7和P9）深处或链间界面，这一拓扑定位直接由基准结构的X射线晶体学研究确认。相比之下，有符号Kyte-Doolittle亲水性总和的分布变化高度依赖位点：HLA-DP位点的变化最大（秩双列相关r = 0.689），标准算法产生正的中位数亲水性总和，而双过滤将其移至近中性水平。
3.5 固定截断分类与物理化学协方差
在初始超过≥10表位截断的配对子集中，I类位点显示高一致性（仅4.42%的HLA-A、7.95%的HLA-B和13.24%的HLA-C高负荷配对过滤后降至截断以下），而II类高负荷子集的降幅显著更高：HLA-DQ为21.84%，HLA-DR为25.99%，HLA-DP为33.74%。Cohen加权kappa系数确认了位点依赖性二分法：HLA-A、HLA-B、HLA-C和HLA-DQ显示“几乎完美”一致性（κw > 0.80），而HLA-DR和HLA-DP仅显示“实质性”一致性（κw < 0.80）。电荷与亲水性协方差分析显示，大多数位点存在强负相关（HLA-A ρ = ?0.842），但HLA-DP表现出最弱相关性（ρ = ?0.400）。1:1计数匹配子采样和偏相关分析确认，即使在相同计数分布下，HLA-DP仍保留最弱的电荷-亲水性关联（ρmatched = ?0.268），表明该位点的静电势和疏水性作为更独立的理论抗原性驱动因素。
四、讨论与结论
讨论部分指出，线性算法假设所有多态性残基具有同等的免疫可及性，而研究结果证明这一简化引入了计算噪声。通过整合AlphaFold 3的高保真结构能力，研究实现了显著的“结构去噪”效应，将空间可及、生物物理超阈值的候选表位与隐蔽、空间屏蔽残基分离。这一发现与Tambur等人强调的差异免疫原性模型一致，为表位免疫原性的经验性临床层级提供了机制性解释。研究还强调了I类与II类之间的显著方法学差异：I类分子的构象刚性导致线性算法与结构建模高度一致，而II类异二聚体约三分之一的“高风险”配对被结构分析重新分类为较低风险。研究严格区分了供者选择与移植后免疫抑制管理两种临床应用场景，强调后者需要更高水平的临床证据门槛。局限性包括：结构库局限于欧洲高频等位基因、静态无肽模型无法捕捉构象可塑性、未建模糖基化等生理变量、未整合HLA-DP表达水平相关遗传变量（如rs9277534）、生物物理阈值校准的部分循环性，以及序列衍生点电荷和亲水性总和作为热力学结合强度替代指标的固有局限性。
结论部分指出，将AlphaFold 3纳入组织相容性评估实现了从线性序列匹配向高分辨率结构分析的转变。通过整合生物物理约束——溶剂可及性、静电势和疏水性——该流程算法性地过滤了理论识别的表位，特别是针对复杂的HLA II类位点。预计算结构库的生成和存储允许在护理点无需高性能计算即可快速评估供者-受者配对。研究人员假设该计算框架提供了HLA错配更全面的生物物理特征描述，但严格需要跨纵向临床移植队列的回顾性验证和后续前瞻性临床试验，以确定这些理论上精炼的结构特征是否与实际同种免疫和移植物存活结局相关。
