《Results in Chemistry》：Hydrophobic chain-length effects on the binding of
N-acyl homoserine lactones with cucurbit[7]uril: Molecular dynamics and free-energy analysis
编辑推荐：
研究人员采用综合分子动力学（MD）模拟与MM/GBSA方法分析了葫芦[7]脲（CB7）与N-酰基高丝氨酸内酯（AHL）的主客体相互作用。结果表明，AHL对CB7的结合亲和力随疏水性酰基链长度的增加而逐步增强。N-辛酰基-L-高丝氨酸内酯（C8）复合物表现出最高
研究人员采用综合分子动力学（MD）模拟与MM/GBSA方法分析了葫芦[7]脲（CB7）与N-酰基高丝氨酸内酯（AHL）的主客体相互作用。结果表明，AHL对CB7的结合亲和力随疏水性酰基链长度的增加而逐步增强。N-辛酰基-L-高丝氨酸内酯（C8）复合物表现出最高的热力学稳定性和最强的结合自由能。范德华（vdW）相互作用与CB7空腔内的疏水堆积被鉴定为复合物形成的主导驱动力。较长的烷基链降低了溶剂化惩罚、改善了空腔占据效率并增强了超分子组装体的构象稳定性。
研究背景方面，超分子化学通过非共价相互作用调控分子识别与自组装，葫芦[n]脲（CB[n]s）因其结构刚性、热稳定性及独特空腔识别特性而占据重要地位。CB7作为研究最广泛的同系物之一，具有最佳空腔尺寸和水溶性，对多种生物相关底物表现出强主客体络合行为。N-酰基-L-高丝氨酸内酯（AHLs）是细菌群体感应（quorum sensing）信号分子，调控生物膜形成、毒力、运动性及致病性等集体微生物行为。AHL的疏水酰基链影响分子识别与膜转运行为，使其成为研究链长依赖性超分子包封的理想模型体系。理解AHL–CB7相互作用的分子基础对于开发通过选择性客体识别和螯合来调控细菌通讯的超分子体系具有重要意义。然而，实验表征通常难以提供构象动力学、能量贡献及分子间相互作用机制的原子级信息。为此，研究人员开展了CB7与烷基链长从戊酰基到辛酰基的系统同系AHL系列的计算研究，旨在阐明渐进性疏水链延伸如何影响超分子结合亲和力、构象稳定性及CB7空腔内分子间相互作用模式，重点关注疏水堆积效率、去溶剂化行为与结构动力学同热力学趋势的关联。
研究人员采用的关键技术方法包括：分子动力学（MD）模拟，使用AMBER18软件包，采用GAFF力场和Espaloma机器学习电荷模型，在TIP3P水盒中进行200 ns生产模拟；MM/GBSA结合自由能计算，从200 ns轨迹中均匀提取2000个快照进行端点自由能分析，并通过块平均法估计统计不确定性；空腔占据分析，以CB7空腔碳原子质心为参考，统计5.0 ?范围内酰基链碳原子数；质心径向分布函数（COM-RDF）分析，评估主客体结构稳定性；以及RMSD和能量轨迹分析评估体系稳定性。样本队列来源为经修饰酸催化缩合 procedure 合成的CB7，AHL分子结构使用ChemDraw Professional 17.0生成，初始主客体复合物使用UCSF Chimera构建。
研究结果方面，在"研究主客体体系的结构特征"部分，通过分子可视化发现AHL烷基链优先朝向CB7疏水空腔内部，而极性内酯和酰胺基团靠近羰基修饰的端口；结构兼容性随酰基链长增加逐步改善，较长链类似物表现出更高效的空腔填充和增强的分子间堆积相互作用。在"N-酰基-L-高丝氨酸内酯–CB7复合物"部分，通过MM/GBSA计算得出结合自由能顺序为C8 > C7 > C6 ≈ C5，C8–CB7复合物表现出最有利的结合自由能（ΔGbind = ?29.53 ± 1.15 kcal/mol）；复合过程主要由有利的范德华相互作用驱动（约?28至?30 kcal/mol），静电贡献较小；溶剂化贡献从C5的+2.37 ± 0.42 kcal/mol显著降至C8的+0.39 ± 0.35 kcal/mol。在"MM/GBSA结合自由能的收敛性"部分，通过将200 ns轨迹分为五个重叠时间窗口重新计算，发现各复合物ΔGbind绝对值在所有窗口中稳定，相对排名在每个窗口中均保持，证明链长依赖性趋势稳健且非采样不足所致。在"与LasR–AHL识别系统及实验CB7数据的比较"部分，发现CB7–AHL包封模式与LasR天然受体结合架构平行，两者均以疏水酰基链包封为主要识别要素，但CB7链长选择性主要由空腔占据和疏水堆积效率决定，而LasR表现出更广泛的配体混杂性；计算结果与实验观察到的CB7以链长依赖方式识别AHL定性一致。在"空腔占据分析"部分，通过统计空腔内酰基链碳原子数，发现空腔占据从C5的3.2 ± 0.4单调增至C8的7.9 ± 0.2，C8酰基链几乎完全填充CB7空腔；空腔占据与ΔGsolv的线性回归R2 = 0.94，与ΔGbind的R2 = 0.91。在"体系稳定性"部分，通过RMSD分析发现C5–CB7的RMSD最高（4.68 ± 0.58 ?），C8–CB7最低（1.27 ± 0.23 ?），表明较长链衍生物形成更构象刚性的主客体组装体；能量轨迹保持稳定无显著漂移。在"主客体结构稳定性：COM-RDF分析"部分，通过质心径向分布函数分析发现C5–CB7呈现宽双峰分布指示部分脱出构型，C7–CB7和C8–CB7呈现尖锐单峰确认完全包封；COM标准差从C5的0.96 ?单调降至C8的0.48 ?。在"主客体化学计量比"部分，指出本研究仅考虑1:1复合物，CB7空腔可容纳两个客体分子形成1:2化学计量比，但在所模拟的稀释1:1条件下1:1结合预期占主导。
讨论部分总结了本研究的局限性：研究采用经典MD和MM/GBSA端点方法，未产生沿结合坐标的完整自由能剖面；模拟从预形成包封复合物出发而非分离的主客体分子；每个复合物仅使用单一手动构建的起始几何构型；未进行伞形采样或PMF计算，因此结果描述的是采样构型内的相对热力学稳定性而非全局自由能最小值或结合/解离态间能垒；C5–CB7的COM-RDF部分客体脱出解释来自分布展宽和轨迹检查而非直接计算的自由能剖面；短链同系物间MM/GBSA结合自由能的统计不确定性与差异相当，C5至C8亲和力增加趋势稳健（尤其C8的更强结合），但C5、C6、C7间排名落在统计不确定性内不应过度解读。未来工作可通过伞形采样或well-tempered metadynamics沿主客体质心距离计算PMF，揭示客体进入的自由能能垒、替代结合模式的稳定性及链长依赖性趋势的起源；从体相水中分离的主客体出发的无偏模拟可检验包封构型是否在平衡态自发达到。
研究结论部分翻译如下：总之，本计算研究提供了葫芦[7]脲（CB7）与疏水酰基链长不同的同系N-酰基-L-高丝氨酸内酯（AHLs）系列之间主客体相互作用的详细分子级理解。分子动力学模拟结合MM/GBSA自由能计算揭示了从戊酰基（C5）到辛酰基（C8）衍生物的结合亲和力和结构稳定性的明确链长依赖性增强。能量分解分析表明复合物形成主要由烷基链在CB7空腔内高效疏水堆积产生的有利范德华相互作用驱动，而静电贡献相对较小。随链长增加溶剂化惩罚的逐步降低表明包封时疏水表面免受水环境屏蔽的改善。与此一致，COM-RDF分析显示从C5到C8主客体分布逐步锐化和增强，COM标准差从C5–CB7的0.96 ?单调降至C8–CB7的0.48 ?，表明较长酰基链的包封越来越紧密。RMSD轨迹进一步证实较长链衍生物在模拟期间形成更构象刚性的主客体组装体。在本研究采样的四个复合物中，C8–CB7体系表现出最有利的MM/GBSA结合自由能、最低的RMSD和最窄的COM分布，因此被鉴定为采样构型中最稳定的复合物。研究人员强调该结论来自经典MD和端点MM/GBSA计算，反映所探索构型内的相对稳定性而非严格确定的全局自由能最小值。总之，本研究结果建立了CB7内疏水链延伸与超分子结合效率之间的直接相关性，并提供了调控CB7基主客体识别的分子因素的机制洞察。这些发现凸显了葫芦脲作为选择性识别和调控生物相关群体感应分子的超分子平台的潜力，同时将增强采样方法（如伞形采样、well-tempered metadynamics和无偏结合模拟）确定为未来工作的自然方向。
