四个手性异构体铕配合物与人血清白蛋白结合选择性的分子机制

《Physical Chemistry Chemical Physics》：Molecular mechanism of binding selectivity in four stereoisomeric europium complexes with human serum albuminOpen Access

　　

一系列立体异构的铕配合物在与主要人体转运蛋白——人血清白蛋白（HSA）结合时表现出立体选择性。研究人员利用分子建模计算，深入研究了两种手性Eu(III)配合物对与HSA中药物位点2（drug site 2, DS2）之间结合相互作用的本质。在全原子分子动力学（molecular dynamics, MD）模拟中进行了集合对接和精修，从而能够评估主要的结合构象。该方法可与实验数据进行比较，并能合理解释相对结合自由能，以及通过选择性范德华（van der Waals, vdW）作用、氢键（hydrogen bonding）和与药物位点2中邻近氨基酸残基（尤其是Arg-410和Lys-414）的阳离子-π（cation–π）相互作用，追溯结合亲和力差异的起源。研究揭示了每种配合物与HSA的结合如何影响蛋白质和Eu(III)配合物的构象偏好，为所观察到的选择性起源提供了深入见解。

论文解读：四种立体异构铕配合物与人血清白蛋白结合选择性的分子机制

一、研究背景与科学问题

人血清白蛋白（human serum albumin, HSA）是血浆中含量最丰富的蛋白质，在调控生理活性化合物的转运与分布中发挥核心作用。HSA由三个同源结构域组成，其中包含多个结合口袋，包括药物位点1（drug site I, DS1）和药物位点2（drug site II, DS2）。鉴于HSA在药物药效学中的重要性，开发用于检测蛋白质结合事件的发光探针一直是研究热点。尽管在表征配体与HSA结合方面取得了显著进展，但对探针在其多个结合口袋中选择性的全面理解仍不完整。具体而言，探针对DS1、DS2及次级区域的差异偏好尚未明确解析，控制这种选择性的结构特征也未得到系统评估。因此，现有探针能否可靠地区分这些结合环境、准确报告位点特异性相互作用，仍是未解之谜。

已有研究发现，向手性Eu(^III)配合物Δ-[EuL^1a]³⁺和Δ-[EuL^1b]³⁺中加入HSA后，其主要的方形反棱柱（square-antiprismatic, SAP）构象转变为扭曲的方形反棱柱（twisted square-antiprismatic, TSAP）构象，而对映体配合物Λ-[EuL^2a]³⁺和Λ-[EuL^2b]³⁺则未发生此变化。为了阐明这些手性非共价Eu配合物与HSA结合亲和力差异的起源，本研究对SAP和TSAP两种构象的配合物与HSA的相互作用进行了建模分析。

二、关键技术方法

研究人员采用了一系列计算化学方法：首先，利用密度泛函理论（density functional theory, DFT）优化Eu(^III)配合物及参照配体布洛芬（ibuprofen）的几何结构并分配部分电荷；其次，通过全原子分子动力学（MD）模拟对无配体的HSA受体构象进行采样，使用Gromos聚类算法提取代表性构象构成受体集合；再次，运用GNINA软件将配体对接至该受体集合，筛选初始结合模式；随后，对HSA-配体复合物进行多副本MD模拟；最后，采用分子力学泊松-玻尔兹曼表面积（molecular mechanics Poisson–Boltzmann surface area, MM/PBSA）方法计算结合自由能，并通过主成分分析（principal component analysis, PCA）和均方根偏差（root-mean-square deviation, RMSD）分析表征蛋白质构象动态与配体柔性变化。

三、主要研究结果

1. 竞争实验与圆偏振发光光谱表明配合物的HSA结合位点

研究人员重新测定了各Eu配合物与HSA的结合常数（log?K），采用比率分析法处理原始数据。HSA药物位点竞争实验表明，强DS1结合剂碘帕胺（iodipamide）无法有效置换结合态Eu配合物，而DS2结合剂布洛芬可导致发射强度下降约18%-20%，且信号未完全恢复。这表明Eu配合物优先结合于DS2，但可能不限于单一药物位点，或布洛芬与Eu配合物可在DS2附近共结合。圆偏振发光（circularly polarised luminescence, CPL）光谱显示，Δ异构体结合HSA后谱带发生明显位移，而Λ异构体则无此现象，进一步证实了蛋白质结合对Δ异构体构象的调控作用。

