基于复合密度泛函理论(r2SCAN-3c)洞察理性设计B、Ge、Mg和Si修饰Al12N12纳米笼：面向高效美金刚(Memantine)纳米载体设计用于阿尔茨海默病治疗

《South African Journal of Chemical Engineering》：Rational design of B, Ge, Mg, and Si-decorated Al??N?? nanocages from composite density functional theory insight: toward efficient memantine nanocarrier design for Alzheimer’s therapy

　　

摘要：美金刚(Memantine, MMT)被广泛用于延缓阿尔茨海默病进展，但其疗效受限于药物在到达中枢神经系统前因胃肠道吸收屏障及血脑屏障(blood–brain barrier, BBB)而部分损失。因此，开发高效的纳米材料基药物递送系统(drug delivery system, DDS)至关重要，Al12N12纳米笼是一种有前景的平台。本研究采用r2SCAN-3c水平的密度泛函理论(density functional theory, DFT)计算，在气相和水相环境中考察了MMT与原始及B、Ge、Mg、Si修饰Al12N12纳米笼的相互作用，以更接近生理条件。通过结构和热力学分析、电子结构表征（前线分子轨道(frontier molecular orbital, FMO)、态密度(density of states, DOS)、紫外–可见吸收(UV–Vis)）及拓扑相互作用分析（分子中的原子量子理论(quantum theory of atoms in molecules, QTAIM)、非共价相互作用–约化密度梯度(non-covalent interaction–reduced density gradient, NCI-RDG)、相互作用区域指示剂(interaction region indicator, IRI)、基于Hirshfeld分区的独立梯度模型(independent gradient model based on Hirshfeld partition, IGMH)、静电势(electrostatic potential, ESP)）评估吸附行为。结果表明，掺杂未显著改变纳米笼几何结构，但减小了带隙并改性了电子性质。关键发现为相依赖的掺杂剂切换效应：Al12N12Ge–MMT在气相中吸附最强（Eads= ?366.23 kJ/mol），而Al12N12Mg–MMT在水相中更优（Eads= ?196.97 kJ/mol），反映了溶剂化诱导的电子重组。美金刚吸附通过部分共价和范德华相互作用以放热方式自发进行，具有有利的水合热力学和适宜的解吸恢复时间。电子分析表明存在电荷重分布及电荷转移驱动的成键特征。上述发现确定原始及修饰Al12N12纳米笼是美金刚药物递送的有前景的理论候选体系，值得进一步实验研究。

论文解读：B、Ge、Mg、Si修饰Al₁₂N₁₂纳米笼负载美金刚的DFT研究

研究背景与意义

阿尔茨海默病（Alzheimer's disease, AD）是一种慢性、进行性、不可逆的神经退行性疾病，中重度AD常用治疗药物为美金刚（Memantine, MMT），其通过拮抗N-甲基-D-天冬氨酸（N-methyl-D-aspartate, NMDA）受体发挥神经保护作用。目前市售MMT为口服制剂，药物在抵达中枢神经系统（central nervous system, CNS）前会因胃肠道吸收障碍及血脑屏障（blood–brain barrier, BBB）的阻挡而大量损失，生物利用度低。纳米材料基药物递送系统（drug delivery system, DDS）因其可跨越BBB、生物相容性好、毒性低及靶向性增强等特点成为潜在解决方案。Al₁₂N₁₂纳米笼为宽禁带半导体，体外实验显示AlN纳米颗粒对正常细胞无毒且引发氧化应激小，具备作为安全DDS的材料基础；通过掺杂B、Ge、Mg、Si等元素可调控其电子性质（如缩小带隙）从而增强吸附性能。已有研究分别报道了B₁₂N₁₂及Al₁₂N₁₂纳米笼负载其他药物的潜力，但尚未有研究在统一计算协议下系统评估B、Ge、Mg、Si修饰Al₁₂N₁₂纳米笼在气相和水相对MMT的负载行为。本研究由Untung Gunawan等人发表于《South African Journal of Chemical Engineering》，采用高精度复合泛函r2SCAN-3c水平的密度泛函理论（density functional theory, DFT）计算，系统筛选四种掺杂剂修饰的Al₁₂N₁₂纳米笼作为MMT纳米载体，阐明掺杂剂身份如何调控电子性质和药物结合强度，并揭示相依赖的掺杂剂优选转换现象，为合理设计神经治疗用纳米载体提供量化与机理依据。

