[论文解读] Ab initio Investigation of Adsorption Characteristics of Bisphosphonates on Hydroxyapatite (001) Surface

Mun-Hyok Ri, Yong-Man Jang+2 more|arXiv (Cornell University)|Apr 26, 2016

Bone health and treatments参考文献 29被引用 3

一句话总结

本项从头算研究利用密度泛函理论，研究了七种双膦酸盐（BPs）在羟磷灰石（001）表面的吸附行为。结果表明，帕米膦酸盐的吸附能最高（-3.33 eV），其次为阿仑膦酸盐（-3.21 eV），这主要归因于氮和羟基基团产生的强氢键作用及电子离域效应。计算结果与实验测得的亲和力排序一致，证实吸附能是决定双膦酸盐与骨矿结合亲和力的关键因素。

ABSTRACT

The structures of some bisphosphonates (clodronate, etidronate, pamidronate, alendronate, risedronate, zoledronate, minodronate) were obtained and analyzed, and their adsorption energies onto hydroxyapatite (001) surface were compared to find out ranking order of binding affinity, which shows that the adsorption energy is the largest for pamidronate, followed by alendronate, zoledronate, clodronate, ibandronate, the lowest for minodronate and etidronate.

研究动机与目标

• 确定关键双膦酸盐（BPs）在羟磷灰石（HAP）表面的相对结合亲和力，HAP是骨骼中主要的矿物成分。
• 阐明控制BPs在HAP(001)表面吸附的原子尺度机制，包括电子和几何因素。
• 将计算得到的吸附能与实验测得的结合亲和力排序进行关联，以验证计算预测的可靠性。
• 分析电荷转移及轨道相互作用（HOMO-LUMO），以揭示BP-HAP复合物的稳定性来源。

提出的方法

• 采用广义梯度近似（GGA）结合PBE泛函，进行从头算密度泛函理论（DFT）计算。
• 构建HAP(001)表面的周期性 slab 模型，BPs吸附于钙终止表面。
• 吸附能按公式 E_ads = E_{BP+HAP} - (E_{BP} + E_{HAP}) 计算，能量收敛精度达 10^{-6} eV。
• 采用Hirshfeld电荷分析研究吸附过程中的电荷重分布，评估电子供给与表面稳定化作用。
• 通过分子轨道分析（HOMO与LUMO）识别羟基（-OH）和氨基（-NH2）等官能团的电子供体/受体角色。
• 对吸附构型进行几何优化与分析，研究BP中羟基与HAP表面氧原子之间的氢键作用。

实验结果

研究问题

• RQ1 双膦酸盐在HAP(001)表面的相对吸附能顺序如何？与实验排序相比有何异同？
• RQ2 官能团（如-OH、-NH2、杂环氮）如何影响BP-HAP复合物的电子与几何结构？
• RQ3 BP中羟基与HAP表面氧原子之间的氢键在吸附稳定化过程中起何作用？
• RQ4 通过Hirshfeld分析量化电荷转移后，如何影响BP-HAP界面的稳定性？
• RQ5 HOMO-LUMO相互作用在多大程度上可解释含氮BPs结合亲和力的增强？

主要发现

• 帕米膦酸盐的吸附能最高（-3.33 eV），其次为阿仑膦酸盐（-3.21 eV）、唑来膦酸盐（-2.92 eV）、氯膦酸盐（-2.35 eV）、伊班膦酸盐（-2.31 eV）、利塞膦酸盐（-2.25 eV）、米诺膦酸盐（-2.17 eV）和依替膦酸盐（-2.17 eV）。
• 计算得到的吸附能排序与核磁共振（NMR）、FPLC及竞争结合实验测得的实验亲和力排序高度吻合，验证了计算模型的可靠性。
• BP的羟基与HAP表面氧原子之间的氢键是吸附构型中关键的稳定作用力。
• Hirshfeld电荷分析显示，吸附过程引起电荷重分布，使HAP表面电荷分布更接近体相晶体的电荷特征，从而实现表面稳定化。
• 含氮BPs（如帕米膦酸盐、阿仑膦酸盐）中的氮原子在吸附过程中表现出电子供体特性，这一结论由HOMO-LUMO分析支持。
• P–C–P主链结构对R2基团变化敏感，其中含氮R2基团显著增强了结合亲和力。