基于三维结构虚拟筛选和分子动力学方法靶向冠状病毒核衣壳蛋白的广谱抑制剂发现

《COVID》：Broad-Spectrum Inhibitor Discovery Targeting Coronavirus Nucleocapsid Proteins via 3D Structure-Based Virtual Screening and Molecular Dynamics Ebtisam Aldaais, Munthir Aldukhi, Hind Alotaibi, Heba Mofleh Alzabni, Subha Yegnaswamy and Nada F. Alahmady

　　

冠状病毒刺突蛋白的快速抗原漂移促使研究人员探索替代抗病毒策略。研究人员靶向保守的核衣壳核衣壳（N）蛋白——该蛋白在RNA结合、基因组包装和复制中发挥核心作用——并开展了一项比较性、跨物种的基于三维结构的计算机模拟评估。一个包含494种化合物（天然、植物化学、合成）的文库被用于通过AutoDock Vina对MERS-CoV的N端RNA结合结构域（NTD；PDB 7DYD）以及SARS-CoV（2CJR）和SARS-CoV-2（8R6E）的C端二聚化结构域（CTD）进行分子对接，这反映了每个病毒均具有高分辨率、功能相关的结构域结构可用性。排名靠前的构象进行了ADME（吸收、分布、代谢、排泄）分析和100 ns的GROMACS分子动力学（MD）模拟。杨梅素-3-O-β-D-吡喃半乳糖苷（myricetin）在所有靶标中显示出最有利的预测对接分数（MERS–NTD：?8.9 kcal/mol；SARS–CTD：?10.1；SARS-CoV-2 CTD：?9.8）。姜黄素（curcumin）显示出中等预测亲和力（?7.1至?8.1），而MCC950则始终取得有利的对接分数（?7.9至?9.0）。ADME结果揭示了一种权衡：糖基化黄酮类化合物提供了丰富的相互作用网络，但违反了口服类药性标准（例如，高拓扑极性表面积TPSA），而MCC950符合Lipinski/Veber准则，支持其转化潜力。MD分析揭示了配体和靶标特异性的稳定性：杨梅素在SARS-CoV-2 CTD中保持持续结合超过100 ns，其根均方偏差（RMSD）低于对照；姜黄素在MERS和SARS结合复合物中表现出短暂稳定性（约30 ns）；MCC950显示出间歇性相互作用。总体而言，这些发现表明，保守的N蛋白RNA结合沟是一个对耐药具有弹性的广谱抗病毒发现靶点。天然黄酮类化合物为优化提供了有前景的骨架，而MCC950因其类药特性值得进一步探索。由于本研究纯粹是计算性的，结果属于假设生成性质，应通过RNA结合破坏测定、抗病毒细胞研究和体内模型进行验证。

## 论文解读：基于三维结构虚拟筛选与分子动力学发现靶向冠状病毒核衣壳蛋白的广谱抑制剂

### 研究背景与问题

冠状病毒，包括SARS-CoV、MERS-CoV和SARS-CoV-2，引发了严重的全球公共卫生危机。现有抗病毒策略多聚焦于刺突蛋白（Spike, S），然而该蛋白突变率高，导致疫苗保护效力下降和抗体逃逸。相比之下，核衣壳（Nucleocapsid, N）蛋白在冠状病毒中高度保守，负责RNA结合、基因组包装和复制调控，且具有较低突变趋势，因此成为开发广谱抗冠状病毒药物的理想靶点。目前尚缺乏系统比较SARS-CoV、MERS-CoV和SARS-CoV-2三者的N蛋白结构域对抑制剂结合响应的研究，多数研究仅评估了有限化合物库。针对这一空白，研究人员开展了一项跨物种、基于三维结构的虚拟筛选研究，并结合分子动力学模拟和药代动力学评估，旨在发现具有广谱活性的N蛋白抑制剂，为后续实验验证提供假说基础。该论文发表于《COVID》期刊。

### 主要关键技术方法

