综述：生物碱抗肠道病毒A71的疗效与机制：基于临床前证据的系统评价与荟萃分析联合网络药理学及分子对接研究

《International Journal of Molecular Sciences》：Efficacy and Mechanisms of Alkaloids Against Enterovirus A71: A Systematic Review and Meta-Analysis of Preclinical Evidence Integrated with Network Pharmacology and Molecular Docking

　　

肠道病毒A71（EV-A71）是手足口病的主要致病原，目前尚无获批的特异性抗病毒疗法。本研究通过整合荟萃分析、网络药理学及分子对接方法，系统评价生物碱抗EV-A71感染的疗效与作用机制。研究共纳入9项动物实验，荟萃分析结果显示，生物碱干预可显著提升感染动物生存率（比值比=30.62，95%置信区间：10.44–89.82），降低临床严重程度评分，减轻体重下降幅度，并减少感染组织内的病毒载量。亚组分析初步提示，喹嗪啶类生物碱及高剂量方案（>5 mg/kg）可能与临床前干预效应存在关联。网络药理学预测得到7种活性生物碱与EV-A71相关基因共有的155个交集靶点，经计算筛选确定丝裂原活化蛋白激酶1（MAPK1）、MAPK3、Jun原癌基因（JUN）、极光激酶B（AURKB）及MAPK8为核心靶点。功能富集分析显示，MAPK、肿瘤坏死因子（TNF）及白细胞介素-17（IL-17）信号通路显著富集。分子对接结果从计算层面支持活性生物碱与核心靶点间存在稳定结合亲和力（?6.7至?9.3 kcal/mol）。上述研究结果提示，生物碱可通过直接抗病毒作用及宿主导向调控机制发挥抗EV-A71效应，为其作为新型抗EV-A71治疗药物候选物提供了理论依据。

1. 1.
   引言
   肠道病毒A71（EV-A71）属于小RNA病毒科，为非包膜正义单链RNA病毒，基因组全长7.4 kb，由5'非翻译区（5' UTR）、单一开放阅读框及3'非翻译区（3' UTR）构成。其基因组通过内部核糖体进入位点（IRES）依赖的翻译机制合成多聚蛋白，经蛋白酶解切割后形成4种结构蛋白（VP1-VP4）及7种非结构蛋白（2A、2B、2C、3A、3B、3C、3Dpol）。EV-A71是手足口病（HFMD）的首要致病原，该病多为儿童自限性轻症感染，表现为发热、皮疹及口腔溃疡，但约1.1%的感染者会快速进展为无菌性脑膜炎（AM）、脑干脑炎（BSE）、急性弛缓性麻痹（AFP）、神经源性肺水肿（NPE）、心肺衰竭甚至死亡等严重神经系统及全身并发症。自1969年首次分离以来，EV-A71已在全球尤其是亚太地区引发多次大流行，构成持续的全球公共卫生威胁。目前国内获批的三价灭活EV-A71疫苗仅对同源毒株具有保护作用，无交叉保护效力，亦无法对已建立的感染产生治疗效果。广谱抗病毒药物利巴韦林抗EV-A71疗效欠佳且儿科不良反应显著，未被批准用于该适应证，现有治疗手段仍以对症支持为主，亟需开发安全有效的抗EV-A71药物。植物源天然产物是抗病毒药物发现的重要来源，其中生物碱因含氮杂环结构具备多样药理活性，已有多项临床前研究证实小檗碱、千金藤素、石蒜碱及其衍生物等天然生物碱在细胞及动物模型中具有显著抗EV-A71活性，作用机制涵盖靶向病毒入侵、抑制吸附与复制、阻断蛋白合成、调控宿主氧化应激与免疫反应及减轻组织损伤等，但现有研究多为小样本单中心临床前试验，疗效异质性较高，尚缺乏系统评价定量整合临床前证据，其多靶点作用机制亦未明确。本研究通过整合网络药理学与荟萃分析，系统评价生物碱抗EV-A71感染的疗效并阐释其多靶点作用机制，为生物碱类抗EV-A71治疗策略的开发与临床转化提供临床前证据与理论基础。
