亮氨酸拉链螺旋间静电相互作用的系统计算机分析

《Molecules》：Systematic In Silico Analysis of Interhelical Electrostatic Interactions in Leucine Zippers

　　

亮氨酸拉链是由两条α-螺旋组成的平行卷曲螺旋（coiled-coil）二聚体，具有重复的七肽（heptad）模式，记为a–g。a和d位的疏水残基构成二聚化核心，而e和g位的带电残基产生螺旋间静电相互作用，影响稳定性与特异性。相反电荷的e–g对通常被认为具有稳定作用，同号电荷对则预期不稳定。然而，其在延伸亮氨酸拉链界面上的集体效应尚未被系统评估。研究人员设计了一个含10个e–g相互作用位点的40残基同源二聚亮氨酸拉链。10个e–g位点被独立赋予中性、成盐桥或排斥三种状态之一，产生59,049个模型。每个模型用FoldX分析相互作用能及各分量。增加盐桥数逐步改善相互作用能，相对中性对平均稳定1.68 kcal/mol。意外地，排斥对相对中性对也略有有利，改善0.48 kcal/mol。能量分解显示有利的溶剂化和范德华贡献超过了静电惩罚。这些效应具有位置依赖性，表明亮氨酸拉链稳定性反映静电、溶剂化、堆积与位置背景的平衡。

研究背景与立项依据

亮氨酸拉链（leucine zipper）作为平行卷曲螺旋（coiled-coil）二聚体，是解析蛋白质-蛋白质相互作用（PPI）界面的经典模型。其七肽重复（heptad）中a、d位疏水残基形成核心，e、g位带电残基通过盐桥（salt bridge）或同号电荷排斥调控稳定性与特异性。传统观点认定异号e–g对稳定、同号对不稳定，但多个静电相互作用在延伸界面上的集体行为、位置依赖性及非静电项（溶剂化、范德华）的相对权重长期缺乏系统定量。已有研究多聚焦单点突变或个别盐桥，难以回答“是否纯加和”“位置是否等价”“排斥对是否真如库仑定律预期般有害”。该研究以人甲状腺胚胎因子（TEF, UniProt Q10587, PDB 4U5T）亮氨酸拉链为模板，延伸为40残基同源二聚体LZ-40，借AlphaFold3建模、FoldX全组合扫描59,049（3¹⁰）个e–g状态变体，建立定量框架。论文发表于《Molecules》。

主要关键技术方法

研究人员以PDB 4U5T中TEF亮氨酸拉链（残基258–294）为起点，复制中央14残基段并末端改造获得LZ-40序列（ALAAENTALRTEVAELRKENTALRTEVAELRAEVGALRNI）；用AlphaFold3预测同源二聚体初始骨架（ipTM>0.81）；在PyMOL中定位10个e–g对（侧链距2.5–3.1 ?）；对每个e–g位赋中性（Ala–Ala, N）、盐桥（Glu–Arg/Arg–Glu, S）、同号排斥（Glu–Glu/Arg–Arg, R）三态，穷举3¹⁰组合；经FoldX 5的RepairPDB、BuildModel、AnalyseComplex以默认参数计算相互作用能（Interaction Energy）及静电、范德华、极性溶剂化、疏水性溶剂化分量；用自编Python脚本做三联比对、增量分析、位置分解与相邻耦合能统计。

研究结果

2.1 e–g静电相互作用突变库构建

研究人员将LZ-40十对e–g（E5_g–R10_e等）分别置N/S/R三态，获59,049模型。全库平均相互作用能?40.69±2.59 kcal/mol；十盐桥模板最稳（?51.01），三排斥七中性最弱（?33.74）。确认侧链几何允许盐桥形成。

2.2 盐桥与离子排斥的全局能量效应

按盐桥数分箱，相互作用能随盐桥数增而线性下降（斜率?1.489 kcal/mol，R²=0.735）；全库中排斥对数增多使能量趋正（弱相关），中性对数增多亦趋正。排除成分约束后取“仅S+N”1024模型，每增一盐桥?1.69 kcal/mol；取“仅R+N”1024模型，每增一排斥对反而?0.29 kcal/mol（R²=0.384）。相对中性基线，盐桥ΔΔG_int=?1.65±0.92，排斥ΔΔG_int=?0.33±0.51 kcal/mol，反驳“同号必损”。

2.3 盐桥与排斥的能量分解基础

以N→S比对：静电?0.56、范德华?0.54、疏溶剂?0.82（有利），极性溶剂化+0.46（去溶剂罚）；净?1.65。以N→R比对：静电+0.11（库仑罚），范德华?0.27、疏溶剂?0.62、极性溶剂化?0.15（同号相斥使带电头外露获水合利），净?0.33。表明排斥对“非静电收益>静电罚”。

2.4 位置依赖性效应

单盐桥在第三位（?3.02）、第四位（?3.17）最强，第二位仅?0.43；换骨架模型位置排序相关性低（r=0.23/?0.14）但“有利趋势”保留。排斥对呈锯齿状，ΔIE +0.20至?0.82；S优于R在所有位成立，第三位差?2.82最大、第二位?0.43最小。位置效应受骨架与侧链重建路径影响。

2.5 相邻e–g非加和耦合

将10位分奇位集A(1,3,5,7,9)与偶位集B(2,4,6,8,10)。邻位双盐桥耦合能?0.30~?0.14（协同，A均?0.238/B均?0.215），远位≈0；邻位双排斥耦合能+0.04~+0.39（抗协同，A均+0.151/B均+0.213）。耦合广分布显示背景依赖。

讨论与结论翻译

研究表明亮氨酸拉码稳定性由多重能量平衡而非单纯静电决定。同号e–g对相对Ala–Ala基线反而有利，因Glu/Arg侧链脂肪链带来疏水掩埋与范德华收益超过库仑罚；盐桥因带相反电荷近端致极性去溶剂罚，排斥对因侧链外展保水合利，缩小二者差距。但排斥对仍弱于盐桥（所有位S>R），替换盐桥为排斥会削弱界面，幅度位置依赖。邻位盐桥协同、邻位排斥抗协同，且耦合受序列背景调制。对合理设计卷曲螺旋的启示是：e–g电荷图案须与溶剂暴露、疏水核堆积、邻位静电网络及螺旋对齐一并考虑；本文定量值非普适常数，而是上下文依赖的设计原则，用于优先筛选候选突变供后续结构与实验评估。

（结论节译）亮氨酸拉链稳定性反映静电、溶剂化、堆积与位置背景之平衡；同号e–g对相对中性Ala基线具净有利性，但逊于盐桥；盐桥贡献具位置依赖与邻位协同；排斥对呈位置锯齿响应与邻位抗协同；上述规律为卷曲螺旋特异性工程提供上下文依赖框架而非单点电荷规则。