《Journal of Genetic Engineering and Biotechnology》：RNAi-based therapeutics targeting the F gene in human metapneumovirus: An
in-silico approach
编辑推荐：
背景：人偏肺病毒（Human metapneumovirus，HMPV）是急性呼吸道感染的主要病原之一，尤其在婴幼儿、老年人及免疫低下人群中危害显著。尽管临床负担沉重，目前尚无获批的抗病毒药物或疫苗。本研究旨在针对HMPV高度保守的融合（Fusion，F）基因
背景：人偏肺病毒（Human metapneumovirus，HMPV）是急性呼吸道感染的主要病原之一，尤其在婴幼儿、老年人及免疫低下人群中危害显著。尽管临床负担沉重，目前尚无获批的抗病毒药物或疫苗。本研究旨在针对HMPV高度保守的融合（Fusion，F）基因，识别具有潜力的RNA干扰（RNA interference，RNAi）治疗候选分子。方法：研究人员采用计算机模拟流程设计并评估靶向HMPV F基因的小干扰RNA（small interfering RNA，siRNA）。使用siDirect 2.1软件初步生成869条siRNA序列，依据序列特征与预测沉默效率逐级筛选至34条，再至10条候选；随后进行分子对接评估候选siRNA与人Argonaute-2（AGO2）蛋白的相互作用，并开展100 ns分子动力学模拟以评价结构稳定性，辅以主成分分析（principal component analysis，PCA）与自由能景观分析。结果：三条siRNA候选（F11、F19、F26）表现出良好的RNAi特征，包括高预测沉默效率与强靶标结合力。其中F11在全球HMPV毒株中序列保守性最高、熔解温度最高、AGO2互作有利且动态行为稳定；PCA显示F11采样的最受限构象子空间具有明确低能阱，被认定为最有前景候选。结论：该研究通过计算分析确定F11为靶向高度保守HMPV F基因的最有前景siRNA候选，具备作为RNAi治疗剂的较强潜力，但尚需体外与体内验证确认其效力与安全性。
人偏肺病毒（HMPV）属单股负链RNA病毒、肺病毒科，自2001年发现以来已被证实为全球婴幼儿及老年人群急性呼吸道感染的重要病原，近年自2024年末起病例数在多国回升，但目前临床仅能支持治疗，无获批抗病毒药或疫苗，核糖韦林与免疫球蛋白仅有轶事证据。鉴于HMPV基因组中F基因在A、B两谱系间氨基酸一致性约95%、且介导病毒包膜与细胞膜融合这一进入必需步骤，研究人员选择其作为广谱、耐用RNAi靶标，以弥补N、P基因作用于进入后复制阶段、G与SH基因变异性大不适于广谱设计的不足。论文发表于《Journal of Genetic Engineering and Biotechnology》，报告了首次针对HMPV F蛋白鉴定高效siRNA的计算框架。
为开展研究，研究人员从NCBI GenBank获取45株全球HMPV毒株F基因编码序列构建进化分析，使用siDirect 2.1依Ui-Tei、Reynolds、Amarzguioui经验规则初筛869条siRNA，以seed-target双链熔解温度（Tm）<15 °C、GC含量≥30%、折叠自由能≥0 kcal/mol、siRNAPred评分>0.8且siPRED抑制率>80%逐级过滤至F11、F19、F26三条；以BLASTN比对人转录组/基因组控脱靶；用RNAstructure测二级结构与双链ΔG，DINAMelt算Tm；HDOCK盲对接人AGO2（PDB 4Z4D）并用PDBsum析互作；GROMACS 2024.4跑100 ns分子动力学（AMBER99SB-ILDN力场、TIP3P水模型），后续做PCA与基于Boltzmann分布的自由能景观（free-energy landscape，FEL）。
研究结果如下。靶标序列获取与siRNA预测：从869条经经验规则留34条，Tm过滤留10条，GC、折叠能、效能过滤终得F11、F19、F26。脱靶控制与进化关系：10条候选BLASTN比对人均E值>0.05，F11在45株中覆盖保守区超50%，系统树bootstrap>70%支持F基因跨株保守。GC与二级结构：候选GC 30%–39%，F28因折叠ΔG=-1.3 kcal/mol排除，F11折叠ΔG=0 kcal/mol无自折叠劣势。热力学与Tm：双链ΔG为-30.0至-39.2 kcal/mol，F11为-38.9 kcal/mol；Tm(Conc.)与Tm(Cp)均>80 °C，F11达84.5 °C与85.5 °C。效能预测：siRNAPred中F11=0.965、F19=0.987、F26=0.988，siPRED抑制率均>92%。AGO2对接：HDOCK得分F11=379.76、F19=321.2、F26=420.66，F11定位于PAZ-MID-PIWI经典沟槽，与Gly178、Arg179、Arg710等保守残基成键。分子动力学：100 ns内F11 backbone RMSD Plateau约0.45–0.55 nm、RMSF均值0.224±0.119 nm为三者最低，Rg后期维持3.00–3.06 nm紧凑态；PCA中F11 PC1+PC2方差98.4 nm2为最小，FEL最低能阱13.7 kJ/mol，构象子空间最受限。讨论部分指出，F靶标兼顾进入必需性与跨谱系保守性，F11覆盖超50%全球株且具最高Tm与最稳AGO2复合体，RNAi平台已有patisiran等临床先例及SARS-CoV-2吸入式siRNA佐证呼吸道可行性；作者同时承认纯计算局限，未做MM-GBSA分解与重复模拟，脱靶与免疫刺激（TLR7/8、RIG-I）仅经序列 motif 排查，需体外 Reporter 敲降、qRT-PCR、空斑减少及动物模型跟进。结论翻译：RNA干扰（RNAi）是细胞天然机制，siRNA引导互补mRNA降解；本计算机模拟研究设计出F11、F19、F26三条靶向HMPV保守F蛋白的siRNA，F11凭热力学、沉默效率、AGO2互作与动态稳定性领先，但结论仅基于计算模型，须经体外体内实验验证方具治疗转化价值，本研究视为后续抗病毒设计的首步。
