从气相复合物到DFT计算：探索茶碱对胞嘧啶和腺嘌呤的分子识别选择性

《Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy》：Probing the selectivity of theophylline toward cytosine and adenine: Insights from gas-phase complexes and DFT calculations

本文通过质谱分辨红外/紫外光谱与DFT计算，在气相中研究了茶碱与胞嘧啶/腺嘌呤二元复合物的结构、成键与稳定性。研究发现二者均形成由强NH···N/NH···O=C双氢键稳定的构象，其中模拟适配体结合模式的构象THEOP·C-1比对应的THEOP·A-2结合能高约20 kJ/mol，为茶碱适配体对胞嘧啶的偏好性提供了分子层面的基础性解释，揭示了“锁钥”模型中关键“栓-槽”（Cut-Pin）相互作用的本质，对理解核酸适配体高选择性分子识别机制具有重要意义。

　　

分子识别是化学与生物学中的核心概念，它调控着酶与底物、抗体与抗原乃至核酸碱基配对的多种生命过程。其基础是分子间的非共价相互作用（Non-Covalent Interactions, NCIs），如氢键和π-π堆积。在这一领域中，适配体（Aptamers）——能够以高亲和力、高特异性折叠并结合特定靶分子的短链RNA或DNA——因其在诊断、治疗和分子传感中的应用潜力而备受关注。茶碱（Theophylline, THEOP），一种结构与咖啡因相似的甲基黄嘌呤类药物，是研究最深入的小分子适配体靶标之一。其RNA适配体对茶碱的选择性甚至可达咖啡因的10,000倍，这种惊人的特异性源于适配体结合口袋中茶碱分子与核苷酸残基之间复杂的非共价相互作用。

尽管“锁钥”模型（Lock-and-Key）常被用来描述这种分子识别，但更精细的“栓-槽”（Cut–Pin）模型能更好地解释其本质。在这个模型中，“槽”（Cut）代表锁（适配体）中的空腔或沟槽区域，而“栓”（Pin）则是钥匙（配体）上与其互补匹配的突起或官能团。每一个“栓-槽”对都代表一组分子间的接触点。理解茶碱与核碱基之间相互作用的分子基础，对于理性设计适配体生物传感器和药物至关重要。然而，在溶液中，溶剂效应和复杂的分子环境会干扰对最基础、最本质相互作用的精确观测。为了剥离这些复杂因素，研究分子识别的内在驱动力，科学家们将目光投向了气相环境。

在最新发表于《Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy》的研究中，Maider Parra Santamaria、Camilla Calabrese、Sonia Melandri、José A. Fernández、Ander Camiruaga和Imanol Usabiaga团队合作，将茶碱与胞嘧啶（Cytosine, C）或腺嘌呤（Adenine, A）的复合物引入气相，利用高分辨光谱与理论计算相结合的手段，直接探究了茶碱对这两种核碱基的内在结合偏好性，为适配体的高选择性识别提供了关键见解。

研究人员主要采用了以下几项关键技术方法：首先，利用结合激光烧蚀的超声速膨胀技术，在气相中生成了茶碱-胞嘧啶（THEOP·C）和茶碱-腺嘌呤（THEOP·A）二元复合物，有效消除了溶剂和基质的干扰。其次，通过质量分辨的共振增强多光子电离（REMPI）光谱和离子损耗红外（IDIR）光谱，获得了复合物的电子光谱和构象选择性的红外光谱。最后，结合密度泛函理论（DFT）计算，在M06-2×/6–311++G(d,p)和B3LYP-GD3BJ/def2TZVP等多个理论水平上，对复合物的构象空间、结合能和振动光谱进行了系统模拟，并与实验光谱对比以指认结构。

3.1. 构象空间的探索

通过DFT计算对THEOP·C和THEOP·A复合物的势能面进行系统搜索，发现存在三类主要相互作用机制：I) 由两个强氢键稳定的复合物；II) 单氢键结构；III) π-π堆积相互作用复合物。研究聚焦于最稳定、选择性更高的I类（双氢键）结构。对于THEOP·C，考虑了胞嘧啶的酮式和烯醇式互变异构体。最稳定的构象（如THEOP·C-1和THEOP·C-2）均由强双氢键稳定。值得注意的是，THEOP·C-1的相互作用模式与已知的茶碱-适配体复合物中观察到的模式相似，即茶碱与胞嘧啶的沃森-克里克（Watson–Crick）面结合。对于THEOP·A，同样有三个由双氢键稳定的最稳定构象（THEOP·A-1, A-2, A-3）。其中，THEOP·A-2也表现出与适配体模型相似的相互作用。理论分析揭示了一个关键发现：在模拟适配体结合模式的构象中，THEOP·C-1的结合能比THEOP·A-2高出约20 kJ/mol（0 K下），表明茶碱与胞嘧啶的沃森-克里克面结合在能量上更具优势。这一能量趋势在不同理论方法和基组下均保持一致。

3.2. THEOP·A和THEOP·C复合物的光谱表征

