《Biology》：The Effects of Bisphosphonates Used in Osteoporosis Treatment on Breast Cancer: Analysis with Integrative Bioinformatics Methods, DFT, ADMET and Molecular Docking Analysis
编辑推荐：
本研究评估了3种双膦酸盐衍生物——阿仑膦酸盐、利塞膦酸盐和唑来膦酸盐——的结构、电子、药代动力学及受体结合特性，以探讨其在骨质疏松症与乳腺癌中的治疗相关性。基于B3LYP/6-31G(d,p)水平的密度泛函理论（DFT）计算显示，利塞膦酸盐具有最高的动力学稳定
本研究评估了3种双膦酸盐衍生物——阿仑膦酸盐、利塞膦酸盐和唑来膦酸盐——的结构、电子、药代动力学及受体结合特性，以探讨其在骨质疏松症与乳腺癌中的治疗相关性。基于B3LYP/6-31G(d,p)水平的密度泛函理论（DFT）计算显示，利塞膦酸盐具有最高的动力学稳定性（ΔE = 6.7468 eV），而唑来膦酸盐表现出更高的化学反应性（ΔE = 2.9669 eV）以及最强的非线性光学响应（β = 1.20 × 10?30 esu）。ADMET（吸收、分布、代谢、排泄与毒性）分析表明，所有化合物均具有可接受的安全性特征，但其高极性与低亲脂性可能限制口服生物利用度。针对11个与乳腺癌及骨代谢相关靶点开展的分子对接分析显示，这些化合物均具有较有利的结合亲和力，其中以唑来膦酸盐和利塞膦酸盐表现尤为突出。唑来膦酸盐与ESR2、VEGFR/KDR、GGPS1和FPPS之间表现出较强相互作用，而利塞膦酸盐则对BRCA2和MMP9显示出显著亲和力。生物信息学分析鉴定出GGPS1、FDPS、TNFSF11、ESR1、MMP9和BRCA2在乳腺癌组织中存在显著失调；生存分析进一步表明，FDPS、MMP9和BRCA2高表达与较差预后相关。网络分析突出显示了与甲羟戊酸代谢、激素信号传导、血管生成、细胞外基质重塑以及RANK/RANKL/OPG轴相关的通路。这些结果支持将双膦酸盐，尤其是唑来膦酸盐，重新定位用于乳腺癌相关骨病治疗的潜力。
该文发表于《Biology》，围绕“骨质疏松治疗药物能否在乳腺癌相关骨病中发挥超出抗骨吸收之外的作用”这一问题展开。研究背景在于：乳腺癌是女性高发恶性肿瘤之一，疾病进展过程中常伴发骨重塑失衡、骨破坏及骨转移相关并发症，不仅降低生活质量，还显著增加骨折风险。与此同时，肿瘤进展与骨代谢之间存在复杂双向调控，尤其涉及破骨细胞活化、骨微环境重塑、激素受体信号、血管生成及细胞外基质降解等过程。现有双膦酸盐已广泛用于抑制骨吸收和改善骨健康，但其在乳腺癌中的潜在分子作用尚未被系统厘清，因此有必要从多靶点、多层次角度重新评估其在乳腺癌相关骨病中的治疗价值。
在这一背景下，研究人员选取阿仑膦酸盐、利塞膦酸盐和唑来膦酸盐3种代表性双膦酸盐，综合应用量子化学、药代毒理预测、分子对接及生物信息学方法，对其结构性质、电子反应性、药代安全特征、与乳腺癌及骨代谢关键蛋白的结合能力，以及相关基因在乳腺癌中的表达和预后意义进行了整合分析。研究结论表明，这3种双膦酸盐均可作用于与乳腺癌进展和骨代谢相关的重要分子网络，其中唑来膦酸盐总体表现最佳：其化学反应性更高、与多个关键靶蛋白结合更强，并在整合分析中显示出较高的再利用潜力。该研究的重要意义在于，从计算层面为双膦酸盐在乳腺癌相关骨病中的药物重定位提供了机制依据，也提示甲羟戊酸通路、RANK/RANKL/OPG轴、激素受体信号和基质重塑通路可能是相关干预重点。
研究人员采用的主要技术方法包括：首先从PubChem获取3种双膦酸盐结构，并使用Gaussian 09在B3LYP/6-31G(d,p)水平进行密度泛函理论（DFT）计算，分析HOMO–LUMO、分子静电势（MEP）和非线性光学（NLO）参数；其次借助ADMETlab 2.0评估化合物吸收、分布、代谢、排泄与毒性特征；再次使用AutoDock Vina对11个来源于RCSB PDB的靶蛋白进行分子对接；最后结合TCGA-BRCA和GTEx队列，经GEPIA完成差异表达分析，经Kaplan–Meier Plotter完成生存分析，并通过STRING与GeneMANIA进行蛋白互作和功能网络分析。
在研究结果部分，论文首先从化学结构差异入手指出，阿仑膦酸盐具有柔性的脂肪族侧链，利塞膦酸盐含杂芳香性吡啶环，唑来膦酸盐则含电子密度更高的咪唑环；这些结构差异决定了三者在电子分布、极性和配位行为上的不同，并为后续反应性和结合能力差异提供了解释基础。
在“3.1. HOMO–LUMO Analysis”中，研究人员通过前线分子轨道（FMO）理论评估化学稳定性与反应性。结果显示，利塞膦酸盐的HOMO–LUMO能隙ΔE最高，提示其动力学稳定性最强、化学反应性最低；唑来膦酸盐ΔE最低，说明其电子柔性更高、化学反应性最强；阿仑膦酸盐处于中间水平。电离势、电子亲和势、化学硬度、柔度和亲电指数等参数进一步支持这一结论。总体而言，该部分说明3种双膦酸盐存在明确的“结构—电子反应性”对应关系，其中唑来膦酸盐更适于与靶标发生电子层面的相互作用。
在“3.2. Molecular Electrostatic Potential (MEP) Analysis”中，研究人员利用分子静电势分析各化合物的亲电/亲核反应中心及潜在结合位点。结果表明，3种双膦酸盐最富电子区域主要位于双膦酸基团及羟基氧原子，显示出较强的氢键受体能力；缺电子区域则集中于各自含氮或相关氢原子部位。利塞膦酸盐表现出最宽的静电势差，提示其表面极化和电荷分离更明显，理论上更易在蛋白结合腔内建立强静电作用和稳定氢键。该部分结果为后续分子对接中结合强弱差异提供了静电学依据。
在“3.3. Non-Linear Optical Properties (NLO)”中，研究人员比较了3种化合物的偶极矩、平均极化率和一级超极化率。结果显示，唑来膦酸盐具有最高偶极矩和最强非线性光学响应，提示其分子极性和电场响应能力最强；利塞膦酸盐则表现出最高电子云可变形性。与尿素标准相比，3种双膦酸盐的非线性光学参数均明显增强。该部分主要说明不同双膦酸盐不仅在药理相关电子性质上存在差异，其分子电学特征也受取代结构显著影响。
在“3.4. ADMET Analysis”中，研究人员评估了3种化合物的理化性质与药代安全性。结果显示，它们分子量较低，具备一定药物可及性；但由于氢键供体/受体数较多、拓扑极性表面积（TPSA）较高且logP为负值，整体呈现高极性、低亲脂性特征。这意味着虽然水溶性较好，但被动跨膜扩散能力有限，可能导致口服生物利用度偏低。分布方面，其血脑屏障（BBB）通透性较低，符合骨靶向用途；代谢方面，对CYP2D6和CYP3A4抑制概率较低，提示药物相互作用风险有限；毒性方面，hERG阻断、AMES致突变性和肝毒性风险均较低。该部分结论是：3种双膦酸盐具有可接受的药代和安全特征，但高极性是限制口服利用度的关键因素，同时也可能促进其在骨组织中的选择性富集。
在“3.5. Molecular Docking Analysis”中，研究人员针对11个与乳腺癌进展和骨代谢相关的蛋白开展对接，结果整体显示3种化合物均具有有利的结合能力。就总体表现而言，唑来膦酸盐最佳，尤其对FPPS、ER-β/ESR2和COL18A1显示最强结合评分，对GGPS1也具有较高亲和力，这与含氮双膦酸盐抑制甲羟戊酸通路酶、抑制破骨细胞骨吸收的已知机制相一致。利塞膦酸盐在BRCA2、MMP9、VEGFR和ABCB1/MDR等靶点上表现突出，提示其可能影响DNA修复、基质降解、血管生成和多药耐药相关过程。阿仑膦酸盐虽然整体评分略低，但对VEGFR、MMP-9和FPPS仍保有较好结合能力。进一步的二维相互作用分析表明，阿仑膦酸盐的结合主要依赖氢键和静电作用；利塞膦酸盐由于芳香环存在，额外增加了π–π堆积和疏水作用，因此结合模式更丰富；唑来膦酸盐则呈现出多模式结合特征，在多个靶点上兼具氢键、静电、范德华力及π型相互作用，显示出更强的多靶点适应能力。该部分核心结论是，唑来膦酸盐和利塞膦酸盐可能比阿仑膦酸盐具有更优的乳腺癌相关分子干预潜力。
在“3.6. Bioinformatics Analysis”中，研究人员整合网络、表达和生存数据，揭示相关靶基因在乳腺癌中的生物学意义。STRING蛋白互作网络显示，所研究基因之间的相互作用数显著高于随机预期，提示这些基因构成了生物学上紧密相关的功能网络。ESR1、ESR2、KDR、MMP9、COL18A1、TNFSF11和TNFRSF11B位于网络中心，说明激素受体信号、血管生成、骨重塑和基质调控在网络组织中具有关键地位。GeneMANIA分析进一步表明，这些基因关系主要由物理相互作用和共表达构成，并延伸至VEGFC、THBS2、CXCL5、CXCL6、TNFSF13、BRCA2和HOXC8等关联基因。表达分析显示，GGPS1、FDPS、TNFSF11、ESR1、MMP9和BRCA2在乳腺癌组织中显著上调，而ESR2、KDR和ABCB1显著下调，COL18A1无显著差异。生存分析显示，FDPS、MMP9和BRCA2高表达与较差总生存相关；GGPS1、TNFSF11/RANKL、TNFRSF11B/OPG、ESR2、ESR1、KDR/VEGFR2和ABCB1较高表达则与较好生存相关。该部分结果支持如下判断：甲羟戊酸通路、激素受体信号、血管生成、细胞外基质重塑及RANK/RANKL/OPG轴共同参与乳腺癌进展，而FDPS、MMP9和BRCA2可视为潜在不良预后标志物。
讨论部分的核心在于整合不同层次证据，建立双膦酸盐可能作用于乳腺癌相关骨病的机制框架。研究人员认为，FDPS和GGPS1作为甲羟戊酸通路关键酶，在乳腺癌中表达异常，且该通路与肿瘤细胞增殖、迁移及抗凋亡有关，因此双膦酸盐通过抑制该通路具有合理的机制基础。MMP9高表达且与不良预后相关，提示其在基质降解、肿瘤侵袭和骨转移中具有重要作用。ESR1与ESR2在网络中的中心地位及其与较好生存的关联，提示激素受体相关机制仍是乳腺癌重要的生物学决定因素。RANK/RANKL/OPG轴的改变则进一步说明骨微环境调控与乳腺癌进展密切耦联。总体上，论文将甲羟戊酸代谢、激素信号、血管生成、基质重塑和骨代谢轴线整合为一个相互关联的病理网络，并据此支持双膦酸盐尤其是唑来膦酸盐的再利用价值。
研究结论可译述为：本研究通过DFT、ADMET、分子对接和生物信息学方法，对阿仑膦酸盐、利塞膦酸盐和唑来膦酸盐进行了整合性计算评估。结果显示，3种双膦酸盐具有不同的电子特性，其中唑来膦酸盐化学反应性最高，利塞膦酸盐动力学稳定性最强。ADMET分析提示其药代动力学和安全性总体可接受，但高极性可能限制口服生物利用度。分子对接显示，这些化合物可与多种参与乳腺癌进展和骨重塑的蛋白形成有利相互作用，其中唑来膦酸盐总体结合能力最强，尤其针对ESR2、VEGFR/KDR、GGPS1和FPPS；利塞膦酸盐则对BRCA2和MMP9具有显著相互作用。生物信息学结果进一步证实，GGPS1、FDPS、TNFSF11、ESR1、MMP9和BRCA2在乳腺癌组织中显著失调，且多种靶基因与患者生存显著相关。STRING与GeneMANIA分析强调了甲羟戊酸通路、激素受体信号、血管生成、细胞外基质重塑及RANK/RANKL/OPG轴在乳腺癌进展中的参与。总体而言，这些发现提示双膦酸盐，特别是唑来膦酸盐，除抑制骨吸收外，还可能通过调控与乳腺癌相关的分子通路发挥生物学作用，但仍需进一步实验和临床研究验证。
