抗癌药物Afinitor中光谱特征与配体-蛋白相互作用的定量关联：一项实验-计算集成研究

《Scientific Reports》：Quantitative correlation of spectroscopic signatures with ligand–protein interactions in anti-cancer drug Afinitor: an integrated experimental–computational study

　　

对抗癌药物的详细分子层面理解对于提高治疗效果和指导合理药物设计至关重要。依维莫司（Everolimus，商品名Afinitor）是哺乳动物雷帕霉素靶蛋白（mammalian target of rapamycin，mTOR）通路的重要临床抑制剂，广泛用于癌症治疗；然而，其光谱特征与配体-蛋白相互作用之间基于定量的关联仍未被充分探索。本研究采用集成实验-计算方法，结合傅里叶变换红外（Fourier transform infrared，FT-IR）光谱、紫外-可见（ultraviolet-visible，UV-Vis）光谱和分子对接模拟。光谱分析用于表征官能团和电子结构，而对接模拟用于研究与FK506结合蛋白（FK506 binding protein，FKBP12）及FKBP-雷帕霉素结合（FKBP–rapamycin binding，FRB）结构域的相互作用。FT-IR分析揭示了高密度含氧官能团，包括羟基和羰基部分，其振动频率表明强极化电子环境。分子对接显示与FKBP12（?9.7和?9.6 kcal·mol?1）及FRB结构域（?8.5和?6.6 kcal·mol?1）的有利结合亲和力。详细相互作用分析表明，这些官能团对应于特定相互作用原子（如O66、O67、O63和O36），与TYR82、THR85和GLU54等关键残基形成可量化氢键（1.7–2.9 ?）和静电相互作用（~4.37 ?）。278 nm处的UV-Vis吸收最大值对应于4.46 eV的最高占据分子轨道-最低未占分子轨道（Highest Occupied Molecular Orbital–Lowest Unoccupied Molecular Orbital，HOMO–LUMO）能隙，表明中等电子极化率，支持结合过程中的电荷重分布。该研究建立了定量且基于机制的构效-光谱-相互作用关系，证明光谱可观测量编码了调控配体-蛋白相互作用倾向的局部电子环境。结合亲和力被证明源于由官能团空间分布实现的多种非共价相互作用协同网络。该集成框架为药物-蛋白相互作用提供了预测性见解，并为理性设计和优化靶向mTOR的治疗药物提供了坚实基础。

该研究发表于《Scientific Reports》。癌症是全球重大健康挑战，发病率持续上升，亟需改进治疗策略并深入理解临床抗癌药物的分子机制。哺乳动物雷帕霉素靶蛋白（mTOR）通路抑制剂是重要的靶向抗癌药，依维莫司（Everolimus，Afinitor）作为其衍生物通过结合FK506结合蛋白12（FKBP12）抑制mTOR活性，阻断细胞周期进展以抑制增殖，广泛用于恶性肿瘤治疗及器官移植排斥预防。尽管依维莫司的药理活性和分子结构已明确，但整合光谱表征与FKBP相关靶点的分子对接相互作用分析的研究有限，光谱识别的官能团在蛋白结合中的具体作用关注不足，建立此类关联可深化药物-靶点识别的分子决定因素理解，因此开展本研究。研究人员通过集成傅里叶变换红外（FT-IR）光谱、紫外-可见（UV-Vis）光谱与分子对接模拟的实验-计算方法，表征依维莫司的官能团、电子结构及与FKBP12、FKBP-雷帕霉素结合（FRB）结构域的相互作用，建立光谱特征与配体-蛋白相互作用的定量关联，揭示其结合机制的结构基础。