2. 2.
   方法
   2.1 荟萃分析
   本研究遵循PRISMA 2020声明开展，采用SYRCLE工具评估动物研究偏倚风险，研究方案已在PROSPERO注册（注册号：CRD420261336466）。
   2.1.1 检索策略
   检索数据库涵盖PubMed、Web of Science核心合集、Embase、Cochrane Library、Scopus、中国知网、维普数据库及万方数据，检索时限为建库至2026年3月，采用医学主题词（MeSH）与自由词结合的方式构建检索式，无语言及发表时间限制，同时手动筛查参考文献以补充潜在合格文献。
   2.1.2 纳入与排除标准
   基于PICOS框架制定纳入标准：研究对象为EV-A71感染的体内动物模型；干预措施为生物碱化合物干预；对照组采用溶剂、生理盐水、磷酸盐缓冲液（PBS）或无干预；结局指标包括主要结局生存率，次要结局临床评分、体重及病毒复制情况；研究设计为对照动物实验。排除标准包括：不符合纳入标准的研究、无法获取完整数据用于荟萃分析的研究、二次文献（综述、摘要等）及无任何预设相关结局指标的研究。
   2.1.3 数据提取与质量评估
   由两名研究者独立采用标准化表格提取数据，内容包括第一作者/发表年份、实验动物特征（物种、品系、年龄）、EV-A71感染模型构建细节（病毒株型、接种剂量、接种途径）、干预措施详情（生物碱类型及名称、给药剂量、频率、时长、途径）、对照组干预方案及结局指标。图表数据优先联系通讯作者获取原始数据，否则采用WebPlotDigitizer v4.7进行量化；组织病毒载量单位统一转换为log₁₀尺度，低于检测限（LOD）的病毒载量采用LOD/2进行插补。偏倚风险评估采用SYRCLE RoB工具，涵盖选择偏倚、实施偏倚、检测偏倚、 attrition偏倚、报告偏倚及其他偏倚，每个条目评为低风险、不清楚风险或高风险，分歧通过讨论或第三方仲裁解决。采用GRADE体系评估方法学质量。
   2.1.4 统计分析
   采用Review Manager 5.4进行数据分析与可视化。二分类结局采用Mantel–Haenszel法合并比值比（OR）及95%置信区间（CI）；连续变量采用均差（MD）及95% CI。P<0.05为差异有统计学意义。采用卡方检验与I²统计量评估异质性，P>0.1且I²<50%时采用固定效应模型，否则采用随机效应模型。预设亚组分析基于生物碱类型（异喹啉类vs喹嗪啶类）及剂量（≤5 mg/kg vs >5 mg/kg）探索异质性来源，采用Stata 17.0进行敏感性分析评估单个研究对合并效应量的影响。因纳入研究不足10项，未采用漏斗图与Egger检验，改用Orwin失安全系数法评估发表偏倚，设定无效效应量为OR=1.1，失安全系数越大表明结果受发表偏倚影响越小。

2.2 网络药理学

2.2.1 化合物靶点预测

经荟萃分析筛选得到小檗碱、千金藤素、吐根碱、哈敏、石蒜碱、苦参碱及氧化苦参碱7种化合物，因石蒜碱衍生物7e与LY-55在公共数据库中无标准化结构标识符与规范SMILES序列，故以原生石蒜碱作为替代化合物进行探索性分析。从PubChem数据库获取各化合物的规范SMILES序列，采用SwissTargetPrediction服务器预测人源潜在靶点，合并去重后获得非冗余靶点列表。

2.2.2 疾病相关靶点识别

