1型糖尿病中β细胞应激相关C19S胰岛素新表位对TCR–HLA-DQ8识别的结构重塑

《International Journal of Molecular Sciences》：Structural Remodeling of TCR–HLA-DQ8 Recognition by a β-Cell Stress-Associated C19S Insulin Neoepitope in Type 1 Diabetes

　　

胰腺胰岛微环境内的炎症和氧化应激（oxidative stress）可以改变胰岛素衍生肽并产生新表位（neoepitopes），这些新表位可能重塑1型糖尿病（type 1 diabetes, T1D）中的自身反应性T细胞识别。其中一种修饰，C19S，代表胰岛素B链第19位半胱氨酸到丝氨酸的替换，最近已被鉴定为人类白细胞抗原II类（human leukocyte antigen class II, HLA-II）相关胰岛素新表位，被自身反应性CD4⁺ T细胞识别。尽管C19S的生物学相关性已确定，但可能区分C19S特异性T细胞受体（T cell receptor, TCR）结合与天然胰岛素识别的分子特征仍不完全明确。在此，研究人员使用比较蛋白质-蛋白质对接、分子动力学（molecular dynamics, MD）模拟、界面接触分析、构象景观分析和结合能计算，来研究TCR对呈现天然胰岛素肽或相应C19S胰岛素肽的人类白细胞抗原DQ8（human leukocyte antigen DQ8, HLA-DQ8）的结合。初始建模表明，两种肽-HLA-DQ8复合物均与TCR结合的三元复合物形成相容。然而，预测含C19S的复合物表现出改变的肽中心动力学、肽骨架呈现的变化以及TCR-肽和TCR-HLA-DQ8接触的重组。比较分子力学泊松-玻尔兹曼表面积（molecular mechanics Poisson–Boltzmann surface area, MM/PBSA）和分子力学广义玻恩表面积（molecular mechanics generalized Born surface area, MM/GBSA）分析进一步表明，在所应用的建模条件下，含C19S的复合物具有独特的计算能量特征，残基水平分解将能量差异定位于选定的界面热点。总之，这些发现为生成关于C19S如何重塑HLA-DQ8呈现的胰岛素识别表面的假设提供了分子框架，对未来T1D中自身反应性CD4⁺ T细胞识别和抗原特异性耐受策略的实验研究具有意义。

**研究背景与问题**
1型糖尿病（type 1 diabetes, T1D）是一种自身免疫性疾病，其中自身反应性T细胞靶向由疾病相关HLA分子呈现的β细胞衍生抗原，导致产生胰岛素的胰腺β细胞发生免疫介导破坏。在HLA II类分子中，HLA-DQ8与T1D易感性强相关，并参与自身反应性抗原识别。尽管天然胰岛素是T1D中的核心自身抗原，但决定选定胰岛素肽在HLA-DQ8上持续存在及其被自身反应性T细胞受体（T cell receptor, TCR）识别的结构决定因素仍不完全明确。最近免疫肽组学和功能研究鉴定出胰岛素B链第19位半胱氨酸到丝氨酸的转化，称为C19S，作为自身免疫性糖尿病中微环境相关的胰岛素新表位。C19S胰岛素肽（包括InsB9-23 (C19S)）在HLA-II肽组中被检测到，并在氧化应激、炎症细胞因子暴露和β细胞功能障碍条件下生成。该修饰与单残基水平的胰岛素自身反应性改变以及被特定结合寄存器的CD4⁺ T细胞识别相关，包括人类疾病背景下的HLA-DQ8限制性反应。然而，C19S相关肽呈现变化如何延伸至TCR面对界面的分子特征尚未完全阐明。因此，需要分析完整的TCR-HLA-DQ8-胰岛素肽三元复合物，以定义C19S如何改变TCR-可及的识别表面。

**研究内容与结论**
研究人员扩展了先前对肽-HLA-DQ8复合物的分析，通过比较蛋白质-蛋白质对接、分子动力学（molecular dynamics, MD）模拟、界面接触分析、构象景观分析和结合能计算，对HLA-DQ8呈现天然胰岛素肽或C19S胰岛素肽的完整TCR-HLA-DQ8-胰岛素肽三元复合物进行了系统比较。研究得出以下结论：C19S复合物在整体三元结构上与天然肽复合物相容，但表现出更有利的计算结合自由能谱；C19S改变了肽中心动力学，导致肽骨架呈现和界面接触网络的重组；C19S与TCR-肽和TCR-HLA-DQ8接触的重新分布相关，并伴有极性相互作用和静电相互作用的预测变化；晚期构象景观分析表明，C19S复合物占据了与天然肽复合物可区分的构象集合。这些发现为C19S如何重塑HLA-DQ8呈现的胰岛素识别表面提供了分子框架，对自身反应性CD4⁺ T细胞识别和T1D中抗原特异性耐受策略的未来实验研究具有重要意义。该论文发表在《International Journal of Molecular Sciences》。