本研究为纯计算性工作，未涉及样本队列。关键技术方法包括：1) **蛋白结构准备**：从RCSB PDB获取MERS-CoV N端结构域（NTD, PDB: 7DYD）、SARS-CoV C端结构域（CTD, PDB: 2CJR）和SARS-CoV-2 CTD（PDB: 8R6E）的高分辨率晶体结构，去除结晶水、离子及非必需分子，添加氢原子，处理替代构象。2) **配体文库构建**：汇集494种化合物（462种天然营养品、18种植物化学物、14种合成化合物），从PubChem获取三维结构并能量优化。3) **分子对接与验证**：采用AutoDock Vina进行对接，对7DYD的共结晶配体进行盲对接验证（RMSD=2.94 ?）。4) **ADME分析**：使用基于RDKit的Python工作流及SwissADME服务器评估药代动力学性质，以Lipinski五规则和Veber标准评估口服生物利用度。5) **分子动力学模拟**：对选定的三种配体（杨梅素、姜黄素、MCC950）与各靶标复合物进行100 ns的GROMACS模拟，采用CHARMM36力场和SPC216水模型。6) **MM-PBSA分析**：使用gmx_MMPBSA计算70-100 ns稳定轨迹的结合自由能。

### 研究结果

**3.1 靶蛋白的结构与理化特性**
通过分析三个靶蛋白的晶体结构和理化参数（分子量、等电点、疏水性等），发现SARS-CoV-2 CTD（8R6E）具有最低的不稳定指数（30.7）和最高的脂肪族指数（53.9），表明其结构更刚性的疏水核心可能有助于抑制剂选择性结合二聚化界面。MERS-CoV NTD具有亲水特性（GRAVY -0.63），而SARS-CoV CTD亲水性更强（GRAVY -0.912）。

**3.2 对接方案验证**
对7DYD的共结晶配体进行重对接，得到RMSD为2.94 ?，在盲对接可接受范围内，表明对接参数有效；2CJR和8R6E因无共结晶配体，采用已验证参数统一应用。

**3.3 虚拟筛选与ADME分析**
494种化合物对接后，植物化学物Tetrandrine在三个靶标中均显示最高亲和力（-11.5至-14.8 kcal/mol），但后续综合ADME和配体多样性，选出代表性配体：杨梅素（天然）、姜黄素（植物化学）、MCC950（合成）。ADME显示杨梅素亲和力最佳（-8.9至-10.1 kcal/mol）但TPSA高（230.74 ?2）违反口服类药性；姜黄素亲和力中等（-7.1至-8.1），参数适中；MCC950持续有利亲和力（-7.9至-9.0）且符合Lipinski/Veber准则，支持转化潜力。

**3.4 配体-蛋白相互作用**
对接分析显示：杨梅素与三个靶标均形成广泛氢键（如Tyr99、Thr264、Arg13）和π-π/π-烷基相互作用，平均得分最高（-9.6 kcal/mol）。MCC950依赖芳香堆积和疏水接触（Trp43、Phe275），平均得分-8.5 kcal/mol。姜黄素氢键有限且存在供体-供体冲突（如SARS-CoV-2 CTD的Arg13），平均得分最低（-7.46 kcal/mol）。

**3.5 分子动力学模拟 - RMSD分析**
100 ns RMSD分析表明：杨梅素在MERS-CoV NTD（7DYD）和SARS-CoV-2 CTD（8R6E）中稳定（0.4-0.5 nm），但在SARS-CoV CTD（2CJR）中波动较大（1.5 nm）。MCC950在所有靶标中均维持低RMSD（0.4-0.7 nm），显示最高稳定性。姜黄素在8R6E中RMSD较高（1.0 nm），稳定性最差。