检索DisGeNET与GeneCards数据库获取EV-A71感染相关疾病靶点，DisGeNET检索词为“Enterovirus A71”，筛选得分≥0.2的靶点；GeneCards检索词为“Enterovirus A71”“Enterovirus infection”“Hand, foot and mouth disease”，筛选相关性得分≥5的靶点，合并去重后得到非冗余EV-A71相关疾病靶点共1768个基因。

2.2.3 交集靶点识别

将化合物相关靶点与疾病相关靶点取交集，获得生物碱抗EV-A71感染的潜在治疗靶点，采用R软件VennDiagram包绘制韦恩图，共得到155个共同靶点用于后续蛋白-蛋白互作（PPI）网络构建与功能富集分析。

2.2.4 PPI网络构建

将共同靶点导入STRING数据库，设定物种为人，最小互作评分设为0.4，隐藏无连接节点，下载TSV格式互作数据，采用Cytoscape 3.10.4进行网络可视化，运用CytoHubba插件计算最大团中心性（MCC）、介数中心性（BC）、接近中心性（CC_net）及度中心性（DC）4种中心性指标以识别网络枢纽基因。

2.2.5 核心靶点筛选与网络分析

基于化合物-靶点互作网络筛选核心靶点：首先通过CytoHubba计算PPI网络的4种中心性指标，结合药理相关性过滤，筛选被至少2种化合物靶向且在网络中发挥核心作用的基因，最终确定MAPK1、MAPK3、JUN、AURKB及MAPK8为核心药理靶点，此处的度值为靶向该基因的活性化合物数量，而非PPI网络中的蛋白互作数量。采用Cytoscape的yFiles Layouts插件进行分层布局，节点大小与度值成正比，按节点类型（化合物、靶点、疾病、通路）区分颜色，构建化合物-靶点-通路网络并进行可视化。

2.2.6 功能富集分析

采用R软件clusterProfiler包（4.8.0版）对155个共同靶点进行基因本体（GO）功能富集分析，涵盖生物过程（BP）、细胞组分（CC_GO）及分子功能（MF）三类，设定物种为人，显著性阈值为P<0.05，采用Benjamini–Hochberg（BH）法进行多重检验校正，每类选取前10个富集条目绘制条形图。采用clusterProfiler包（4.10.0版）进行京都基因与基因组百科全书（KEGG）通路富集分析，选取前20个富集通路绘制气泡图。通路可视化优先选取同时满足以下标准的通路：调整后P<0.05、与EV-A71感染生物学相关（病毒感染、免疫应答、炎症或凋亡）及至少连接1个核心靶点。

2.3 分子对接验证

从RCSB蛋白质数据库获取核心靶蛋白MAPK1、MAPK3、JUN、AURKB、MAPK8的三维结构，从PubChem数据库下载活性化合物的3D结构。采用PyMOL 2.5去除蛋白结构中的水分子与原配体，添加氢原子与电荷，根据原配体位置或靶蛋白活性中心定义对接口袋。采用CB-Dock服务器进行分子对接，结合亲和力以kcal/mol为单位，≤?6.0 kcal/mol定义为稳定结合。采用PyMOL 2.5与Maestro 12.8进行3D结合构象与2D相互作用图的可视化分析，因氧化苦参碱仅靶向1个核心靶点（MAPK1）且预测亲和力（?5.8 kcal/mol）略高于?6.0 kcal/mol阈值，故未纳入对接分析。