2. 计算相对结合自由能再现实验趋势

MM/PBSA计算的结合自由能趋势与实验结果基本一致。布洛芬结合最强，Δ-[EuL^1a]³⁺、Δ-[EuL^1b]³⁺和Λ-[EuL^2b]³⁺形成一组亲和力略低但相近的配合物，Λ-[EuL^2a]³⁺结合最弱。对于Δ-[EuL^1a/1b]³⁺，TSAP构象的计算亲和力高于SAP构象；而对于Λ-[EuL^2b]³⁺，则是SAP构象亲和力更高。

3. 配合物与DS2的分子相互作用揭示SAP/TSAP的不同结合模式

对于Δ-[EuL^1b]³⁺-TSAP，研究人员观察到独特的相互作用模式：Arg-410的胍基侧链与发色团的吡嗪环之间存在阳离子-π相互作用，Lys-414与发色团的苯基取代基也存在阳离子-π相互作用，此外还有疏水残基（Phe-403、Tyr-411、Leu-430）的稳定作用。相比之下，Δ-[EuL^1b]³⁺-SAP和Λ-[EuL^2b]³⁺-SAP/TSAP的结合模式高度相似，均缺乏Arg-410与吡嗪之间的阳离子-π相互作用，而是由Arg-410与手性悬臂的苯环形成阳离子-π相互作用。ProLIF相互作用指纹图谱的Tanimoto相似性分析定量支持了这一结论：Λ-[EuL^2b]³⁺的SAP和TSAP构象指纹相似度高达0.71，而Δ-[EuL^1b]³⁺的SAP和TSAP构象指纹相似度仅为0.33。

4. PCA相关性揭示不同结合对HSA运动的影响

通过对无配体HSA的四条1 μs轨迹进行PCA，前三个主成分（PC1-3）分别解释了26.1%、16.3%和10.4%的方差。将配体结合态模拟投影到同一PC空间后发现，Δ-[EuL^1b]³⁺结合态沿PC2方向呈现明显的SAP/TSAP分离，而Λ-[EuL^2b]³⁺结合态的分离则分布在PC1至PC3上。这表明不同立体化学构型的铕配合物结合会改变描述HSA骨架构象差异的主导运动模式。

5. HSA结合稳定了Eu配合物的构象

RMSD分析表明，Δ-[EuL^1b]³⁺-TSAP和Λ-[EuL^2b]³⁺-SAP在HSA结合态下的RMSD分布峰更尖锐，表明其构象稳定性较未结合态增加；而Δ-[EuL^1b]³⁺-SAP在结合态下稳定性反而降低，Λ-[EuL^2b]³⁺-TSAP则在结合前后无明显差异。这解释了实验中Δ异构体结合HSA后发生SAP→TSAP构象变化的观察结果。

四、结论

综合来看，这些特定的蛋白质-配合物相互作用的差异，不仅解释了实验测得的结合自由能趋势，也解释了加入HSA后Δ和Λ异构体CPL光谱的差异。从计算角度，相互作用分析、PCA和RMSD比较结果表明，配体结合对蛋白质构象行为和配体柔性的影响取决于探针的手性和构象。在蛋白质层面，与HSA药物位点2中Arg-410和Lys-414的选择性接触及阳离子-π相互作用，以及不同异构体配合物的结合，与HSA不同的构象偏好相关联。在配体层面，不同异构体在游离态和HSA结合态之间的柔性变化不同，这与实验中观察到的CPL光谱变化一致。总体而言，这些发现支持一个模型：探针构象影响蛋白质响应和结合行为，而计算结果允许在不同构象状态之间进行比较，而非直接模拟SAP到TSAP的构象转换。该论文发表于《Physical Chemistry Chemical Physics》。