X连锁高IgM综合征中CD40L错义变异A123E、S222F和G257R的计算分析预测其改变了CD40结合与同源三聚体稳定性

《Frontiers in Immunology》：Computational analysis of CD40L missense variants A123E, S222F, and G257R predicts altered CD40 binding and trimeric stability in X-linked hyper-IgM syndrome

　　

CD40L（CD154）是肿瘤坏死因子（TNF）超家族的同源三聚体蛋白成员，表达于活化的CD4⁺ T细胞上，其与CD40的结合对于抗体的类别转换重组（CSR）、体细胞高频突变以及B细胞的终末分化是必不可少的。基因CD40LG的功能丧失突变会导致X连锁高IgM综合征1型（HIGM1），其临床特征为反复感染、缺乏IgG、IgA和IgE，以及IgM水平接近正常或升高。研究人员采用计算策略描述了天然蛋白及三种此前在分子水平上未被表征的错义变异（A123E、S222F和G257R）所对应的CD40L-CD40复合物的结构与动态特征。每个系统在300 K和0.15 M NaCl的显式溶剂中均进行了两次独立的300 ns模拟。此外，研究人员利用连续溶剂模型评估了该复合物的能量学与静电学特征，如分子力学-泊松-玻尔兹曼表面积（MM-PBSA）和自适应泊松-玻尔兹曼求解器（APBS）。研究发现，A123E导致了中等程度但能量学上一致的结合亲和力降低，这与部分功能丧失相符。S222F导致了结合自由能的最大降低（根据MM-PBSA为最弱的结合物），同时伴随着CD40L-CD40接触面积的增加，这与几何上扩张但能量学上受破坏的界面相符，在此界面上，芳香族侧链抵消了矛盾且有利的静电重组。G257R维持了接近天然的结合自由能，但将其三聚体接触面积减少了20-22%，并增加了CD40L的回转半径，这与同源三聚体的结构松弛相符，研究人员预测这可能使结合能力与信号传导能力解偶联，对免疫突触处的CD40受体流动性和机械转导产生潜在影响，这些影响仍待实验验证。这些发现可能为HIGM1的临床异质性提供了结构基础，并展示了计算方法在解释CD40LG变异（包括在拉丁美洲患者中识别出的新变异，这些地区极少能进行实验验证）方面所能提供的贡献。

在适应性免疫中，CD40L-CD40（CD154）轴是最关键的共刺激检查点之一。CD40L属于肿瘤坏死因子（TNF）超家族，而CD40属于TNF受体超家族；它们的相互作用能够激活经典和非经典的核因子κB（NF-κB）通路，以及丝裂原活化蛋白激酶（MAPK）、磷脂酰肌醇3-激酶/蛋白激酶B（PI3K/AKT）和Janus激酶3/信号转导及转录激活因子（JAK3/STAT）信号传导。在B淋巴细胞中，这种相互作用驱动生发中心的形成、免疫球蛋白类别转换重组（CSR）、亲和力成熟以及长寿浆细胞和记忆B细胞的生成。X连锁高IgM综合征（XHIGM）由染色体Xq26.3上的CD40LG基因突变引起，其特征是T-B淋巴细胞串扰缺陷和类别转换重组失败。在拉丁美洲，由于难以获得专门诊断和功能验证检测手段，先天性免疫缺陷易被漏诊。导致HIGM1的CD40L变异其功能后果极少得到表征，因此迫切需要计算结构方法来补充传统实验手段。本研究选取了三种在文献或临床记录中报道过的错义突变：A123E、S222F和G257R。这些变异均与XHIGM表型相关，研究人员假设它们通过各自不同的机制影响结合亲和力、三聚体稳定性或信号传导能力以削弱CD40L功能。该研究论文发表在《Frontiers in Immunology》上，为理解HIGM1相关突变的分子发病机制提供了新视角，并突出了计算结构生物学在资源匮乏地区对于变异分类和分子诊断的应用价值。

研究采用PDB数据库中3QD6分辨率为3.5 ?的晶体结构作为起始结构，并利用AlphaFold2重建了缺失的短环。对天然和突变复合物系统在显式溶剂中进行了两次独立的300纳秒分子动力学（MD）模拟。主要分析技术包括：均方根位移（RMSD）、均方根涨落（RMSF）、回转半径计算、接触图谱分析、结合自由能计算（基于MM-PBSA方法）以及静电学计算（基于APBS方法）。研究中的样本队列主要来源于既往文献报道的临床病例，其中包括由DiSanto等人及Lopez-Granados等人报道的患者，以及此前在墨西哥HIGM1患者队列中发现并表现为CD40L表面表达严重降低的G257R突变病例。

结构稳定性与全局构象变化
通过RMSD分析确认所有系统在两次300 ns模拟的前50 ns内均达到平衡。RMSF分析显示，突变体在180至190号残基附近展现了局部动力学的重组。回转半径分析表明，与天然蛋白及其他突变体相比，G257R在两次独立重复中均表现出显著升高的回转半径，表明该突变引发了三聚体配体的全局结构扩张，而A123E和S222F未产生此效应，这凸显了G257R对同源三聚体组织架构的特定影响。