**3.6 后MD评估（RMSF、Rg、SASA）**
- **RMSF**：杨梅素在7DYD和2CJR中平均波动低（0.2-0.25 nm），但存在局部尖峰；姜黄素波动最高（0.3-0.5 nm）；MCC950在8R6E中波动最低（0.2 nm），表现最佳。
- **Rg（回转半径）**：MCC950在所有复合物中展现最佳紧凑性（1.47-1.51 nm）；杨梅素中等；姜黄素诱导最大膨胀（8R6E中达1.90 nm），稳定性差。
- **SASA（溶剂可及表面积）**：MCC950最小化溶剂暴露（82-85 nm2），杨梅素中等（81-88 nm2），姜黄素最高（83-93 nm2），反映其较差的结构稳定性。

**3.7 MM-PBSA分析**
结合自由能分解显示，所有配体主要依赖范德华力和非极性溶剂化能驱动结合，极性溶剂化能对结合不利。MCC950在多个复合物中具有最有利的结合自由能趋势，与动力学稳定性一致。

### 讨论与结论

讨论部分指出，N蛋白结构域的保守性与域特异性差异（如疏水性、柔韧性）支持开发泛冠状病毒抑制剂，但需针对不同谱系进行优化。虚拟筛选揭示了植物化学物（如粉防己碱）和天然多酚（如安石榴苷）具有多靶标强亲和力，而MCC950在类药性和动力学稳定性上表现突出。MD模拟强调，杨梅素提供强初始结合但存在靶标依赖适应性，MCC950则赋予长期稳定性和热力学优势，姜黄素因持续高波动而不适合作为直接N蛋白抑制剂。

**结论翻译**：这项基于结构的虚拟筛选研究强调了核衣壳（N）蛋白作为SARS-CoV、MERS-CoV和SARS-CoV-2中一个具有前景但尚未充分探索的治疗靶点，为广谱冠状病毒疗法的开发奠定了基础。研究发现，天然黄酮类化合物，特别是杨梅素，展现出卓越的预测结合亲和力（对接分数-8.9至-10.1 kcal/mol），并在整个100 ns分子动力学模拟中于SARS-CoV-2 C端结构域（CTD）内保持持续相互作用，表现出独特稳定性。相比之下，姜黄素在MERS-CoV和SARS-CoV复合物中仅维持约30 ns的短暂稳定，而大多数其他配体更早解离，这一现象可能归因于N蛋白结构域的高内在可动性。同时，合成化合物MCC950成为稳健候选物，兼具有利对接分数（-7.9至-9.0 kcal/mol）和类药药代动力学性质（符合Lipinski五规则和Veber生物利用度参数），其在炎症性疾病临床试验中的前期评价进一步增强其转化潜力。高亲和力化合物在所有三种冠状病毒中一致靶向保守的RNA结合沟，支持其作为广谱可成药口袋的作用，关键相互作用包括与带正电残基（如Arg13、Arg16、Arg92、Lys61）的氢键以及与芳香残基（如Tyr99、Tyr110、Pro12）的π-π堆积，提示竞争性抑制机制可阻断RNA结合和病毒组装。这些残基在冠状病毒间的保守性表明，N蛋白导向的抑制剂可能提供耐药弹性的抗病毒策略，尤其与受快速突变漂移限制的刺突蛋白靶向疗法相比。因此，这些发现倡导双轨发现策略：通过结构优化推进杨梅素和姜黄素等天然骨架以克服其药代动力学局限（如高拓扑极性表面积或次优生物利用度），同时利用MCC950等合成分子因有利ADME特征和已建立的安性数据实现快速转化进程。姜黄素、黄酮类化合物和MCC950的临床前及临床安全性数据增强了它们重定位的潜力，其多药理学能力可同时靶向病毒蛋白并调节宿主免疫反应。总体而言，这些结果为后续实验验证（包括确认N蛋白抑制的体外测定和评估疗效与生物利用度的体内研究）及药物化学优化以增强结合持久性和药代动力学性质提供了坚实基础。通过推进这些方向，本研究为加速开发广谱冠状病毒疗法铺平了道路，通过计算见解和转化策略结合来应对未满足的紧急医疗需求。