主要关键技术方法

研究人员采用Avogadro构建MMT分子及原始Al₁₂N₁₂纳米笼三维结构，通过GFN-xTB方法进行分子对接预筛选最稳定吸附构型，以ORCA 6.0.1程序在r2SCAN-3c复合泛函水平（含几何抵消校正gCP及D4色散校正，mTZVPP基组）对原始及B、Ge、Mg、Si表面修饰（adatom decoration）的Al₁₂N₁₂纳米笼及MMT吸附复合物分别在气相和SMD（solvent model density）水相模型中进行几何优化与频率计算（T=310.15 K，人体体温）。吸附能（E_ads）、吉布斯自由能变（ΔG）、焓变（ΔH）、熵变（ΔS）按标准公式计算；恢复时间（τ）由过渡态理论估算（尝试频率ν=10¹²s^?1）；前线分子轨道（frontier molecular orbital, FMO）分析计算电离能（IP）、电子亲和能（EA）、化学势（μ）、化学硬度（η）与化学软度（S）；采用Multiwfn 3.8进行态密度（density of states, DOS）、分子中的原子量子理论（quantum theory of atoms in molecules, QTAIM）、非共价相互作用–约化密度梯度（non-covalent interaction–reduced density gradient, NCI-RDG）、相互作用区域指示器（interaction region indicator, IRI）、基于Hirshfeld分区的独立梯度模型（independent gradient model based on Hirshfeld partition, IGMH）及静电势（electrostatic potential, ESP）分析；TD-DFT计算UV–Vis光谱；所有结构经频率验证为零虚频的真极小点。

研究结果

3.1. Optimization structure（优化结构与掺杂选择依据）

研究人员选用r2SCAN-3c复合泛函因其对非共价作用及几何参数精度高且含BSSE校正。选取B（ⅢA族）、Ge与Si（ⅣA族）、Mg（ⅡA族）作掺杂剂兼顾化学多样性、前期计算先例及生物相容性文献支持。采用表面修饰（adatom decoration）而非取代掺杂以保持笼化学计量与结构完整性。

3.2. Complex formation（复合体形成与几何结构）

xTB预筛选确定MMT以其仲胺N原子靠近纳米笼Al原子形成最稳定吸附取向（Al←N配位），因Al原子位于LUMO分布区具亲电性，MMT的N具孤对电子。优化后原始与修饰纳米笼Al–N键长（b₆₄≈1.80–1.85 ?，b₆₆≈1.78–1.79 ?）及六元环角（~126°）与文献吻合；掺杂引起局部晶格畸变及六元环角变化（B、Mg增大，Ge、Si减小），掺杂剂—最近N距离顺序为B< />< />< />_(MMT)键长在气相为~2.007–2.011 ?（Mg最短），水相为~1.965–2.007 ?（Ge最短），吸附引起六元环角及掺杂剂—N距微调但笼骨架基本保持不变。

3.3. Adsorption energy and thermodynamic parameters（吸附能与热力学参数）

所有体系E_ads为负值且绝对值>77.19 kJ/mol，属化学吸附范畴，过程自发放热。气相吸附稳定性顺序：Al₁₂N₁₂Ge–MMT（E_ads=?366.23 kJ/mol）> Al₁₂N₁₂Si–MMT > Al₁₂N₁₂–MMT > Al₁₂N₁₂Mg–MMT > Al₁₂N₁₂B–MMT；水相顺序：Al₁₂N₁₂Mg–MMT（E_ads=?196.97 kJ/mol）> Al₁₂N₁₂Ge–MMT > Al₁₂N₁₂Si–MMT > Al₁₂N₁₂–MMT > Al₁₂N₁₂B–MMT。除Mg体系外水相E_ads绝对值降低。ΔG、ΔH均为负（自发、放热），ΔS为负表明吸附为焔驱动且体系有序度增加（MMT自由度受限）。