1. 3.
   结果
   3.1 荟萃分析
   经数据库初检共获得227篇文献，去重后剩余148篇，经标题与摘要筛选后对21篇全文进行评估，最终9项研究符合纳入标准。
   3.1.1 研究选择与特征
   纳入研究发表于2011至2025年，均采用小鼠构建EV-A71感染动物模型，涵盖ICR小鼠、KM小鼠、AG129小鼠及BALB/c小鼠，年龄为出生24小时内新生鼠、7–14日龄乳鼠至5周至3月龄成年鼠。所有研究均通过腹腔注射诱导EV-A71感染，病毒株包括EV-A71 MP4、C4、MP10、GZ-CII、H-MA及Gluc-EV-A71等不同亚型，接种滴度为1×10⁶ 50%组织培养感染剂量（TCID₅₀）至1.2×10⁸空斑形成单位（PFU）。干预措施均为天然生物碱，分为异喹啉类生物碱（千金藤素、小檗碱、吐根碱、石蒜碱衍生物7e、石蒜碱衍生物LY-55、石蒜碱）、吲哚类生物碱（哈敏）及喹嗪啶类生物碱（氧化苦参碱、苦参碱），给药剂量范围为0.20 mg/kg至20 mg/kg，主要给药途径为腹腔注射，仅1项研究采用灌胃给药。干预方案包括感染前预防干预、感染后治疗干预及感染前联合持续感染后给药，干预时长为4至14天。所有研究均报告EV-A71感染相关结局，生存率、临床症状评分及病毒载量为最常评估的指标。
   3.1.2 偏倚风险评估
   SYRCLE偏倚风险评估显示，纳入研究的方法学报告完整性普遍不足：无研究详细描述随机序列生成的具体方法，仅4项研究明确报告组间基线特征可比性，所有研究均未报告分配隐藏、随机饲养及人员盲法；检测偏倚方面，仅2项研究报告随机结局评估，仅1项研究实施结局评估者盲法；所有研究均完整报告结局数据，仅1项研究存在选择性结局报告缺陷，无研究在任何领域被评为高风险。GRADE评估显示，主要结局的证据确定性为中等，次要结局因严重偏倚风险、显著异质性及间接性，证据确定性评级为低至极低，反映次要发现具有探索性与假设生成性质。
   3.1.3 定量合成（荟萃分析）
   需注意，合并分析涵盖了药理特性各异的不同类别生物碱，较大的OR需谨慎解读，其为跨多种化合物的聚合指标，不应被视为统一的类效应证据。生存分析纳入8项研究，结果显示生物碱干预较对照组显著提升EV-A71感染动物的生存率[OR=30.62，95%CI：10.44–89.82，P<0.00001]，异质性I²=0%，P=0.79，GRADE证据确定性为中等。该OR提示汇总分析中存在强统计学关联，但其效应幅度需谨慎解读，可能受纳入小鼠模型中样本量小及零事件研究的夸大影响，不可直接外推为人类EV-A71感染的临床疗效。临床结局分析纳入8项研究，结果显示生物碱干预可降低EV-A71感染动物的临床严重程度评分[MD=?2.86，95%CI：?3.55–?2.17，P<0.00001]，异质性I²=97%，P<0.00001，GRADE证据确定性为低；4项研究数据显示，生物碱干预可减轻EV-A71诱导的体重下降[MD=2.68，95%CI：2.00–3.37，P<0.00001]，异质性I²=89%，P<0.00001，GRADE证据确定性为极低。病毒复制分析纳入3项研究，结果显示生物碱干预可降低感染动物脑组织病毒载量[MD=?2.38，95%CI：?4.16–?0.60，P=0.009]，异质性I²=96%，P<0.00001，GRADE证据确定性为极低；6项研究数据显示，生物碱干预可降低后肢肌肉组织病毒载量[MD=?3.67，95%CI：?6.12–?1.22，P=0.003]，异质性I²=100%，P<0.00001，GRADE证据确定性为低；2项研究数据显示，生物碱干预对VP1蛋白表达无显著影响[MD=?4.93，95%CI：?14.28–4.41，P=0.30]，异质性I²=94%，P<0.00001，GRADE证据确定性为极低。
   3.1.4 亚组分析
   鉴于研究间存在高度异质性，为探索最优干预策略，按生物碱类型（异喹啉类vs喹嗪啶类）及剂量（低剂量≤5 mg/kg vs高剂量>5 mg/kg）进行预设亚组分析。结果显示，生物碱类型与剂量可部分解释临床评分与体重指标的异质性，喹嗪啶类生物碱与高剂量方案在减轻体重下降与降低临床严重程度方面，较异喹啉类生物碱与低剂量方案呈现更优疗效趋势。但亚组划分后所有次要结局仍存在显著异质性，提示动物品系、病毒亚型等其他未测量因素也可能导致变异。
   3.1.5 发表偏倚与敏感性分析
   敏感性分析通过逐一剔除单项研究评估合并结果的稳定性，所有结局的合并效应量及其95%CI均保持统计学显著性且未跨越无效线，证实整体结论稳健。因纳入研究仅9项，未达到漏斗图与Egger检验推荐的最低10项标准，故采用Orwin失安全系数法评估发表偏倚影响，设定临界效应量（dc）为1.1，计算得失安全系数Nfs=279，提示需要279项未检索到的OR=1.1的未发表研究才能将合并效应量（OR=30.62）降低至1.1的微小水平，表明本研究主要结论在一定程度上可抵抗潜在发表偏倚，但因纳入研究数量有限，仍不能完全排除发表偏倚的可能。

