《ACS Omega》：Microscopic Insights into a Ligand Escape Pathway and Energetics in Leucine–Isoleucine–Valine Binding Protein
编辑推荐：
本文针对以往对LIV结合蛋白（LIVBP）构象波动关注较多、但对配体逃逸的微观机制与能量学细节研究不足的问题，研究人员通过μs级原子分子动力学与well-tempered元动力学模拟，揭示了异亮氨酸在LIVBP构象选择中的作用，并首次从微观层面阐明了配体逃逸的路径、能量学与事件序列。结果表明，配体逃逸为两步过程（域分离先于配体逃逸，但反向过程亦存在），且水分子通过在两域间形成桥接氢键，在配体解离以及无配体时稳定闭合构象中起关键作用。该研究为理解PBPs的配体逃逸机制提供了新的能量与结构视角，对后续针对分支链氨基酸运输的药物设计具有启示意义。
在革兰氏阴性细菌中，周质结合蛋白（Periplasmic binding proteins, PBPs）是一个庞大的受体与转运蛋白家族，它们在细胞运输、化学趋向性和群体感应等关键生理过程中扮演着核心角色。这类蛋白具有独特的双域结构，能够通过铰链运动发生大幅度的构象转变，其工作模式常被比喻为“捕蝇草”。其中，亮氨酸-异亮氨酸-缬氨酸结合蛋白（Leucine–Isoleucine–Valine binding protein, LIVBP）专一性感知分支链氨基酸，并经历从开放构象到闭合构象的转变，这一过程传统上被视为配体诱导的构象变化。尽管已有大量研究聚焦于LIVBP的构象波动，但配体如何从结合口袋中“逃逸”出去的微观细节、能量变化及事件顺序，却如同一个尚未被充分探索的“黑箱”。配体的解离路径是怎样的？是蛋白质先“张开嘴”（域分离），配体再“跳出来”，还是相反？在没有配体的情况下，蛋白质的闭合构象为何能保持稳定？这些问题的答案，对于从原子层面理解PBPs的工作机制乃至针对相关通路的药物设计都至关重要。
为了回答这些问题，研究人员在《ACS Omega》上发表了一项研究，他们综合运用了长时间尺度的全原子分子动力学（Molecular Dynamics, MD）模拟和增强采样技术——well-tempered元动力学（metadynamics）模拟，对LIVBP与配体异亮氨酸（Isoleucine, Ile）的相互作用进行了深入探究。研究人员从蛋白质数据库（PDB ID: 1Z17）中获取了配体（异亮氨酸）结合形式的初始原子坐标，并利用CHARMM-GUI服务器生成GROMACS模拟所需文件。系统被置于充满TIP3P水分子和0.15 M Na+和Cl–离子的立方盒子中，采用CHARMM36m力场进行模拟。在完成能量最小化和平衡后，他们进行了两次独立的微秒级无偏分子动力学模拟。为了探究配体逃逸的自由能垒，研究进行了五组独立的单变量元动力学模拟，选取的集体变量是Glu-328的Cα原子与配体质心之间的距离。此外，为了理解域开启与配体逃逸之间的协同作用，他们还进行了一项双变量元动力学模拟，同时以蛋白质-配体距离和域间质心距离作为集体变量，从而构建了二维自由能景观。分析工作主要使用GROMACS后处理工具和Python脚本完成，可视化则借助VMD软件。
结果
1. 从无偏分子动力学轨迹观察构象波动与配体逃逸
在长达1微秒的无偏MD模拟中，研究人员监测了蛋白质-配体距离、域间距离和域间角度的变化。在一条轨迹中，他们观察到域分离发生在配体逃逸之前约40纳秒。在域分离后、配体逃逸前的这段时间里，配体依然通过氢键与结构域2（Domain-2）上的氨基酸（如Ser-79和Thr-102）结合，而与其在结构域1（Domain-1）上的结合伙伴（如Tyr-202）的相互作用则被破坏。有趣的是，在该轨迹中，配体在逃逸后约500纳秒时曾重新进入结合裂隙并与部分氨基酸重新形成氢键，但这次再结合并未能诱导蛋白质重新闭合，蛋白质仍保持开放状态，配体最终再次逃逸。这些观察表明，配体与结构域2的相互作用可能更为关键，并且配体结合并不总是能触发“捕蝇草”式的闭合运动。
2. 配体逃逸的两条微观路径
通过元动力学模拟，研究揭示了配体逃逸的两条主要路径。在路径I中，结构域首先发生分离（即蛋白质“先开口”），随后配体才从结合裂隙中解离。在此路径的中间态，配体主要与结构域2结合，并被部分水化。在路径II中，结构域仅轻微开启，配体的完全逃逸与结构域的全开过程近乎同步发生。在此路径的中间态，配体仍与三个（S79, T102, E226）而非全部的四个关键氨基酸侧链保持氢键相互作用，并且通过一个水分子介导的桥接氢键与Tyr-202相连。两条路径的最终阶段，配体都完全被水分子溶剂化，从而完成逃逸。
3. 二维自由能景观揭示的稳定态与转变路径
通过双变量元动力学模拟构建的二维自由能景观，清晰地展示了几個稳定的能量洼地。状态1是配体结合的本土闭合构象。状态2对应一个开放构象，但配体仍附着在结构域2上，这与无偏模拟中观察到的中间态一致，也是从状态1到配体完全解离的主要能量路径上的关键中间体。状态4是配体完全解离、蛋白质呈开放构象的状态。特别值得注意的是状态5，它代表了无配体的闭合构象，挑战了“配体是维持闭合构象所必需”的传统观点。自由能景观表明，从状态1（闭合结合态）到状态2（开放-配体部分结合态）的转变，即域分离先于配体逃逸的路径，是主导的逃逸途径。
4. 水分子在稳定无配体闭合构象中的关键作用
为了探究无配体时闭合构象如何保持稳定，研究人员比较了超开放、配体结合闭合以及无配体闭合三种构象的模拟轨迹。径向分布函数分析显示，在无配体的闭合构象中，结合裂隙区域存在大量水分子，而在配体结合的闭合构象中，该区域相对“干燥”。进一步的微观结构分析发现，在无配体情况下，多个水分子会渗透到裂隙中，在原本由配体介导的两个结构域的关键极性氨基酸（如Glu-226, Tyr-202, Ser-79, Thr-102）之间形成复杂的、扩展的水介导桥接氢键网络。这些水桥如同“分子胶水”，将两个结构域锚定在一起，从而稳定了闭合构象。氢键时间关联函数和氢键数量分布的分析量化了这一作用：带负电的Glu-226与水形成的氢键数量最多、寿命最长，其他几个残基也形成了寿命显著长于本体水中氢键的稳定相互作用，共同维持了无配体闭合构象的稳定。
结论与讨论
本研究通过综合计算模拟手段，深入揭示了LIV结合蛋白的构象波动及其与配体相互作用的微观机制，特别是在配体逃逸动力学方面取得了新的认识。主要结论如下：首先，配体逃逸是一个多步骤过程，存在至少两条微观路径。主导路径（路径I）是域分离先于配体完全逃逸，这与传统“诱导契合”模型的顺序预期相符，但研究也观察到配体逃逸与域开启近乎同步的路径（路径II），表明该过程具有一定灵活性。其次，研究挑战了“配体是维持蛋白质闭合构象所必需”的经典观点。模拟结果表明，即使没有配体，闭合构象也能通过水分子在双域间形成的、稳定的桥接氢键网络而保持稳定。反之，配体的存在也并非总能将蛋白质“锁”在闭合状态，无偏模拟中观察到的自发域分离证实了这一点。这凸显了水在蛋白质构象稳定中的主动作用，而不仅仅是作为被动的溶剂。
这些发现从能量和原子细节层面，完善了对PBP家族蛋白，特别是LIVBP，其功能循环的理解。研究所识别出的两条配体逃逸路径、与结构域2结合的中间态、配体在蛋白质外围的亚稳态位点，以及由水桥稳定的无配体闭合态，共同描绘了一幅复杂的、多态的能量景观图。这幅图谱表明，配体的结合与解离并非简单的“开-关”切换，而是涉及一系列构象和溶剂化中间态的连续动态过程。
这项工作的重要意义在于，它不仅回答了关于LIVBP配体逃逸路径和能量学的具体科学问题，其揭示的机制特征——特别是水介导的稳定作用和无配体闭合态的可行性——为未来的药物设计提供了新的思路。例如，可以设想设计小分子来特异性稳定某条逃逸路径上的非生产性中间态，或扰动稳定闭合构象的结构化水网络，从而阻断分支链氨基酸的运输，这可能在针对依赖此类氨基酸的细菌生长抑制策略中具有应用潜力。尽管该研究通过模拟获得了自由能垒，但要将其转化为实际的动力学速率常数（如koff）仍面临挑战，因为能垒较宽且摩擦效应显著。未来的研究可以引入更复杂的集体变量（如结合口袋的水化状态），并利用恒定pH分子动力学等方法，进一步探究局部静电环境波动对已识别逃逸路径相对稳定性的影响，并将此框架拓展至其他PBP家族成员，以检验所发现机制的普适性。
