基于FRET的氮掺杂碳量子点-金纳米复合材料用于脱碱基DNA检测的分子动力学研究

本文通过分子动力学（MD）模拟，深入探讨了氮掺杂碳量子点（N-CQDs）与半胱胺修饰的金纳米颗粒（Cyst/AuNPs）组成的荧光共振能量转移（FRET）纳米复合材料在检测脱碱基DNA（Abasic DNA）中的应用。研究揭示了正常DNA与脱碱基DNA在纳米颗粒表面吸附行为的显著差异：完整DNA通过静电作用强烈吸附于Cyst/AuNPs，抑制FRET；而脱碱基DNA被AP核酸内切酶1（APE1）切割后亲和力降低，使得N-CQDs接近AuNPs激活FRET导致荧光猝灭。该工作为开发高灵敏度、无标记的DNA损伤生物传感器提供了分子水平的机理见解和设计原则。

分子动力学模拟方法

碳量子点结构采用基于六边形石墨烯片层堆叠的准球形模型，通过VMD软件的CD Builder工具构建，直径约2.91纳米。模型包含九层结构，表面功能化修饰有氨基（NH₂）、羧基（COOH）、羟基（OH）、羰基（C=O）及石墨氮等基团，功能化密度设为0.8。力场参数通过高斯16软件在B3LYP/6-31++G(d,p)水平进行几何优化，并采用Merz-Kollman方案计算部分电荷，非键参数基于广义AMBER力场（GAFF）设定。

DNA模型通过NAFlex网络服务器生成，脱碱基DNA序列为5′-GCTAC(X)GATCG-3′（X代表脱碱基位点），参照晶体结构PDB ID: 4IEM构建，模拟APE1酶切产物。正常DNA采用更长序列5′-GCTACCGATCGGCTACCGATCG-3′作为对照。标准核苷酸参数来自AMBER parmBSC1库，脱碱基位点参数通过ANTECHAMBER生成。

金纳米颗粒模型为直径2纳米的截角面心立方（FCC）晶格，包含309个金原子，暴露Au(111)和Au(100)晶面。表面修饰42个半胱胺（Cyst）分子，覆盖密度为5.7×10^–10mol/cm²。配体几何结构同样通过高斯16优化，电荷分配采用HF/6-31++G(d,p)水平的Merz-Kollman ESP拟合方法。

模拟体系包含两个N-CQDs（荧光供体）和两个Cyst/AuNPs（FRET受体），并加入2-3个DNA分子以模拟竞争吸附场景。体系在16.5纳米立方盒子内显式溶剂化，采用TIP3P水模型，添加22个Mg²⁺离子中和电荷。能量最小化后，进行5纳秒NPT平衡，随后进行100纳秒生产模拟。温度耦合采用V-rescale方法，压力耦合采用Parrinello-Rahman方法，步长为2飞秒，长程静电作用通过粒子网格Ewald（PME）方法处理。

分析指标包括均方根波动（RMSF）、溶剂可及表面积（SASA）、接触表面积（CSA）以及结合自由能（MM-PBSA）。RMSF计算原子相对于时间平均坐标的波动，CSA通过孤立组分与复合物SASA差值评估界面接触程度，结合自由能则分解为分子力学能、溶剂化自由能和熵贡献。

氮掺杂碳量子点在水溶液中的结构与动力学

氮掺杂碳量子点（N-CQDs）在水性环境中展现出独特的结构与动态特性，对其在生物传感中的应用至关重要。模拟的N-CQDs模型基于实验合成的绿色N-CQDs构建，表面修饰有NH₂、COOH、OH、C=O及石墨氮等官能团，其氮、氧含量与实验值高度吻合（氮含量模拟64.7% vs 实验64.5%）。密度分布分析显示N-CQDs具有层状石墨烯片结构，半径约2.92纳米，层间距0.36纳米，与石墨(002)晶面及已有模拟和实验数据一致。

