中鏈酰基輔酶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的基因-基因互作網絡，以探索其參與的生物學通路和共表達關系。

從數據庫中共獲得934個nsSNP，經過去重處理。利用七種生物信息學工具進行篩選，最終鑒定出16個被所有工具一致預測為有害的nsSNP。這些工具的具體預測結果如圖3所示，其中SIFT預測302個有害，PolyPhen-2預測256個，Meta-SNP預測315個，PhD-SNP預測244個，PANTHER預測298個，SNAP預測287個，MutationAssessor預測64個為高風險。

使用MutPred對上述16個nsSNP進行分析，其中12個（P62T, G99V, T193A, G195R, R206C, R206H, C244R, S245L, G267R, R281T, V373A, I375T）的預測分值大于0.7，表明它們有較高可能性改變蛋白質功能或結構，例如可能影響跨膜區域、催化位點或金屬結合等。進一步的I-Mutant分析顯示，這12個nsSNP均會導致ACADM蛋白穩定性下降。

通過DeepREx-WS分析，發現G99、T193、G195、R206、R281、V373和I375這8個位點具有高度進化保守性，且多為埋藏殘基，提示它們對蛋白質結構和功能至關重要。對其中8個最有害的nsSNP進行三維結構建模并與野生型比對，計算RMSD。結果顯示，R206H（rs200724875）和R281T（rs121434282）兩個突變具有最大的結構偏差，其RMSD值分別為1.26 ?和1.29 ?，超過了預設的1.25 ?閾值，表明它們引起了顯著的結構變化。

為探究突變對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的結合，可能阻礙電子傳遞過程。

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之間維持的氫鍵數量整體上多于兩個突變體復合物，進一步證實突變削弱了蛋白質間的相互作用網絡。

預測分析顯示，ACADM蛋白可能存在多個磷酸化位點（如Ser、Thr、Tyr），但未發現明顯的甲基化位點。此外，還預測了潛在的泛素化位點（主要位于Lys）和一個可能的糖基化位點（N22）。這些PTM位點可能參與調節ACADM的活性、定位或穩定性。

通過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的機制研究、藥物靶點發現以及個性化醫療策略的開發提供了重要的理論依據和方向。后續需要通過實驗生物學手段進一步驗證這些計算預測，并在動物模型中進行功能研究，以全面評估這些突變在疾病發生發展中的實際貢獻。