蒙脱石作为唑来膦酸化合物递送载体的分子模拟

《Surfaces》：Molecular Modeling of Montmorillonite as a Delivery Carrier for Zoledronic Compounds

　　

唑来膦酸化合物（zoledronic compounds）被广泛用作抗骨吸收和抗肿瘤药物，以预防骨丢失和治疗多发性骨髓瘤；然而，其药代动力学限制促使研究人员寻找能够改善其控释和生物利用度的替代递送系统。在本研究中，研究人员通过结合经验力场（FF）、密度泛函理论（DFT）和分子动力学（MD）模拟，研究了唑来膦酸（ZOL）及其Ca²⁺和Zn²⁺盐在蒙脱石（montmorillonite）层间空间中的稳定性、结构行为和吸附性质。该研究的主要目标是评估蒙脱石作为潜在药物递送系统（DDS）的适用性。具体而言，研究人员分析了唑来膦酸[1-(2-羟基-2-膦酸基-2-膦酰基乙基)-1H-咪唑-3-鎓]（ZOL）及其Ca²⁺和Zn²⁺盐的晶体多晶型结构。通过两种方法（FF和DFT）计算得到的晶体结构与已知实验数据吻合良好。此外，ZOL嵌入蒙脱石受限层间空间在能量上是有利的。研究人员还评估了几种层间阳离子（Na⁺、Ca²⁺和Zn²⁺）。MD模拟表明，ZOL在蒙脱石层间空间内采用稳定的受限构型，咪唑基团表现出较小的扭转角。此外，在模拟的酸性介质中，ZOL的脱附在能量上是有利的。研究人员的计算预测，这种粘土矿物在唑来膦酸化合物的控释方面具有巨大潜力。

**论文解读：蒙脱石作为唑来膦酸化合物递送载体的分子模拟研究**

**一、研究背景与问题**

唑来膦酸（ZOL）是一种氨基双膦酸盐，广泛用于治疗骨质疏松症、多发性骨髓瘤及乳腺癌、前列腺癌等恶性肿瘤。其作用机制包括抑制骨吸收和抗肿瘤活性。然而，ZOL的临床性能受到不利药代动力学特征的限制：静脉给药后血浆浓度迅速下降，导致生物利用度低。因此，开发能够提高ZOL生物利用度并实现持续释放的新型递送系统成为迫切需求。ZOL的分子结构包含两个膦酸基团和一个咪唑环，具有多电离状态和对二价/三价金属阳离子的强配位能力，同时与无机基质（如羟基磷灰石）有强亲和力。这些特性提示，杂化材料可能成为控释和维持治疗浓度的理想平台。

层状硅酸盐矿物，如蒙脱石（montmorillonite），因其2:1层状结构、大比表面积、溶胀性和高阳离子交换容量，在药物负载和递送中具有吸引力。蒙脱石层间空间的过量负电荷由可交换阳离子（如Na⁺、Ca²⁺）补偿，这些阳离子可通过离子交换和静电相互作用被有机阳离子或两性离子分子取代。药物与蒙脱石的相互作用可提高其溶解度和生物利用度。已有研究表明，蒙脱石可作为多种活性分子的载体，包括抗糖尿病药、抗炎药、抗生素和神经保护剂。

ZOL的强金属结合能力和多阴离子特性使其特别适合嵌入无机基质以实现控释。研究人员提出将蒙脱石用作ZOL的递送载体，通过将ZOL嵌入其层间空间，该载体可保护药物在胃酸环境中，从而探索口服给药的可能性。然而，这些药物-粘土复合物的物理化学性质尚未完全了解。因此，分子模拟方法成为在原子尺度探索这些杂化系统性质的有价值工具。

**二、研究内容与结论**

本研究采用经验力场（FF）、密度泛函理论（DFT）和分子动力学（MD）模拟，系统研究了ZOL及其Ca²⁺和Zn²⁺盐的晶体结构，以及它们在蒙脱石模型（SMEC）中的插层和脱附行为。研究人员构建了含Na⁺、Ca²⁺或Zn²⁺层间阳离子的蒙脱石超胞，并模拟了ZOL以中性、两性离子和阳离子形式嵌入层间空间的过程，以及ZOL与层间阳离子形成配位复合物的插层。计算结果表明，ZOL的插层在能量上是有利的，且层间阳离子类型影响插层结构和吸附能。MD模拟显示ZOL在层间空间中保持稳定构型，扩散系数较低。脱附模拟表明，ZOL在纯水介质中脱附能量不利，但在模拟的酸性NaCl介质中脱附能量有利，提示酸性条件可促进释放。

该研究的重要意义在于：从理论层面验证了蒙脱石作为ZOL控释载体的可行性，为实验制备ZOL-蒙脱石复合材料提供了原子尺度的结构预测和能量依据，尤其是XRD衍射图谱中d(001)间距的变化可作为实验监测的基准。论文发表在《Surfaces》。

**三、主要技术方法**