3.2 网络药理学分析

基于荟萃分析提示的7种生物碱潜在抗EV-A71疗效，开展网络药理学分析以探索其作用的多靶点机制。从GeneCards数据库获取1768个EV-A71相关疾病靶点，经SwissTargetPrediction预测得到444个化合物相关靶点，交集分析得到155个共同药物-疾病靶点。构建化合物-靶点-疾病网络显示，7种化合物通过155个共同靶点与EV-A71感染核心节点相连。构建靶点-通路网络揭示核心靶点与关键信号通路的关联，桑基图直观展示了7种化合物、4个核心靶点（JUN、MAPK1、MAPK3、AURKB）与10条关键通路的核心调控轴。通过4种中心性指标（MCC、BC、CC_net、DC）的交集筛选，确定MAPK1、MAPK3、JUN、AURKB及MAPK8为核心靶点，PPI网络显示MAPK1、MAPK3与JUN等枢纽基因间存在紧密互作。GO富集分析显示，生物过程主要富集于免疫应答、炎症反应与细胞增殖，分子功能主要与蛋白结合及激酶活性相关。KEGG通路富集分析显示，核心靶点主要参与MAPK信号通路、肿瘤坏死因子（TNF）信号通路、白细胞介素-17（IL-17）信号通路及癌症通路等与病毒感染及炎症密切相关的通路。

3.3 分子对接分析

为验证网络药理学预测的核心靶点与代表性生物碱的结合可行性，开展分子对接分析。小檗碱与MAPK1结合亲和力最强，达?9.3 kcal/mol，与活性口袋关键残基LYS-54及LYS-114形成稳定氢键，并伴有广泛疏水相互作用；小檗碱与MAPK3结合亲和力为?9.1 kcal/mol，与THR-185形成关键氢键，可紧密嵌入结合域；苦参碱与JUN结合亲和力为?6.8 kcal/mol，与ASP-464形成稳定氢键，并通过疏水接触匹配JUN的疏水空腔；石蒜碱与MAPK1稳定结合，亲和力为?7.6 kcal/mol，与ASP-179及SER-208形成氢键，增强活性位点内的极性相互作用；吐根碱与AURKB结合亲和力为?7.9 kcal/mol，与ATP结合口袋残基形成多重氢键与疏水相互作用，支持复合物稳定形成；哈敏与MAPK8结合亲和力为?6.7 kcal/mol，与ASN-114形成氢键及稳定疏水相互作用，确保牢固结合。2D相互作用图显示，所有化合物均可与靶蛋白口袋关键氨基酸残基形成氢键与疏水相互作用，为预测的结合互作提供了计算层面的结构支持。