3.4. Frontier molecular orbital (FMO)（前线分子轨道分析）

原始Al₁₂N₁₂带隙最大（气相3.12 eV，水相3.48 eV），掺杂显著缩窄带隙——气相Ge最小（E_g=0.32 eV），其次B(1.02 eV)、Si(1.16 eV)、Mg(1.50 eV)；水相顺序变为Si(1.15 eV)_g及最高电子亲和能（EA=3.33 eV）与较低电离能（IP=3.66 eV），化学软度最大（S=3.10 eV^?1），利于电荷转移；水相最优Mg体系虽带隙不是最小，但水合稳定了其电子态。吸附后各体系E_g进一步略降，HOMO向掺杂剂定域、LUMO于Al/掺杂剂离域，MMT的N→纳米笼发生电荷转移。

3.5. Density of states (DOS)（态密度分析）

TDOS显示掺杂及吸附均改变价带–导带间带隙宽度印证FMO结果。重叠布居态密度（overlapping population density of states, OPDOS）在价带区微正（>0），表明弱共价相互作用。吸附体系中TDOS主要由纳米笼片段贡献，说明纳米笼兼具亲核与亲电特性，FMO主要分布于纳米笼。

3.6. UV-Vis analysis（紫外–可见吸收分析）

TD-DFT显示MMT吸附引起吸收带红移（E_g减小），红移幅度与吸附强度相关——气相Ge体系λ_max从598 nm移至791 nm（Δλ=193 nm）最为显著，印证其强电子软化与电荷重分布；水相Mg体系从272 nm移至464 nm（Δλ=192 nm）红移明显且振子强度维持较高，反映水相下Mg–MMT复合体电荷转移增强。溶剂使部分体系发生蓝移，体现溶剂化对电子跃迁的调制。

3.7. Quantum theory of atom in molecules (QTAIM)（QTAIM拓扑分析）

Al–N_(MMT)键临界点上?²ρ>0且H_BCP<0，属部分共价相互作用（混合范德华+弱共价特征），|V/G|≈1.05–1.08为中等强度非真共价键。次要氢键BCP呈典型弱范德华（|V/G|<1）。各体系相互作用类型一致，吸附强弱源自掺杂调控的电子结构而非键型改变。

3.8. Noncovalent interaction–reduced density gradient (NCI-RDG)（NCI-RDG分析）

三维等值面显示纳米笼与MMT间主要为绿色区域（范德华作用主导），蓝色点状区（Al–N界面处弱静电/氢键），无强红色排斥区，验证了QTAIM判定之混合作用机制及热力学有利性，气/水相NCI模式相同。

3.9. Interaction region indicator (IRI)（IRI分析）

IRI等值面同样显示绿区（范德华主导）伴Al–N界面微小蓝区（部分共价特征），无红区，进一步确认混合作用本质及掺杂不改变基本作用类型。

3.10. Independent gradient model on Hirshfeld partition (IGMH)（IGMH分析）

IGMH散点图与等值面呈现内部蓝渐变至外绿之分层作用图像——内层部分共价特征、外层范德华主导，与QTAIM/NCI/IRI结论自洽，明确MMT–纳米笼吸附为"部分共价+范德华（主导）+静电/氢键（次要）"的多重协同作用。

3.11. Charge transfer and electrostatic potential (ESP)（电荷转移与ESP分析）

Mulliken分析显示孤立MMT的N原子电荷?0.498 a.u.，吸附后积累至?1.08~?1.13 a.u.（ΔQ≈?0.58~?0.63 a.u.），证实N→纳米笼显著电荷转移；Ge掺杂剂正电荷最大（+0.436气相，+0.584水相）对应最强气相吸附，Mg在水相正电荷减小但仍获最优水相吸附能，体现相依赖电荷调制。ESP图中MMT呈正电势红区、纳米笼呈正负混合蓝红区，吸附时间匹配——MMT缺电子区朝向纳米笼富电子区，佐证电子互补驱动吸附。

3.12. Solvent effect（溶剂效应）

