运用计算生物学技术探索ACADM基因有害非同义单核苷酸多态性及其在中链酰基辅酶A脱氢酶缺乏症中的作用

《Genetics Research》：Exploring Deleterious Nonsynonymous SNPs in the ACADM Gene: Insights Into Medium-Chain Acyl-CoA Dehydrogenase Deficiency (MCADD) via In Silico Analysis

这篇计算生物学综述通过综合生物信息学分析，系统探究了中链酰基辅酶A脱氢酶（MCAD）编码基因ACADM中非同义单核苷酸多态性（nsSNP）的功能影响。文章运用SIFT、PolyPhen-2等多种工具筛选出对蛋白质结构和稳定性有潜在破坏作用的nsSNP，并通过分子对接、分子动力学模拟等手段，重点揭示了rs200724875（R206H）和rs121434282（R281T）两个关键突变如何通过削弱ACADM蛋白与电子转移黄素蛋白（ETF）的相互作用，影响脂肪酸β-氧化，从而阐释了其导致MCADD的潜在分子机制，为精准医疗提供了靶点线索。

　　

¹. 引言

中链酰基辅酶A脱氢酶缺乏症（MCADD）是一种常染色体隐性遗传的潜在致死性代谢疾病，与ACADM基因的致病性突变密切相关。该基因编码的MCAD酶是脂肪酸β-氧化的关键催化剂，在空腹或代谢应激时为机体供能。功能丧失性突变会扰乱此过程，导致能量短缺，引发低酮性低血糖、呕吐、嗜睡、癫痫甚至猝死等临床症状。单核苷酸多态性（SNP）是人类基因组变异的主要形式，其中非同义SNP（nsSNP）能直接改变蛋白质的氨基酸序列，从而深刻影响其结构与功能，并与个体疾病易感性相关。本研究旨在通过计算生物学（in silico）方法，系统分析ACADM基因中的nsSNP，预测并理解它们对蛋白质功能的潜在影响，从而阐明遗传变异与代谢途径之间的复杂关系，为精准医疗提供分子层面的见解。

². 方法论

研究采用了多步骤的综合计算分析流程。首先，从dbSNP和Ensembl数据库获取了ACADM基因的所有nsSNP，经过去重后得到待分析的变异列表。随后，使用SIFT、PolyPhen-2、Meta-SNP、PhD-SNP、PANTHER、SNAP和MutationAssessor这七种生物信息学工具对获取的nsSNP进行扫描，筛选出被所有工具一致预测为有害的变异。对这些潜在有害nsSNP，进一步利用MutPred 1.2评估其对蛋白质结构和功能的影响，并使用I-Mutant 2.0分析其对蛋白质稳定性的影响。通过DeepREx-WS进行进化保守性分析，以确定突变位点的重要性。

为深入理解突变的结构后果，研究利用Robetta服务器对筛选出的最有害nsSNP（及野生型）进行了三维结构建模，并通过TM-align计算模型与野生型结构之间的均方根偏差（RMSD）和TM分数，以量化结构差异。其中，结构偏差最大的两个突变（R206H和R281T）还使用AlphaFold2进行了重新建模，并利用PyMOL进行可视化展示。

考虑到ACADM的功能依赖于与电子转移黄素蛋白（ETF）的相互作用以完成电子传递，研究对野生型及突变体ACADM与ETF（PDB ID: 1EFV）进行了分子对接分析，使用ClusPro v2.0评估结合能和簇大小，并通过PDBsum服务器分析复合物间的氢键、盐桥和非键接触等相互作用细节。为验证对接结果并评估复合物的动态稳定性，对野生型、R206H和R281T突变体分别与ETF形成的复合物进行了时长为200纳秒的分子动力学（MD）模拟。模拟使用Amber22软件包和ff19SB力场，分析了体系的均方根偏差（RMSD）、均方根涨落（RMSF）、回转半径（R_g）和氢键数量等指标。

此外，研究还利用GPS-MSP、NetPhos 3.1、GPS 6.0、RUBI、GPS-Uber和NetOGlyc4.0等工具预测了ACADM蛋白潜在的翻译后修饰（PTM）位点，包括甲基化、磷酸化、泛素化和糖基化。最后，通过GeneMANIA和STRING数据库构建了ACADM的基因-基因互作网络，以探索其参与的生物学通路和共表达关系。

³. 结果

^3.1. 有害nsSNP的识别与筛选

从数据库中共获得934个nsSNP，经过去重处理。利用七种生物信息学工具进行筛选，最终鉴定出16个被所有工具一致预测为有害的nsSNP。这些工具的具体预测结果如图3所示，其中SIFT预测302个有害，PolyPhen-2预测256个，Meta-SNP预测315个，PhD-SNP预测244个，PANTHER预测298个，SNAP预测287个，MutationAssessor预测64个为高风险。

^3.2. 突变对结构和功能的影响及蛋白稳定性评估

使用MutPred对上述16个nsSNP进行分析，其中12个（P62T, G99V, T193A, G195R, R206C, R206H, C244R, S245L, G267R, R281T, V373A, I375T）的预测分值大于0.7，表明它们有较高可能性改变蛋白质功能或结构，例如可能影响跨膜区域、催化位点或金属结合等。进一步的I-Mutant分析显示，这12个nsSNP均会导致ACADM蛋白稳定性下降。