实验上，由于DNA/RNA碱基的激发态寿命短，研究利用茶碱作为发色团，通过REMPI光谱探测复合物。THEOP·A和THEOP·C的REMPI光谱呈现为无结构的宽吸收，与茶碱单体清晰的振动结构形成对比。通过在不同电离波长下获取IDIR光谱，并与理论模拟的IR谱图对比，成功指认了复合物中存在的构象。对于THEOP·C，实验IDIR光谱在3400 cm^-1以上显示出对应于自由NH/OH伸缩振动的离散吸收峰，在3400 cm^-1以下则是一个宽而强的吸收带，这是强耦合氢键的典型特征。光谱分析表明，至少有两种构象（THEOP·C-1和THEOP·C-2）共存。THEOP·C-1中未参与成键的胞嘧啶-NH₂基团的反伸缩模式在3535 cm^-1处产生强峰，而THEOP·C-2的自由-NH₂对称和反对称振动模式则分别对应于3460 cm^-1和3575 cm^-1处的吸收。对于THEOP·A，其IDIR光谱同样可清晰分为两个区域。最低能量构象THEOP·A-1准确地再现了3458 cm^-1和3566 cm^-1处的振动带，分别对应腺嘌呤自由-NH₂基团的对称和反对称N-H伸缩模式。其余谱带则与THEOP·A-2的预测吻合。光谱证据证实了计算预测的多种稳定构象在气相中同时存在。

3.3. THEOP·A和THEOP·C之间“栓-槽”相互作用的比较分析

对比THEOP·C和THEOP·A的结构指认发现，在无外界干扰的气相环境中，这些分子的聚集主要由强双氢键的形成驱动。尽管并非实验条件下最稳定的排列方式，但具有生物学相关性的相互作用（即茶碱与核碱基沃森-克里克面的相互作用）被证实存在，且与茶碱-适配体复合物中的相互作用模式一致。这凸显了分子识别的一个基本原则：结构互补性。如同栓精确插入对应的槽，茶碱分子的突出官能团（“栓”），特别是其O=C⁶和NH基团，与胞嘧啶和腺嘌呤沃森-克里克面上的互补氢键受体和供体（“槽”）完美契合。

研究进一步分析了实验光谱中不同构象特征谱带的相对强度。对于两个体系，光谱中均观察到分别对应于自由NH₂基团和氢键键合NH₂基团的两个特征谱带。这两个谱带的强度比可能源于多种因素，包括振子强度、实验条件（如红外光子吸收导致的复合物解离）以及构象的相对稳定性（种群）。分析指出，THEOP·C-1构象的信号强度相对更高，而这正是两个体系中更稳定的物种。结合理论计算显示THEOP·C-1比THEOP·A-2结合能更高，可以定性地得出结论：结合能与构象特异性信号的相对强度存在相关性。尽管茶碱与胞嘧啶或腺嘌呤结合时观察到相同类型的官能团相互作用，但THEOP·C-1复合物在自由环境中的更高丰度，为两种核碱基之间的选择性提供了证据。

总结与讨论部分指出，本研究通过质谱分辨激光光谱和DFT计算，表征了茶碱与胞嘧啶、腺嘌呤核碱基二聚体的形成。研究成功分离出THEOP·C和THEOP·A复合物各两种稳定构象，其形成均由非常强的NH···N/NH···O=C双氢键驱动。其中一种分离出的构象（THEOP·C-1和THEOP·A-2）呈现出与相应核碱基沃森-克里克面的相互作用，这与适配体中的结合模式一致。在孤立条件下，计算预测和振动谱带相对强度的实验趋势均表明，尽管氢键相互作用相似，但THEOP·C-1的结合能高于THEOP·A-2二聚体。这或许是RNA适配体对茶碱表现出高选择性和特异性的原因之一，因为其结合口袋中存在与胞嘧啶的相互作用。与DNA适配体-茶碱复合物结合口袋中相互作用的比较表明，即使与某个核碱基沃森-克里克面相关的结合能较低，通过周围环境中额外的非共价相互作用的协同效应以及结合腔的动态特性，仍然可以实现相当甚至更高的分子选择性。这些发现意味着，虽然强双氢键的形成是THEOP·C和THEOP·A复合物聚集的主要驱动力，但分子识别的高选择性不能仅用这些相互作用来解释。次要的非共价相互作用、协同效应以及结合口袋的动态特性都共同决定了最终的结合偏好性和稳定性。

本研究成功揭示了“栓-槽”相互作用的本质属性，表明茶碱-胞嘧啶（THEOP·C-1）这一相互作用模式在本质上更稳定、更容易形成，从而将相互作用强度确定为影响分子识别的变量之一。这项工作突出了孤立复合物光谱学与尖端量子化学方法结合，在阐明生物分子识别背后分子原理方面的强大能力。将此类研究扩展到更大体系，并系统比较所有可能的DNA/RNA碱基与黄嘌呤衍生物（包括茶碱）之间的相互作用，将为支配适配体分子识别的非共价作用力提供更深入的见解。