1. 4.
   讨论
   4.1 证据总结
   本研究首次整合网络药理学、分子对接与临床前荟萃分析，系统评价生物碱抗EV-A71感染的干预效应。荟萃分析纳入9项研究，涵盖9种生物碱或其衍生物，核心证据表明天然生物碱干预可显著提升EV-A71感染动物的生存率（OR=30.62，95%CI：10.44–89.82，P<0.00001），该结局无异质性（I²=0%），GRADE证据确定性为中等。需注意，生存率的合并OR值较大，应谨慎解读，致死感染模型的小动物研究效应估计可能存在夸大，尤其当对照组出现零事件数据或接近完全死亡率时，该结果不应被解读为人类EV-A71感染的同等获益幅度，而应视为特定实验条件下的干预有效性证据。生物碱干预还可显著降低临床严重程度评分、减轻EV-A71诱导的体重下降、降低脑与后肢肌肉组织病毒载量，提示其在体内可有效抑制病毒复制。预设亚组分析初步显示喹嗪啶类生物碱与高剂量方案（>5 mg/kg）具有更优的保护效应，为未来药物开发提供了明确方向。多项次要结局存在显著异质性，可能源于小鼠品系、年龄、病毒株、感染模型、生物碱类型、剂量、给药途径及干预时机等差异，因此这些合并估计值应主要视为假设生成而非确证性证据。方法学质量评估显示纳入研究整体质量中等，随机化、盲法等关键方法学细节报告不足，次要结局证据确定性评级为低至中等。Orwin失安全系数分析（n=279）与敏感性分析证实，生物碱改善感染动物生存率的核心结论稳健，受发表偏倚影响较小。网络药理学预测得到7种生物碱与EV-A71感染共有的155个交集靶点，核心靶点（MAPK1、MAPK3、JUN、AURKB、MAPK8）富集于MAPK、TNF及IL-17信号通路，这些通路在病毒复制与宿主炎症反应中发挥关键作用。分子对接进一步预测活性化合物与核心靶点间存在稳定结合亲和力（?6.7至?9.3 kcal/mol），为假设的化合物-靶点互作提供了计算层面的结构支持，但这些靶点仍需通过基因敲低、通路抑制、蛋白表达分析及病毒复制检测等实验验证其功能相关性。

4.2 机制概述

4.2.1 靶向EV-A71生命周期的直接抗病毒机制

多种生物碱通过干扰EV-A71复制周期的不同阶段发挥直接抗病毒作用。双苄基异喹啉类生物碱千金藤素通过中和内溶酶体酸化阻断病毒进入，将病毒粒子隔离于溶酶体内，阻止脱壳，该机制作用于感染早期，不依赖于病毒结合或内化过程。吐根碱在感染后阶段通过抑制IRES驱动的病毒多聚蛋白翻译，直接靶向病毒翻译装置，发挥复制抑制作用。石蒜碱及其衍生物（原生石蒜碱、7e、LY-55）通过阻断多聚蛋白延伸干扰病毒蛋白合成；此外，LY-55通过抑制JNK磷酸化减弱病毒诱导的自噬，进而限制病毒复制。上述直接作用机制靶向病毒进入、翻译与复制环节，实现对EV-A71的强效抑制。

4.2.2 宿主导向的保护机制：抗炎与抗氧化效应

除直接抗病毒作用外，多种生物碱通过调控炎症反应、氧化应激与免疫功能发挥宿主保护作用。小檗碱激活Keap1-Nrf2抗氧化通路，减少活性氧（ROS）生成，减轻EV-A71诱导的神经损伤；哈敏抑制核因子κB（NF-κB）活化，降低ROS生成，从而抑制病毒复制并减轻炎症反应，两种化合物均可降低感染组织中肿瘤坏死因子-α（TNF-α）、白细胞介素-1β（IL-1β）、白细胞介素-6（IL-6）等促炎细胞因子水平。喹嗪啶类的苦参碱与氧化苦参碱通过恢复T细胞群（CD3⁺、CD4⁺、CD8⁺）增强适应性免疫，提升生存率。这些宿主导向机制通过减轻组织损伤、控制过度炎症及促进免疫介导的病毒清除，与直接抗病毒效应产生协同作用。