**主要关键技术方法**
研究人员使用以下主要技术方法（基于PDB ID: 6XCP和1JK8的晶体结构模板，无具体样本队列来源）：
1. 蛋白质-蛋白质对接（LightDock框架）：用于模拟TCR与HLA-DQ8-肽复合物的结合取向。
2. 全原子分子动力学（MD）模拟（GROMACS，显式SPC水溶剂，生理离子强度0.15 M NaCl，500 ns生产模拟）：用于比较WT和C19S三元复合物的动态行为。
3. 界面接触、氢键和盐桥分析：计算轨迹中的接触占有率、氢键持续性和盐桥持久性，以比较TCR-肽和TCR-HLA-DQ8相互作用模式。
4. 构象景观分析（主成分分析PCA、聚类和自由能景观）：用于评估晚期模拟窗口中WT和C19S三元复合物的构象集合差异。
5. 结合能计算（MM/PBSA和MM/GBSA，gmx_MMPBSA包）：用于比较TCR与WT或C19S肽-HLA-DQ8复合物结合的计算相对能量谱，并进行残基水平能量分解。

**研究结果**
2.1. C19S保留与TCR相容的HLA-DQ8呈现（相对于WT胰岛素肽）
通过蛋白质-蛋白质对接和MD模拟，研究人员发现天然肽（WT）和C19S肽均能形成稳定的TCR-HLA-DQ8-肽三元复合物结构。C19S复合物在模拟前后均保持了整体三元架构，表明半胱氨酸到丝氨酸的替换未破坏TCR-HLA-DQ8-肽的整体组装，为后续比较提供了基础。

2.2. C19S与更有利的计算结合自由能谱相关
MM/PBSA分析表明，C19S复合物显示出比WT复合物更有利的计算结合自由能估算。成分水平能量分解显示，C19S替换与范德华、静电、极性溶剂化和非极性溶剂化贡献的重新分布相关，而非均匀变化。MM/GBSA分析进一步验证了这一趋势，残基水平分解将能量差异定位于选定界面热点。

2.3. C19S改变TCR-HLA-DQ8复合物中的肽中心动力学
骨架RMSD分析显示，C19S复合物在模拟中表现出与WT复合物不同的轨迹行为。残基水平RMSF分析表明，最明显的扰动集中在胰岛素肽上，并延伸至TCR-HLA-DQ8组装的受体面对区域。晚期肽构象比较显示，C19S肽在多个位置的骨架呈现与WT不同。

2.4. C19S与TCR面对接触架构的改变相关
接触占有率分析显示，C19S复合物中TCR-肽界面的接触模式与WT显著不同，暗示替换肽通过不同的局部接触网络被结合。TCR-HLA接触图也显示差异，表明C19S效应扩展至更广泛的HLA-DQ8受体面对表面，导致协调的界面重组。

2.5. C19S与极性相互作用和静电相互作用的预测重新分布相关
TCR-肽氢键分析显示，C19S复合物中氢键模式重新分布。盐桥分析进一步区分了C19S界面，发现一组重构的静电相互作用，其中TCR链A Arg52-肽Glu374盐桥在C19S轨迹中突出。残基水平能量分解独立突出了同一界面区域，与氢键和盐桥分析一致。

2.6. C19S与晚期模拟窗口中的构象分离相关
联合主成分分析（PCA）将WT和C19S晚期构象集合沿主要构象坐标分离。PCA空间聚类表明两种复合物占据不同构象区域。自由能景观分析进一步显示，WT和C19S在PCA空间中占据不同的低能量盆地，表明C19S复合物采样的构象集合与WT复合物可区分。

**讨论与结论**
讨论部分指出，新表位形成是连接β细胞应激、抗原呈现改变和T1D中免疫耐受失败的机制。C19S胰岛素肽在人类和小鼠HLA-II肽组中被发现，并由自身反应性CD4⁺ T细胞识别，提示该修饰参与疾病相关抗原呈现途径。本研究通过计算建模提供了TCR-可及的HLA-DQ8识别表面的结构框架。结果支持一个模型，即C19S与TCR-HLA-DQ8-胰岛素肽界面的协调重塑相关，包括肽呈现、受体面对接触和界面能量学的改变。这些发现与MHC-II抗原性受肽序列、结合寄存器、肽-MHC稳定性、HLA-DM编辑和TCR-可及表面几何形状影响的观点一致。C19S可能通过改变肽表面氢键势、重新分布局部静电相互作用并微妙地改变界面堆积，同时维持整体结合几何形状，来影响TCR界面。横向比较分析的一致性表明，C19S效应通过协调的界面重塑而非孤立数值变化产生。研究局限性包括：基于计算模型，仅代表一种TCR-HLA-DQ8-肽背景，未包含肽生成、内体加工、HLA-DM编辑、膜定位、CD4共受体或下游信号。未来需通过高分辨率结构、生物物理测量和功能实验验证预测。

**研究结论部分翻译**：
总体而言，本研究提供了一个分子建模框架，将C19S胰岛素新表位呈现与TCR-HLA-DQ8界面架构的预测改变联系起来。结果表明，C19S可能通过肽呈现、受体面对接触、极性和静电相互作用、构象采样以及能量热点分布的协调变化，影响自身反应性识别的结构决定因素。这些计算观察结果补充了近期关于T1D中C19S特异性CD4⁺ T细胞应答的免疫肽组学和功能证据，并通过纳入模拟的TCR结合识别界面，扩展了先前的肽-HLA-DQ8模拟。所得框架为自身免疫性糖尿病中C19S特异性自身反应性识别的结构、生化和功能验证提供了可检验的假设。