^3.3. 进化保守性与结构建模

通过DeepREx-WS分析，发现G99、T193、G195、R206、R281、V373和I375这8个位点具有高度进化保守性，且多为埋藏残基，提示它们对蛋白质结构和功能至关重要。对其中8个最有害的nsSNP进行三维结构建模并与野生型比对，计算RMSD。结果显示，R206H（rs200724875）和R281T（rs121434282）两个突变具有最大的结构偏差，其RMSD值分别为1.26 ?和1.29 ?，超过了预设的1.25 ?阈值，表明它们引起了显著的结构变化。

^3.4. 分子对接分析

为探究突变对ACADM功能的关键环节——与ETF蛋白相互作用的影响，进行了分子对接。结果显示，野生型ACADM与ETF的结合能最低（-1341.4），且拥有最大的簇成员数（94），表明两者结合最稳定。相比之下，R206H-ETF和R281T-ETF复合物的结合能较高（分别为-1075.8和-1080.1），簇成员数更少（分别为61和73），提示结合亲和力下降。PDBsum的相互作用细节分析进一步证实，野生型复合物具有最多的氢键（24个）、非键接触（267个）和盐桥（5个），而两个突变体复合物的各类相互作用数量均减少，其中R281T突变体的相互作用最少。这表明R206H和R281T突变削弱了ACADM与ETF的结合，可能阻碍电子传递过程。

^3.5. 分子动力学模拟分析

200纳秒的MD模拟从动态角度验证了对接结果。RMSD分析表明，野生型ACADM-ETF复合物在整个模拟过程中最为稳定，RMSD值在3.8-4.2 ?之间小幅波动。而R281T-ETF复合物的RMSD在模拟前期急剧上升至17 ?，后期虽有下降但仍不稳定；R206H-ETF复合物的RMSD也升高至7 ?左右，表现出持续性不稳定。

RMSF分析反映了残基水平的柔性。野生型复合物的大部分残基涨落较低（<3 ?），而两个突变体复合物，特别是R281T，在多个区域显示出更高的残基涨落（最高达12-18 ?），表明局部结构更不稳定。

回转半径（R_g）分析揭示了复合物的结构紧密度。野生型复合物的R_g值最低（约37.2 ?），结构最紧凑。R281T和R206H突变体复合物的R_g值显著更高（分别约52.4 ?和47.4 ?），表明结构更为松散和扩展。

氢键分析显示，在模拟过程中，野生型复合物与ETF之间维持的氢键数量整体上多于两个突变体复合物，进一步证实突变削弱了蛋白质间的相互作用网络。

^3.6. 翻译后修饰位点预测

预测分析显示，ACADM蛋白可能存在多个磷酸化位点（如Ser、Thr、Tyr），但未发现明显的甲基化位点。此外，还预测了潜在的泛素化位点（主要位于Lys）和一个可能的糖基化位点（N22）。这些PTM位点可能参与调节ACADM的活性、定位或稳定性。

^3.7. 基因-基因互作网络

通过GeneMANIA和STRING分析发现，ACADM与多个基因存在密切关联。它可能与PPARA、ECHS1、ETFB等基因发生物理相互作用；与PPARA、ACADL、ACADS、ETFA、ETFB等基因存在共表达关系；并与ACADS、ACAD8、IVD等基因共享蛋白质结构域。这些互作网络提示ACADM被嵌入一个复杂的代谢调控体系中，其功能异常可能产生广泛影响。

⁴. 讨论

本研究通过系统的计算分析，揭示了ACADM基因中多个nsSNP的潜在有害性，并重点阐明了R206H（rs200724875）和R281T（rs121434282）两个突变导致MCADD的分子机制。这些突变位于进化上高度保守的残基上，可能破坏ACADM蛋白的稳定性。更重要的是，分子对接和MD模拟一致表明，这两个突变会显著削弱ACADM与其关键伙伴ETF蛋白的相互作用，表现为结合亲和力降低、复合物动态稳定性下降、结构更加松散以及界面氢键网络减少。这种相互作用的破坏很可能损害了脂肪酸β-氧化过程中的电子传递效率，从而导致MCADD的病理表现。其中，R281T（rs121434282）在以往的临床研究中已被报道与MCADD相关，本研究的计算结果为其致病性提供了结构生物学层面的解释。而关于R206H（rs200724875）的研究较少，本研究首次通过计算手段提示其潜在危害。研究还预测了可能影响ACADM功能的PTM位点，并构建了其基因互作网络，为了解该基因在更广泛代谢背景下的作用提供了线索。这些发现为理解MCADD的遗传基础、开发精准诊疗策略奠定了基础，但后续仍需通过湿实验（如酶活测定、细胞模型等）对这些计算预测进行验证。

⁵. 结论

综上所述，这项计算生物学研究成功鉴定出ACADM基因中一系列可能有害的nsSNP，并深入揭示了其中R206H和R281T两个突变通过破坏蛋白质稳定性及其与ETF的关键相互作用，进而可能导致MCADD的分子机制。这些发现强调了特定遗传变异在代谢疾病中的重要作用，为未来针对MCADD的机制研究、药物靶点发现以及个性化医疗策略的开发提供了重要的理论依据和方向。后续需要通过实验生物学手段进一步验证这些计算预测，并在动物模型中进行功能研究，以全面评估这些突变在疾病发生发展中的实际贡献。