作者开展研究用到的主要关键技术方法包括：采用FT-IR光谱（Perkin Elmer光谱仪，4000–450 cm^?1范围）和UV-Vis光谱（PerkinElmer Lambda 35，190–1100 nm范围，乙醇溶剂）进行光谱表征；从蛋白质数据库（Protein Data Bank，PDB）获取FKBP12（PDB ID：1FKL）和FRB结构域（PDB ID：1AUE）的高分辨率晶体结构，经BIOVIA Discovery Studio Visualizer预处理（去结晶水、加极性氢、分配Kollman电荷）后存为PDBQT格式；依维莫司结构取自PubChem并能量最小化优化后转PDBQT格式；使用AutoDock Tools 1.5.6（Lamarckian遗传算法）进行分子对接，定义20 ?×20 ?×20 ?网格盒覆盖活性位点，通过重对接验证（均方根偏差Root Mean Square Deviation，RMSD<2.0 ?）确认可靠性；用BIOVIA Discovery Studio Visualizer分析复合物氢键、静电、疏水及π–π堆积相互作用并生成二维/三维互作图。

研究结果如下：

FT-IR Spectrum analysis of Afinitor。通过FT-IR光谱分析依维莫司的振动频率，识别3449、3375、3354 cm^?1的OH伸缩振动，3350、2698、2653 cm^?1的CH伸缩，2902 cm^?1的CH₃不对称伸缩、3285 cm^?1的CH₃对称伸缩、2977 cm^?1的CH₃次甲基伸缩，2955、2915、2877、2782 cm^?1的CH亚甲基伸缩，1578、1654、1648 cm^?1的C=O伸缩，1361 cm^?1的νC-C，1560 cm^?1的βCH、1475 cm^?1的βHCH弯曲，1449 cm^?1的CH₃面内弯曲（CH₃ib、CH₃ir），1419 cm^?1的βCH₃ob和βCH₃or，1388、1282 cm^?1的CH₂摇摆，1341、1159 cm^?1的CH₂扭曲，1308 cm^?1的βCH₂（CH₂剪式），1288 cm^?1的βNCH₂，1252 cm^?1的βCOH，1266、1225 cm^?1的C-N伸缩，1210、1095、1068、1018、948 cm^?1及1044、988 cm^?1的νC-C，915、734、701、657 cm^?1的环三角变形、对称/不对称变形及扭转，892、833、769 cm^?1的ωCH（面外弯曲），877 cm^?1的βCH₂摇摆，808 cm^?1的环三角变形，597、536 cm^?1的ωCC，631 cm^?1的ωCCN，581、518 cm^?1的τCCOH，552、475、462 cm^?1的βOCO、βCCO及环弯曲扭转，表明药物含高密度含氧官能团且电子环境强极化。

Electron absorption spectra。通过UV-Vis光谱获得依维莫司在乙醇中的吸收光谱，最大吸收波长278 nm，依据E=hc/λ及E(eV)=1240/λ(nm)计算得最高占据分子轨道-最低未占分子轨道（HOMO–LUMO）能隙为4.46 eV，表明中等电子极化率，支持结合时的电荷重分布，分子轨道跃迁对应HOMO至LUMO的电子激发。

Molecular docking analysis。通过分子对接模拟依维莫司与FKBP12（PDB ID：1FKL）及FRB结构域（PDB ID：1AUE）的相互作用，最佳结合亲和力为FKBP12（?9.7和?9.6 kcal·mol^?1）、FRB（?8.5和?6.6 kcal·mol^?1），显示对FKBP12更强偏好；Ramachandran图验证蛋白结构质量（残基多在允许区）。FKBP12对接中，O66与ALA84:HN、THR85:HN、THR85:HG1形成3个常规氢键，O67与TYR82:HH、O63与THR85:HG1各形成1个常规氢键，O37、O39与GLU54:CA形成碳氢键，O36与GLU54:OE2形成静电作用，O40、O11、O67、O41形成π-供体氢键与HIS87相互作用；FRB结构域对接中，N4与GLU2060:OE2静电作用，H134、O39、H145、H150、O66、O41、O11、O61、H54、H152、H153分别与GLU2060:OE2、ARG2087:HE、GLU2053:OE2等形成常规氢键，结合稳定性由氢键、静电、疏水接触协同维持，符合依维莫司已知机制。