SMD水相所有ΔE_sol<0（自发溶解有利）。Ge体系溶剂化能最小（?50.66 kcal/mol）但气相吸附最强；Mg体系溶剂化能最大（?77.05 kcal/mol）且水相吸附最优，说明水相性能取决于溶剂化稳定效应，产生掺杂剂优选反转（气相Ge→水相Mg）。

3.13. Recovery time（恢复时间/解吸时间）

基于ΔG估算τ，气相Ge体系τ≈1.26×10³⁰s（极长滞留），水相Mg体系τ≈8.28×10²s（数百秒适中），其余体系较短（秒或亚秒级）。表明修饰Al₁₂N₁₂可在保护转运（适度结合）与靶位释放（可解吸）间取得平衡，Mg在水相具较适宜的载–释动力学窗口。

3.14. Comparison with the prior study（与既往研究对比）

相比文献中B₁₂N₁₂基纳米笼负载MMT最强E_ads=?216.87 kJ/mol，本研究Al₁₂N₁₂Ge–MMT气相E_ads=?366.23 kJ/mol约强1.7倍，水相Al₁₂N₁₂Mg–MMT（?196.97 kJ/mol）亦优于已报道B₁₂N₁₂体系，显示Al₁₂N₁₂平台具更优递送潜力。

讨论与结论总结

研究人员在讨论中指出：本工作首次在统一r2SCAN-3c协议下系统考察了B、Ge、Mg、Si修饰Al₁₂N₁₂纳米笼负载MMT的气相与水相行为，发现掺杂不破坏笼骨架但显著缩窄带隙并调控电子性质；MMT通过Al←N配位以部分共价+范德华（主导）混合机制吸附，过程自发放热且伴随显著N→笼电荷转移；关键现象为相依赖掺杂剂切换——气相Ge修饰吸附最强（E_ads=?366.23 kJ/mol），水相Mg修饰最优（E_ads=?196.97 kJ/mol），源于溶剂化对不同掺杂中心电子态的差异稳定，强调生理相关条件下须做水相评估。恢复时间分析显示合适的解吸窗口，支持其作为MMT纳米载体的理论可行性。研究局限包括：xTB预筛非全局极小、SMD隐式溶剂未捕捉BBB特性/第一水合层/离子配对、未做Boys–Bernardi BSSE校正（gCP已控且趋势不受影响）、未涉及合成可行性/体内毒性/BBB穿透验证。结论翻译如下：

在本计算研究中，研究人员采用r2SCAN-3c水平DFT计算系统分析了原始及B、Ge、Mg、Si修饰Al₁₂N₁₂纳米笼与美金刚在气相和水相中的相互作用特征及吸附行为。结果表明，掺杂不显著改变纳米笼结构框架但持续调控电子性质（所有修饰变体带隙减小）。美金刚吸附通过具部分共价和范德华性质的混合相互作用进行，作用强度适中且反映电荷转移驱动的结合。所有体系吸附过程自发且放热。Al₁₂N₁₂Ge–MMT表现出最强气相结合，而Al₁₂N₁₂Mg–MMT在水相中最有利——这种相依赖掺杂剂切换反映了溶剂化诱导的电子重组，即水环境优先稳定Mg中心电子态。该发现证明气相电子优化不保证溶液相性能，是合理药物递送设计的关键认识。所有体系对美金刚结合表现出有利热力学特征（自发溶剂化、吸附致电子结构改变及光谱位移）。DFT计算确定原始及修饰Al₁₂N₁₂纳米笼是美金刚递送有前景的理论候选体系，值得开展实验研究。但本研究纯属计算并仅提供理论筛选结果，不能推断临床适用性。须承认局限性：SMD隐式溶剂无法捕捉BBB属性（紧密连接、转运蛋白、各向异性膜）、首层水合动力学或生理介质中离子配对效应；未评估纳米笼合成可行性、生物相容性或毒性，未评价体内稳定性或生物分布。应通过QM/MM分子动力学与显式溶剂细化计算、专用BBB渗透膜模拟、纳米笼合成及生理流体稳定性表征、体外药物载–释动力学、细胞毒性及体内分布与毒理研究予以验证，方可将本研究发现转化为AD治疗的实用制剂。