《BioMed Research International》：Predictive Analysis of Brain-Derived Neurotrophic Factor and Apolipoprotein E SNPs in Alzheimer’s Pathogenesis
编辑推荐：
摘要：脑源性神经营养因子（Brain-derived neurotrophic factor，BDNF）与载脂蛋白E（apolipoprotein E，APOE）是调控神经元功能及认知健康的关键分子，其基因变异尤其是非同义单核苷酸多态性（nonsynonymo
摘要：脑源性神经营养因子（Brain-derived neurotrophic factor，BDNF）与载脂蛋白E（apolipoprotein E，APOE）是调控神经元功能及认知健康的关键分子，其基因变异尤其是非同义单核苷酸多态性（nonsynonymous single-nucleotide polymorphisms，nsSNPs）已被证实与阿尔茨海默病（Alzheimer’s disease，AD）存在关联。研究人员采用计算生物学方法，对BDNF与APOE的nsSNPs潜在功能影响进行预测，以探索其在AD发病中的作用机制。研究从dbSNP数据库共检索到3590个BDNF SNPs与27830个APOE SNPs，经编码区定位与次要等位基因频率（minor allele frequency，MAF）≥0.001的质量筛选后，对33个BDNF nsSNPs与95个APOE nsSNPs开展系统性分析。致病性评估采用SIFT与PolyPhen-2算法，功能影响通过CADD评分进行评价，蛋白质稳定性效应使用MUpro与I-Mutant工具预测，翻译后修饰通过GPS预测系统分析，二级结构改变采用GOR4评估，三维结构模型通过SWISS-MODEL构建并经拉氏图（Ramachandran plot）验证。结果显示，三个变异在两种算法中均呈现一致的致病性预测结果：rs1048218（BDNF Q75H）、rs7412（APOE R176C）与rs769455（APOE R163C）。蛋白质稳定性分析表明，rs1048218（ΔΔG：-1.001至-2.08 kcal/mol）与rs7412（ΔΔG：-0.859至-0.07 kcal/mol）均表现出一致的去稳定效应，而rs769455的预测结果在不同算法间存在差异。所有变异的翻译后修饰位点均保持保守，二级结构分析显示α螺旋轻微减少（0.31%–0.81%），随机卷曲相应增加，三维建模显示尽管存在局部结构扰动，整体蛋白折叠得以保留，模型质量指标符合要求（MolProbity评分≤1.39，拉氏图有利区域占比>91%）。计算分析表明，BDNF与APOE的部分nsSNPs可能对蛋白功能与稳定性产生负面影响，即便其结构与翻译后修饰特征保持完整。这些计算预测结果仍需进一步的实验验证，以深入理解其在AD发病机制中的作用。
阿尔茨海默病（Alzheimer’s disease，AD）是一种以脑细胞损伤为核心的神经退行性疾病，也是最常见的痴呆类型，主要影响60岁及以上人群。其核心病理特征包括β淀粉样蛋白（amyloid-β，Aβ）斑块沉积与含tau蛋白的神经原纤维缠结（neurofibrillary tangles，NFTs），导致记忆丧失、认知下降、神经炎症及神经元慢性丢失。AD的发病受显著遗传因素影响，其中脑源性神经营养因子（brain-derived neurotrophic factor，BDNF）与载脂蛋白E（apolipoprotein E，APOE）被确认为两大关键遗传决定因子。早发性AD多与淀粉样前体蛋白（amyloid precursor protein，APP）、早老素1（presenilin 1，PSEN1）及早老素2（PSEN2）基因突变相关，而迟发性AD（约占95%）的主要遗传风险来自APOE ε4等位基因。BDNF支持神经元存活、分化及突触重塑，其水平降低与Aβ积累、tau磷酸化及神经元死亡密切相关；APOE则在脂质转运与稳态中发挥关键作用，其氨基酸替换可改变与受体、脂质及Aβ的结合亲和力，从而影响AD进程。然而，BDNF与APOE的具体功能受损机制尚未完全阐明，非同义单核苷酸多态性（nonsynonymous single-nucleotide polymorphisms，nsSNPs）作为可导致氨基酸替换的重要遗传变异，其对蛋白结构与功能的影响亟待解析。在此背景下，研究人员开展了一项计算生物学研究，旨在预测BDNF与APOE中nsSNPs的功能影响及其在AD发病中的潜在作用，相关成果发表于《BioMed Research International》。
研究采用的主要关键技术方法包括：从dbSNP数据库靶向检索BDNF与APOE的SNP数据，经编码区定位及MAF≥0.001的纳入标准与内含子、同义SNP及MAF<0.001的排除标准筛选后，获得高质量nsSNPs数据集；使用SIFT与PolyPhen-2算法进行致病性预测并筛选一致性高危变异；通过CADD框架评估功能影响；利用MUpro与I-Mutant预测蛋白稳定性变化；借助GPS工具分析翻译后修饰位点的保守性；采用GOR4算法预测二级结构变化；通过SWISS-MODEL构建三维同源模型并经拉氏图验证模型质量。
研究结果如下：
功能性影响评估通过CADD评分完成，筛选标准为PHRED评分≥20。结果显示，10个BDNF变异评分范围为20.4至31.0，12个APOE变异评分范围为23.3至56.0，提示其具有较强的有害功能影响潜力。
非同义SNP鉴定与致病性预测阶段，从原始数据库中筛选出505个BDNF nsSNPs与375个APOE nsSNPs，经严格筛选后最终纳入33个BDNF nsSNPs与95个APOE nsSNPs进行分析。SIFT预测BDNF有4个、APOE有17个nsSNPs为有害变异，PolyPhen-2则将3个BDNF变异与4个APOE变异归类为可能有害。其中，rs1048218（BDNF Q75H）、rs7412（APOE R176C）与rs769455（APOE R163C）在两种算法中均表现出一致的高致病性预测结果。
蛋白质稳定性分析显示，rs1048218（BDNF Q75H）在两个预测工具中均表现出显著的去稳定效应（ΔΔG：-1.001至-2.08 kcal/mol）；rs7412（APOE R176C）亦呈去稳定趋势，但I-Mutant预测效应较温和（ΔΔG：-0.07 kcal/mol）；rs769455（APOE R163C）在两个工具中的预测结果存在冲突，提示其效应可能为中性或微弱。
翻译后修饰分析通过GPS-MSP完成，结果表明所有已鉴定的磷酸化位点在野生型与各变异蛋白序列中均保持高度保守，未因突变而发生丢失，提示翻译后调控机制未受显著影响。
二级结构影响评估采用GOR4算法，结果显示各变异仅引起极轻微的α螺旋减少（0.31%–0.81%），伴随随机卷曲的代偿性增加，未出现新的二级结构元件，整体蛋白骨架得以保留。
三维结构建模与验证通过SWISS-MODEL完成，模型质量指标良好（MolProbity评分≤1.39），拉氏图分析显示APOE变异的立体化学质量优于BDNF变异，但所有模型均符合结构可靠性阈值，突变主要引起局部构象扰动而非全局折叠破坏。
讨论部分指出，rs1048218（BDNF Q75H）的谷氨酰胺替换为组氨酸可能破坏局部氢键网络，导致BDNF稳定性下降，进而影响其神经营养信号传递，这与AD患者中BDNF水平降低的临床观察相符。rs7412（APOE ε2等位基因）虽在计算预测中被判定为有害，但流行病学证据显示其为保护性因素，这可能与该变异改变了APOE与受体的结合特性，从而增强Aβ清除或减少tau病变有关，体现了单一算法预测的局限性。rs769455（APOE R163C）位于APOE结构域铰链区，其稳定性预测冲突反映了该变异效应的复杂性，近期研究提示其可能是非洲血统人群中的AD风险修饰因子。研究同时强调，这些变异对蛋白整体折叠影响有限，提示蛋白结构具有较强韧性，这为针对BDNF或APOE通路的精准治疗提供了可能。
结论部分总结认为，研究人员通过多工具整合分析，确定了rs1048218（BDNF）、rs7412与rs769455（APOE）三个高置信度潜在致病性变异。这些变异在多种计算流程中均被优先筛选出来，表现为一致的致病性预测、不同程度的热力学去稳定效应、翻译后修饰位点保守、仅引起细微结构扰动且模型质量可靠。这些发现为后续开展功能验证研究提供了靶点，有助于深入解析BDNF介导的神经营养信号通路与APOE介导的脂质代谢通路在AD中的作用，为精准医学策略的制定奠定基础。