结构中心形成疏水核心，而外层官能团高度水化，形成空间定义的疏水间隙，影响溶剂化及分子相互作用。在生理pH下，表面氨基和羧基分别以质子化（NH₃⁺）和去质子化（COO^–）形式存在，产生电荷异质性表面，增强溶解性并促进与生物分子的静电相互作用。

RMSF分析揭示从核心到表面的灵活性梯度：内部石墨层保持刚性，而表面原子波动显著（红色区域），表明表面层可作为分子识别的自适应位点。SASA图谱进一步证实氨基和羧基等官能团具有最高溶剂可及性，位于表面且可直接与水分子及分析物相互作用。这种结构特性使N-CQDs兼具柔性功能化外壳（利于生物分子识别）和稳定内核（维持结构完整性），特别适用于FRET生物传感器，其中一致的荧光发射和可控的表面相互作用对检测DNA损伤标志物至关重要。

DNA长度对FRET激活的影响

DNA长度通过调节纳米颗粒吸附动力学显著影响FRET效率。完整双链DNA（dsDNA）因其延伸构象和高负电荷密度，与带正电的Cyst/AuNPs产生强静电作用，空间上分离N-CQDs与AuNPs，抑制FRET。相反，APE1在脱碱基位点切割产生的短DNA片段对Cyst/AuNPs亲和力降低，允许N-CQDs接近AuNPs表面，激活FRET诱导的荧光猝灭。

短DNA（序列5′-GCTAC(X)GATCG-3′）与纳米颗粒的相互作用动力学显示：Cyst/AuNPs在约18纳秒时通过长程静电作用与DNA骨架结合，随后N-CQDs在27纳秒开始竞争结合，在40纳秒形成稳定复合物，最小距离达0.16纳米。氢键分析表明，核碱基与N-CQDs表面官能团（如羧酸盐和氨基）形成特异性氢键，稳定复合物。

而对于长DNA（序列5′-GCTACCGATCGGCTACCGATCG-3′），Cyst/AuNPs在18纳秒结合DNA后，N-CQDs在31纳秒附着于DNA，但供体-受体最小距离保持在约1.1纳米，显著影响FRET效率。氢键数量随时间增加反映复合物稳定性增强。长DNA因其更大表面积和更高电荷密度，与Cyst/AuNPs产生更强静电和范德华作用，从而阻断N-CQDs接近AuNPs，抑制FRET。这种DNA长度依赖的吸附行为突出了纳米颗粒结合调控FRET信号输出的关键机制。

脱碱基DNA与纳米颗粒相互作用的核苷酸贡献

为阐明DNA各组分对纳米颗粒吸附的贡献，分析了脱碱基DNA与N-CQDs和Cyst/AuNPs的结合动力学。RMSF图谱显示，与纳米颗粒表面直接接触的DNA区域波动性降低，表明局部结合稳定，而远端区域保持动态。

在N-CQDs–脱碱基DNA复合物中，核碱基和磷酸骨架残基均参与氢键和范德华作用。尽管骨架原子更接近N-CQDs表面，但核碱基形成更多氢键。这表明互补结合机制：带负电的DNA骨架介导初始静电吸引，而核碱基通过与N-CQDs表面官能团（如羧酸盐和氨基）的氢键稳定复合物。

类似地，脱碱基DNA与Cyst/AuNPs的结合涉及与金核的直接接触及与表面半胱胺配体的氢键。近界面核碱基形成更稳定作用，而远端片段保持构象灵活性。骨架和碱基均贡献氢键网络，骨架结合距离略近。这种协同结合模型强调静电作用和氢键共同稳定DNA–AuNP复合物。

N-CQDs与Cyst/AuNPs界面分析显示，N-CQDs上带负电的羧酸盐基团（COO^–）与Cyst分子上质子化氨基（NH₃⁺）间的静电作用主导相互作用，羰基（C=O）亦有贡献。其中COO^–形成最多氢键，凸显其在纳米颗粒间粘附的关键角色。这些相互作用对将供体-受体距离缩短至Förster半径内、激活高效FRET至关重要。