突变残基的局部接触环境
A123E使得与突变位点接触的残基数量增加，带负电的谷氨酸侧链周围聚集了更多残基，包括额外的碱性接触。S222F同样增加了接触残基的数量，并向疏水接触发生转变，苯丙氨酸侧链占据了更大的局部体积。G257R展现了三种变异中最大的局部接触重排，精氨酸残基占据了更大的局部体积，与天然甘氨酸相比，其伴随着接触壳的扩张、疏水接触的减少以及亲水和碱性接触的增加。

CD40L-CD40结合界面的接触图谱分析
在结合界面处，最频繁的接触涉及来自CD40L的碱性精氨酸和来自CD40的酸性谷氨酸（如复制1中的ARG-200和GLU-74）。天然系统显示出一个主要和一个次要接触对。A123E增加了精氨酸介导的次级接触频率。S222F在多个精氨酸残基间的接触频率分布更为均匀，其碱性-酸性接触的比例低于天然系统。G257R则显示出较低的次级接触频率，但同时保留了主要的ARG/GLU接触，导致高频残基对减少。

接触面积分析
在三聚体界面处，天然蛋白、A123E和S222F维持了相当的接触面积（约1947-2042 ?²）。而G257R在三聚体接触面积上表现出显著减少，与天然蛋白相比减少了20-22%，这反映了该界面的结构变异性。在CD40L-CD40分子间界面，S222F显示出一致的接触面积增加，表明苯丙氨酸在几何上扩张了界面但并未改善能量质量；G257R则表现出居中的数值，其特异性地损害了三聚体界面，而使受体结合表面积在很大程度上保持不变。

结合自由能计算
MM-PBSA计算显示，天然复合物的ΔG值在两次重复中分别为5.70和22.31 kcal/mol。A123E的ΔG值分别为18.67和28.57 kcal/mol，表现出亲和力的中等降低。S222F显示出所有系统中最高的ΔG值（分别为41.78和36.46 kcal/mol），是结合力最弱的突变体。G257R的ΔG值（18.28和16.73 kcal/mol）则接近天然值。APBS静电学计算表明，G257R具有所有系统中最高的静电亲和力，这主要归因于带正电的精氨酸残基增加了与CD40的静电互补性。将MM-PBSA和APBS结果结合起来分析，A123E表现出一致的结合亲和力降低；S222F在MM-PBSA中结合力最弱但在APBS中有所改善，表明庞大的侧链带来的范德华力和溶剂化惩罚抵消了有利的静电重组；G257R则在保持近天然结合能量的同时产生了最广泛的结构重组。

在讨论部分，研究人员总结了每种突变预测的致病机制。A123E产生了有限的扰动，其三聚体接触面积未受影响且结合亲和力仅有中等降低，这与其部分或衰减的CD40L功能丧失相符。即使CD40L-CD40相互作用仅发生部分损害，也会削弱CSR的效率，因为B细胞活化所需的CD40信号阈值对相互作用质量和持续时间高度敏感。S222F通过庞大的芳香侧链引起空间位阻，在几何上扩张了界面但在能量上降低了结合质量。其改善的静电互补性可能减慢CD40解离，从而限制了有效CD40聚集和信号传导所需的侧向受体流动性。G257R则维持了接近天然的MM-PBSA结合自由能和最有利的APBS静电学，却产生了最广泛的结构重组。其通过三种机制实现解偶联：首先，三聚体接触面积减少和回转半径增加表明同源三聚体发生结构松弛，即使单个接触能量有利，也无法驱动精确的受体聚集；其次，增强的静电互补性会减慢CD40解离，损害动态CD40聚集；最后，CD40可作为机械传感器与CD40L形成 Catch bond，G257R产生的结构扩张和三聚体紊乱会损害跨界面的力传递。针对拉丁美洲患者中发现的未表征变异G257R，本研究提供了首个机制描述，这突出了计算方法在实验验证受限环境中的独特价值。

结论部分指出，研究人员表征了三种与HIGM1相关的CD40LG错义突变。每种突变均通过不同的机制扰动CD40L-CD40复合物：A123E产生中等且能量一致的结合亲和力降低，与部分功能丧失相容；S222F通过芳香侧链的空间干扰引起结合自由能的最大降低，与实验记录的CD40结合能力丧失相符；G257R维持了近天然的结合自由能，但将三聚体接触面积减少了20-22%并扩张了回转半径，这种结构特征与同源三聚体松弛相符，预测可能使结合能力与信号传导能力解偶联。这些发现表明，HIGM1相关的CD40LG突变不必统一废止受体结合，而是可能通过破坏 productive 信号传导所需的更高阶结构和动态特性来发挥作用，这解释了HIGM1的临床异质性。同时，该计算框架对于在实验验证受限的拉丁美洲及其他代表性不足的人群中进行变异解读和分子诊断具有重要意义。