基于NMR光谱的透明质酸溶液构象分析：链长依赖性构象差异及β糖苷键构象平衡研究

《Carbohydrate Polymer Technologies and Applications》：Solution conformational analysis of hyaluronan with different chain lengths by NMR spectroscopy

本研究针对透明质酸(HA)溶液构象研究中存在的争议，通过NMR光谱技术系统比较了八糖(HA8)和二十糖(HA20)的构象特征。研究发现β(1→3)糖苷键呈现单一主导构象，而β(1→4)糖苷键存在构象平衡，首次为MD模拟预测的次要反式构象提供了实验证据。该研究为不同链长糖链构象分析提供了方法学参考，对理解HA生物功能具有重要意义。

　　

透明质酸（Hyaluronic Acid, HA）是一种广泛存在于脊椎动物细胞外基质中的重要糖胺聚糖（Glycosaminoglycan, GAG），由D-葡萄糖醛酸（GlcA）和N-乙酰-D-葡萄糖胺（GlcNAc）通过交替的β(1→3)和β(1→4)糖苷键连接而成。尽管其一级结构相对简单，但HA在溶液中的高级结构和动态行为却异常复杂，至今尚未被完全理解。这种复杂性主要源于其糖苷键的构象灵活性，以及环境因素（如水分、阳离子）对其构象的影响。多年来，关于HA溶液构象的研究报告了看似矛盾的证据，特别是关于两种糖苷键的结构和动态行为。此外，以往的研究多集中于较短的HA寡糖（如八糖，HA₈）的核磁共振（NMR）分析，而针对短链与长链HA构象比较的研究则相对匮乏。因此，深入探究不同链长HA的溶液构象，对于准确理解其生理功能和作用机制至关重要。本研究旨在通过先进的NMR技术，详细比较HA八糖（HA₈）和二十糖（HA₂₀）的构象特征，重点关注糖苷键和N-乙酰基的构象，以期揭示链长对HA构象的影响，并为分子动力学（MD）模拟的预测提供关键的实验验证。相关研究成果发表在《Carbohydrate Polymer Technologies and Applications》上。

为开展此项研究，研究人员主要运用了以下关键技术方法：首先，通过酶解（睾丸透明质酸酶）和高效阴离子交换色谱（SAX-HPLC）结合尺寸排阻色谱，从400 kDa的HA多糖中分离并纯化了结构明确的HA₈和HA₂₀样品。其次，利用高场核磁共振波谱仪（600 MHz），在严格控制条件（pH 7.0, 5°C和25°C）下，对样品进行了一系列二维NMR实验。核心实验技术包括：1H,13C-HSQC-NOESY（异核单量子相干-核奥弗豪泽效应光谱）系列实验，通过混合时间（τ_m）变化（20-300 ms）并结合PANIC（峰值振幅归一化改进交叉弛豫）方法，精确测量了关键的原子间质子距离；1H,13C-HSQC-ROESY（旋转框架奥弗豪泽效应光谱）实验，用于鉴别和排除自旋扩散的影响，确保NOE距离测量的可靠性；以及J-HMBC（异核多键相关）和恒时J-HMBC实验，用于定量测量跨越糖苷键和N-乙酰基内部的三键碳氢耦合常数（³J_CH）。最后，将实验测得的质子距离和耦合常数与基于已知模型结构（如PDB 2BVK）和参数化Karplus方程计算的理论值进行比较，通过最小化均方根偏差（RMSD）来确定最可能的二面角构象，并对可能存在构象平衡的情况进行了群体加权分析。

3.1. HA样品的分离

研究首先通过睾丸透明质酸酶降解高分子量HA（400 kDa），然后利用制备型阴离子交换色谱（SAX-HPLC）成功分离出具有偶数个单糖残基的HA₈和HA₂₀。经过脱盐和离子交换处理后，获得高纯度的钠盐形式样品，用于后续NMR分析。

3.2. 原子间质子距离和J-耦合常数的测量

为了克服传统1H,1H-NOESY谱图中信号重叠严重的问题，本研究采用了1H,13C-HSQC-NOESY技术。该技术利用碳维的化学位移分散度，使得关键的残基间NOE信号（如G₁-N₃, G₁-NH, N₁-G₄等）能够清晰分辨，从而实现了无重叠干扰的原子间质子距离的精确测量。通过分析不同混合时间下的NOE信号强度构建曲线，并结合HSQC-ROESY谱图排除自旋扩散的干扰，研究人员获得了用于构象分析的高质量距离约束。同时，通过J-HMBC实验测量了与糖苷键二面角（φ, ψ）和N-乙酰基扭转角（θ₁）相关的³J_CH耦合常数。