正常与脱碱基DNA在纳米材料上的结合特性比较

正常DNA与脱碱基DNA的结构完整性差异显著影响纳米颗粒吸附行为。正常DNA保持延伸构象和增强的碱基配对稳定性，而脱碱基DNA片段更短且更灵活。

正常DNA与Cyst/AuNPs复合物中，核碱基和磷酸骨架均通过氢键和静电作用参与结合，但骨架贡献更显著，形成更多氢键且平均距离更短。这种强锚定效应大幅限制N-CQDs接近AuNPs，抑制FRET激活。相比之下，脱碱基DNA系统中，因DNA长度较短和纳米颗粒表面覆盖度较低，N-CQDs可竞争性结合。

正常DNA与N-CQDs作用主要依赖骨架静电作用，核碱基参与较少，平均距离更长、氢键数更少，结合较弱。正常DNA的广泛表面覆盖屏蔽Cyst/AuNPs，阻止供体-受体接近。这些发现与实验报道一致，即长dsDNA可在金纳米颗粒表面形成稳定涂层，调节表面可及性和聚集行为。模拟进一步表明，DNA长度和结构完整性直接调控供体（N-CQDs）与受体（Cyst/AuNPs）间的竞争相互作用，为靶向脱碱基DNA损伤的FRET生物传感器设计提供机理基础。

结构与能量相互作用对FRET效率的影响

结构参数如SASA、CSA和结合自由能对调控DNA–纳米颗粒相互作用、复合物稳定性及供体-受体空间邻近性至关重要，进而影响FRET效率。

SASA映射显示Cyst分子、DNA磷酸骨架及N-CQDs官能团周围溶剂暴露度高。定量上，DNA的SASA最高（约185纳米²），其次为Cyst/AuNPs（约78纳米²）和N-CQDs（约67纳米²），表明DNA主导表面水化和界面反应性。CSA分析表明脱碱基DNA与Cyst/AuNPs形成最广泛接触（约17纳米²），而与N-CQDs和N-CQDs–AuNPs接触较弱（分别约5和8纳米²），凸显AuNPs作为DNA锚定支架并促进FRET重排的作用。

结合自由能分解进一步阐明作用强度和驱动力。脱碱基DNA–Cyst/AuNP复合物中，结合主要由范德华和静电作用稳定，极性溶剂化轻微不利于复合，总体负结合自由能证实热力学可行性。脱碱基DNA–N-CQD系统趋势类似但作用较弱，更易受溶剂化惩罚，净自由能仍为负值，表明生理条件下形成中等稳定复合物。N-CQDs–Cyst/AuNP系统显示带相反电荷官能团（如N-CQDs上COO^–与Cyst上NH₃⁺）间强静电吸引，辅以范德华接触，仅受适度溶剂化效应抵消，自由能景观高度有利。

与正常DNA系统比较揭示关键差异：正常DNA的SASA略低（约172纳米²），但与AuNPs和N-CQDs的接触更广（CSA分别约11和6纳米²），导致更宽表面覆盖。这种广泛吸附稳定纳米颗粒表面，阻碍空间重组所需供体-受体邻近，从而抑制FRET。能量上，正常DNA与Cyst/AuNPs保持强静电和范德华作用，但与N-CQDs亲和力较弱。这些发现与此前研究一致，即长完整DNA可钝化纳米颗粒，减少颗粒间作用，无损伤时抑制FRET。相反，APE1切割的短脱碱基DNA占据更小表面积、流动性更高，允许更近供体-受体间距，促进FRET信号激活。

脱碱基DNA检测的传感机制

N-CQDs与Cyst/AuNPs间的FRET传感机制关键取决于供体-受体距离，该距离受脱碱基DNA与正常DNA存在调控。N-CQDs发射强绿色荧光，Cyst/AuNPs有效吸收绿光，形成理想FRET对。颗粒间距离减小后，从N-CQDs到Cyst/AuNPs的非辐射能量转移导致荧光猝灭，为检测DNA损伤提供灵敏的距离依赖信号。