4.2.3 网络药理学与分子对接预测的机制

网络药理学预测7种生物碱（小檗碱、千金藤素、吐根碱、哈敏、石蒜碱、苦参碱、氧化苦参碱）与EV-A71感染共享155个重叠靶点，核心靶点包括MAPK1、MAPK3、JUN、AURKB及MAPK8，显著富集于MAPK、TNF及IL-17信号通路。MAPK通路作为细胞增殖、分化与炎症的关键调控因子，与LY-55抑制c-Jun氨基末端激酶（JNK，MAPK家族成员）磷酸化的实验证据相印证；TNF与IL-17通路与小檗碱、哈敏干预后促炎细胞因子降低的观察结果一致。分子对接预测代表性生物碱与核心靶点间具有良好的结合亲和力，为这些互作的结构可行性提供了计算证据，但上述发现应视为假设生成而非机制验证，需通过实验研究确认。值得注意的是，AURKB（极光激酶B）在EV-A71感染中的实验表征仍属空白，为未来研究提供了新的假设方向。计算预测与实验证据的一致性凸显了生物碱抗EV-A71效应的多靶点、多通路特征。

4.3 优势与局限性

本研究首次整合网络药理学分析与临床前证据荟萃分析，全面评价生物碱抗EV-A71感染的干预效应，该方法学不仅实现了疗效数据的定量合成，还促进了潜在多靶点机制的系统性预测，衔接了计算预测与实验证据。分子对接的加入进一步探索了活性化合物与核心靶点（MAPK1、MAPK3、JUN、AURKB、MAPK8）的结合亲和力，为预测的化合物-靶点互作提供了结构支持。研究系统纳入9项合格动物实验，涵盖异喹啉类、吲哚类及喹嗪啶类7种核心生物碱，实现了不同生物碱亚型抗EV-A71证据的全面整合。基于生物碱类型与剂量的预设亚组分析探索了异质性的潜在来源，为最优干预策略提供了见解。研究严格遵循PRISMA指南，采用SYRCLE工具与GRADE体系分别评估方法学质量与证据确定性，并在PROSPERO注册研究方案，同时对关键结局开展发表偏倚与敏感性分析，提升了研究结果的稳健性与可靠性。