3.3. N-乙酰基的构象

通过对NH-N₁、NH-N₂、NH-N₃质子距离以及³J_H2-NH和³J_H2-CO耦合常数的分析，研究发现HA₈和HA₂₀中GlcNAc残基的N-乙酰基均主要采取反式（anti）构象（θ₁接近180°）。HA₈和HA₂₀的全局最小值分别位于(165±6)°和(-167±13)°，表明链长对N-乙酰基的构象没有显著影响，其主导构象与之前的MD模拟和NMR研究结果一致。

3.4. β(1→3)糖苷键的构象

通过分析G₁-N₃、G₁-N₂、G₁-NH等质子距离以及³J_N(C3)-G(H1)耦合常数，并结合模型计算，研究人员确定了β(1→3)糖苷键的构象。结果表明，对于HA₈和HA₂₀，该键均呈现单一的主导顺式（syn）构象。HA₈的全局最小二面角为(φ₁₃, ψ₁₃) = ((46±6)°, (15±9)°)，而HA₂₀为((51±7)°, (35±18)°)。这些数值与多数MD模拟预测的范围（φ₁₃ ~46-52°, ψ₁₃ ~ -10° 至 25°）吻合良好，首次为MD模拟建立的β(1→3)糖苷键主导构象提供了坚实的实验支持。分析未发现次要构象存在的显著证据。

3.5. β(1→4)糖苷键的构象

对β(1→4)糖苷键的分析（基于N₁-G₄、N₁-G₃、N₁-G₅距离和³J_G(C4)-N(H1)耦合常数）显示，实验数据无法用单一构象模型很好地解释。特别是N₁-G₃和N₁-G₅的实测距离显著短于基于单一顺式构象模型的计算值。这强烈暗示该糖苷键存在构象平衡。研究人员进而引入了近期MD模拟预测的一个主导顺式构象（Conformation I: (φ₁₄, ψ₁₄) ≈ (50°, -3°)，ΔG=0 kcal/mol）和两个次要反式（anti）构象（Conformation II: (35°, 167°), ΔG=1.31 kcal/mol; Conformation II‘: (60°, -153°), ΔG=1.42 kcal/mol）进行加权计算。发现当HA₂₀中构象I、II、II’的占比分别为83%、9%、8%时，计算出的加权平均距离和耦合常数与实验数据高度吻合。对于HA₈，最佳拟合比例为构象I占92%，构象II‘占8%（构象II占比为0）。这是首次通过NMR实验证实了β(1→4)糖苷键上存在次要反式构象，为MD模拟的预测提供了关键实验证据。这些反式构象被认为与钠离子诱导形成的发夹状（hairpin-like）结构有关。

3.6. HA₈与HA₂₀构象比较

总体而言，HA₈和HA₂₀在N-乙酰基和糖苷键的构象上差异较小。主要区别在于：HA₂₀的β(1→3)糖苷键的ψ角比HA₈大约20°；更重要的是，HA₂₀中β(1→4)糖苷键的反式构象总占比（17%）高于HA₈（8%）。这意味着平均每个HA₂₀分子中约有1.5个β(1→4)键处于反式构象，而HA₈中约为0.24个。这种构象差异可能导致长链HA形成更多局部的发夹状转折，使其溶液结构不如基于短链模型预测的那样伸展。

本研究通过高分辨NMR光谱技术，首次在实验上证实了HA溶液中β(1→4)糖苷键存在构象平衡，包含一个主导的顺式构象和两个次要的反式构象，并精确量化了它们在HA₈和HA₂₀中的分布比例。同时，研究也为β(1→3)糖苷键的主导顺式构象提供了坚实的实验支持。研究表明，链长对HA的局部构象影响相对有限，但长链HA中β(1→4)键反式构象比例的升高可能预示着其整体构象更具柔性和复杂性，更容易形成发夹状结构。这一发现对于理解HA在不同链长下的生物物理性质（如链刚度、流体力学体积）及其与蛋白质的相互作用具有重要意义。本研究建立的基于HSQC-NOESY和J-HMBC的构象分析方法，为研究其他糖胺聚糖乃至复杂多糖的溶液构象，特别是构象平衡问题，提供了可靠的方法学范例。未来研究可进一步探讨不同离子种类、浓度以及更长的HA链对其构象平衡的影响，从而更全面地揭示HA结构与功能之间的关系。