Functional group analysis of Afinitor。依维莫司含大环内酯环（macrocyclic lactone ring）关键结构，羟基（–OH）通过氢键增强与mTOR活性位点氨基酸残基结合、提高水溶性及参与葡糖醛酸化代谢，羰基（C=O）因极性参与氢键、具亲电性、维持分子稳定及可代谢氧化，两者均可通过化学修饰优化药效；FT-IR识别的含氧官能团（羟基、羰基）与对接中稳定配体-蛋白复合物的氧原子（O66、O67、O63、O36）直接参与氢键和静电作用相对应，建立光谱特征与相互作用位点的关联。

Quantitative correlation between spectroscopic signatures and ligand–protein interactions。振动频率作为相互作用潜力间接描述符：Δν∝ρ_e∝E_int（Δν为振动频率偏差，ρ_e为局部电子密度，E_int为相互作用强度），OH（3449–3354 cm^?1）和C=O（1654–1648 cm^?1）的红移表明极化增强氢键能力；原子水平对应光谱识别官能团的氧原子（O66、O67、O63、O36）与TYR82、THR85、GLU54形成1.7–2.9 ?氢键及~4.37 ?静电作用；结合自由能ΔG_bind=∑E_i^H-bond+∑E_j^{electrostatic}+∑E_k^vdw，源于多非共价相互作用协同网络而非单一作用，FKBP12的?9.7/?9.6 kcal·mol^?1亲和力符合多点稳定机制；UV-Vis的4.46 eV HOMO–LUMO能隙表明中等电子极化率促进电荷重分布，增强氢键与静电贡献；集成框架为光谱（FT-IR、UV-Vis）→电子结构与官能团极性→相互作用倾向→对接几何与能量，证明光谱可观测量可预测配体-蛋白相互作用行为。

讨论部分总结：研究整合光谱表征与分子对接，FT-IR识别的含氧官能团（羟基、羰基）直接参与对接中FKBP12关键残基（TYR82、THR85、GLU54）的氢键与静电作用，建立定量构效-光谱-相互作用关系，振动频率反映局部电子密度与相互作用潜力，HOMO–LUMO能隙关联电子极化与结合电荷重分布，结合亲和力源于多非共价协同网络；局限在于仅用静态分子对接未涵盖分子动力学与自由能计算，缺NMR与生物实验验证；未来需扩展分子动力学、结构-活性关系及功能基团修饰与生化验证以优化mTOR靶向药物。

结论部分翻译：本研究通过集成FT-IR光谱、UV-Vis光谱和分子对接模拟，对抗癌药物依维莫司（Afinitor）进行了全面的实验-计算研究，以在分子层面阐明其配体-蛋白相互作用机制。光谱分析揭示了高密度含氧官能团，尤其是羟基和羰基部分，其振动特征反映了强极化电子环境。这些光谱特征被证明直接与对接模拟中识别的原子特异性相互作用位点对应，其中氧原子（如O66、O67、O63和O36）与TYR82、THR85和GLU54等关键残基形成量化氢键（1.7–2.9 ?）和静电相互作用（~4.37 ?）。重要的是，本研究建立了定量构效-光谱-相互作用关系，证明振动频率作为局部电子密度和相互作用潜力的间接描述符，为理解光谱识别官能团如何调控配体-蛋白结合行为提供了机制基础。观察到的结合亲和力（FKBP12为?9.7和?9.6 kcal·mol^?1）源于多种非共价相互作用的协同网络而非单一主导相互作用，突出了多点稳定在分子识别中的作用。此外，UV-Vis衍生的HOMO–LUMO能隙（4.46 eV）表明中等电子极化率，促进结合过程中的电荷重分布，从而支持观察到的静电和氢键相互作用，这确立了电子结构作为连接光谱可观测量与相互作用能量的额外环节。总体而言，本研究超越传统对接分析，引入了基于机制且定量支持的框架，桥接实验光谱与计算相互作用分析。提出的构效-光谱-相互作用范式证明光谱测量可预测配体-蛋白结合行为，为理性设计和优化靶向mTOR的治疗药物提供坚实基础。