本研究存在若干局限性：首先，网络药理学结果为基于公共数据库的计算预测，仅能提供合理的机制假设而非直接实验证据；且网络药理学分析采用人源靶点数据库，而荟萃合成的疗效证据完全来自小鼠模型，尽管多数信号通路在人鼠间进化保守，但物种不匹配引入了间接性来源，可能影响预测靶点的生物学相关性，需通过适当的动物与人源系统进行实验验证。此外，两种半合成石蒜碱衍生物（7e与LY-55）因公共数据库中缺乏规范SMILES标识符，被排除在网络药理学分析之外。已识别核心靶点与通路的作用需通过体外与体内功能实验进一步验证，预测结果受限于底层数据库的完整性与准确性。其次，分子对接基于刚性对接模拟，无法完全捕捉构象动态变化及细胞环境对结合的影响，可能降低结合亲和力预测的准确性。网络药理学与分子对接均为基于可用数据库与静态蛋白结构的计算方法，未考虑化合物代谢、药代动力学行为、蛋白构象动态及体内复杂生物环境，因此预测的化合物-靶点互作应视为未来实验研究的假设而非确证性作用机制。第三，纳入的9项研究均为临床前动物实验，整体方法学质量中等，多数研究未充分报告随机化、盲法等关键细节，可能导致效应量高估，整体证据确定性评级为低至中等。第四，多项核心结局存在显著统计异质性，虽通过亚组分析识别了部分异质性来源，但受限于纳入研究数量（n=9），预设亚组分析仅能对每类至少包含两项研究的变量（生物碱类型与剂量）进行探索，基于动物品系、给药时机或病毒株的亚组分析因数据稀疏而无法开展。次要结局的高度未解析异质性意味着其合并估计值为假设生成性质。另一关键 caveat 为纳入生物碱的药理多样性：尽管按生物碱类别（异喹啉类vs喹嗪啶类）开展了亚组分析，但单个化合物的研究数量有限（n=1–4），无法进行稳健的单化合物荟萃分析，因此合并效应量不应被解读为所有生物碱具有相同疗效，而应视为对不同化合物群体的总体信号探索。未来研究应在标准化实验条件下开展不同生物碱的头对头比较，以实现直接疗效比较并识别最具潜力的先导化合物。干预时机也是重要的临床异质性来源，纳入研究涵盖病毒攻击前的预防性干预、感染后早期干预及联合干预方案，代表不同的生物学与临床场景，因此合并估计值不应仅被解读为症状发作后的治疗疗效证据，而应视为当前可用临床前干预研究的整体合成。此外，研究数量不足10项，未进行meta回归分析，无法全面探索其他潜在混杂因素的影响，观察到的异质性降低了合并效应估计的确定性。第五，所有研究均在小鼠模型中开展，绝大多数采用腹腔给药，与临床场景存在显著差异，限制了转化相关性。纳入研究多采用10–20 mg/kg腹腔给药，与潜在临床给药途径差异显著，无法直接外推至人类暴露水平，药代动力学特征、剂量换算、口服生物利用度、代谢及安全性谱需在临床转化前进一步研究。最后，尽管研究结果提示生物碱抗EV-A71感染的疗效，但纳入研究的局限性凸显了开展更严谨的临床前研究以强化潜在临床转化证据的必要性。

1. 5.
   结论
   本研究首次整合荟萃分析、网络药理学与分子对接，提供了探索性证据以评价生物碱抗EV-A71感染的潜在干预效应。对9项临床前研究的荟萃分析提示，生物碱干预可显著提升生存率（OR=30.62），降低临床严重程度，减轻体重下降，并减少脑与后肢肌肉组织病毒载量，主要结局的GRADE证据确定性为中等，且对发表偏倚具有稳健性。亚组分析确定喹嗪啶类生物碱（苦参碱、氧化苦参碱）与高剂量方案（>5 mg/kg）在减轻体重下降与临床症状方面更具优势。合并荟萃分析结果提示不同生物碱群体中存在有前景的疗效信号，但这些化合物的药理多样性排除了统一类效应的确证性结论，结果应视为假设生成，为未来特定先导化合物的筛选提供方向。
   基于现有实验文献与本研究的计算预测，研究人员提出以下多机制假设以供未来验证：直接抗病毒效应包括靶向病毒进入（千金藤素）、IRES驱动翻译（吐根碱）、多聚蛋白延伸/自噬（石蒜碱及其衍生物）；宿主导向保护效应包括Nrf2介导的抗氧化通路（小檗碱）、NF-κB抑制（哈敏）及T细胞恢复（苦参碱、氧化苦参碱）。网络药理学计算预测得到155个共同靶点，核心靶点（MAPK1、MAPK3、MAPK8、JUN）及信号通路（MAPK、TNF通路）为未来机制分析实验提供了明确方向。预测得到的计算结合亲和力（?6.7至?9.3 kcal/mol）为假设的机制轴提供了计算支持，但所有计算预测均应视为假设生成性质，需通过严谨的实验分析方可得出确证性的作用机制结论。
   本研究填补了天然产物抗EV-A71研究缺乏定量系统评价的空白，为抗病毒天然产物领域的未来荟萃分析建立了标准化方法学框架；确定的核心靶点（MAPK1、MAPK3、MAPK8、JUN）与信号通路（MAPK、TNF通路）为后续机制实验提供了清晰方向；不同生物碱亚型的疗效差异为新型抗EV-A71药物的先导化合物筛选与优化提供了重要线索。