FRET效率（η_FRET）由公式η_FRET= R₀⁶/ (R₀⁶+ R⁶)描述，其中R为供体-受体距离，R₀为Förster半径（效率50%时的距离）。尽管N-CQD/Cyst-AuNP对的R₀未实验测定，MD模拟揭示了相对距离变化及其对FRET激活的影响。

静电势（ESP）分析显示，脱碱基DNA被APE1切割后，短片段对Cyst/AuNPs亲和力降低，允许带负电N-CQDs竞争结合纳米颗粒表面。所得纳米复合材料呈高度紧凑构型，供体-受体距离缩短（<1纳米），最小距离R_A–D达0.16纳米，使颗粒处于FRET激活区，导致高效荧光猝灭。相反，正常DNA作为刚性间隔物，维持N-CQDs与Cyst/AuNPs间更大分离，空间分布反映较大供体-受体距离（多>1纳米），最小R_A–D为1.1纳米，显著降低FRET效率，保持N-CQDs荧光。

这些观察揭示基于DNA结构完整性的独特传感机制：脱碱基DNA允许供体-受体邻近并激活高效FRET，而正常DNA维持分离阻止能量转移。N-CQD/Cyst-AuNP纳米复合材料的FRET响应由DNA完整性而非对脱碱基dsDNA的绝对选择性主导。脱碱基损伤引入局部双链变形和增加灵活性，减少纳米颗粒表面覆盖并促进供体-受体接近；而完整dsDNA保持刚性骨架，抑制纳米颗粒重排和FRET激活。正常DNA系统中颗粒间距离的更宽分布进一步抑制FRET，增强选择性。完整DNA骨架与Cyst/AuNPs间的强静电吸引稳定颗粒分离，而切割后脱碱基DNA缺乏连续骨架结构促进重排和FRET。这些发现表明N-CQD/Cyst-AuNP纳米复合材料可基于距离依赖FRET猝灭灵敏区分正常与脱碱基DNA。需强调，此FRET激活机制是设计特异性的，依赖于故意配对带正电或电荷异质N-CQDs与带负电金纳米颗粒对应物。替代CQDs表面电荷配置将改变DNA、N-CQDs和AuNPs间竞争结合层次，可能抑制供体-受体邻近和FRET效率。因此，本研究提出的传感机制应视为合理设计框架而非适用于所有CQD表面化学的普适行为。

结论

本研究通过原子尺度模拟揭示了FRET启用传感平台用于选择性脱碱基DNA检测的分子机制，利用氮掺杂碳量子点（N-CQDs）和半胱胺功能化金纳米颗粒（Cyst/AuNPs）。研究阐明DNA完整性差异（脱碱基vs正常）如何调控纳米颗粒吸附、颗粒间邻近性及最终FRET激活。脱碱基位点被APE1切割降低DNA对Cyst/AuNPs亲和力，使N-CQDs接近并诱导荧光猝灭；相反，完整DNA强吸附于纳米颗粒表面，维持空间分离并保留荧光。结构与能量分析强调静电作用和溶剂可及性在调节复合物稳定性及传感性能中的关键作用。脱碱基DNA片段促进更紧供体-受体间距（<1纳米），产生高效FRET，而正常DNA通过广泛表面覆盖抑制此作用。这些发现为FRET信号调制提供机理见解，并为开发选择性DNA损伤生物传感器确立设计原则。从实践角度，模拟为实验实施提供明确指导：应采用带正电或电荷异质N-CQDs与带负电AuNPs组合，以最大化脱碱基切割后电驱纳米颗粒重排和FRET激活。所提出的FRET启用N-CQD/AuNP平台为监测基因组完整性及诊断疾病相关DNA损伤提供有前景途径。