研究人员采用三种主要计算方法：（1）基于经验力场的INTERFACE力场（FF）进行晶体结构优化和插层模拟，使用Forcite代码；（2）基于密度泛函理论（DFT）的CASTEP和DMol3代码，采用GGA-PBE交换关联泛函和色散校正，用于精确计算电子结构和能量；（3）分子动力学（MD）模拟采用NVT系综、1 ns时长、Nosé恒温器，用于探索构象空间和评估扩散行为。此外，采用蒙特卡洛模拟退火方法生成有机分子插层构型。所有计算均在周期性边界条件下进行，模型包括NaSMEC、CaSMEC和ZnSMEC三种蒙脱石模型，以及从晶体数据（CCDC）中提取的ZOL多晶型及其Ca²⁺和Zn²⁺盐结构。

**四、研究结果**

**4.1 唑来膦酸晶体结构的优化**
研究人员通过比较INTERFACE FF和DFT计算与实验X射线衍射数据，验证了FF对双膦酸盐的适用性。两种方法均能较好重现ZOL三斜和单斜多晶型以及钙盐、锌盐的晶胞参数，偏差小于12%。DFT计算表明三斜多晶型比单斜多晶型更稳定（ΔE = 0.07 eV），且ZOL以两性离子形式存在，分子通过强氢键形成层状结构。在Ca²⁺和Zn²⁺盐中，阳离子呈八面体配位，与两个ZOL配体和两个水分子结合。

**4.2 ZOL嵌入蒙脱石层间空间**
研究人员将ZOL（中性、两性离子或阳离子形式）置于CaSMEC、NaSMEC或ZnSMEC层间中心进行优化。插层后d(001)间距显著增加（例如CaSMEC从14.07 ?增至14.48–16.0 ?），且层间水分子与矿物表面及ZOL形成氢键网络。ZOL的咪唑环保持近似垂直于甲基双膦酸基团的构型，但层间限制导致轻微扭转。吸附能计算显示，所有插层过程在能量上均有利：两性离子ZOL嵌入CaSMEC的吸附能为-4.96 eV（DFT）/-4.50 eV（FF）；阳离子交换（NaSMEC中ZOL⁺替换Na⁺）的吸附能为-7.57 eV；ZOL与Ca²⁺或Zn²⁺形成配位复合物插层时吸附能约为-1.2 eV，略低但仍有负值，表明插层稳定。

**4.3 ZOL的脱附**
研究人员模拟了ZOL从CaSMEC脱附至纯水介质和模拟NaCl水介质的过程。在纯水中脱附能Ed=+0.916 eV（能量不利）；在NaCl水介质中通过阳离子交换机制（ZOL⁺释出，Na⁺进入层间）的脱附能Ed=-0.726 eV（能量有利），表明酸性或含盐环境有利于ZOL释放。

**4.4 分子动力学**
MD模拟（1 ns）显示，ZOL在CaSMEC层间空间内达到平衡后保持稳定构型，咪唑环与矿物表面(001)垂直并通过N-H与表面氧形成氢键。ZOL分子的扩散系数较低（各向同性扩散系数4.26×10^-4 ?² ps^-1），水分子的扩散系数较高（2.36×10^-3 ?² ps^-1），且主要在xy平面内。ZOL在层间中心位置附近分布，C、N、P原子的浓度分布显示双峰，反映咪唑环的扭转和P-C-P平面的取向变化。

**五、讨论与结论**

讨论部分指出，INTERFACE FF经验证可有效模拟双膦酸盐及其与粘土矿物的相互作用，且与DFT结果一致，但DFT计算更为精确但计算成本高。虽然双膦酸基团存在多种构象（P=O基团syn/anti），但能量差异较小，不影响与粘土矿物的相互作用。分子模拟预测蒙脱石能通过能量有利过程插层ZOL，且酸性条件下有利于释放。这些计算仅描述初始和最终状态，未考虑温度效应或复杂生理环境，但可为口服给药实验提供指导。计算的d(001)间距增加可作为XRD实验的基准。

研究结论部分：ZOL嵌入蒙脱石导致d(001)间距显著增加，无论ZOL是中性、两性离子还是阳离子形式，插层过程在能量上均有利，阳离子交换机制同样有利。该理论预测ZOL插层到蒙脱石中是容易实现的，并能有效将ZOL稳定在无定形态，使其比晶态更易溶、生物利用度更高。这些预测对未来的实验工作非常有用，特别是通过XRD监测插层实验以及设计用于制药应用的药物-粘土复合材料。此外，虽然ZOL在纯水介质中脱附能量不利，但在低pH水介质中释放高度可能。这些预测表明这些材料可作为ZOL向环境控释的优秀候选者。最后，MD模拟证明分子动力学是搜索和采样最小能量结构的稳健方法，有效避免了静态优化中常见的局部极小值，并定性显示插层ZOL的扩散系数低，反映了药物在蒙脱石层间空间的高稳定性。总体而言，这项理论工作为未来专注于制备ZOL-蒙脱石递送系统的实验研究奠定了基础，该递送系统在治疗应用中具有巨大前景。