PINK1突变（T313M和L347P）的动态结构分析揭示了线粒体自噬介导的常染色体隐性帕金森病（ARPD）中一个关键的分子扰动，即磷酸化丝氨酸65泛素识别点

《Frontiers in Molecular Neuroscience》：Dynamic structural profiling of PINK1 mutations (T313M and L347P) reveals a vital molecular perturbation namely phospho-Serine 65 ubiquitin recognition point in mitophagy mediated autosomal recessive Parkinson’s disease (ARPD)

　　

常染色体隐性帕金森病（ARPD）源于PINK1/Parkin通路功能障碍导致的线粒体自噬受损，其中PINK1介导的泛素在Ser65位点的磷酸化对于通路激活至关重要。然而，实验局限性掩盖了疾病相关突变对内在无序区域（IDRs）和翻译后修饰（PTM）动态的影响。方法：采用整合计算流程筛选825个PINK1错义变异，在激酶结构域中鉴定出两个高置信度有害突变T313M和L347P。变异优先级通过保守残基、IDR和PTM分析进行补充，随后进行蛋白质-蛋白质对接、分子动力学模拟、MM/PBSA结合自由能计算、主成分分析（PCA）和动态交叉相关矩阵（DCCM）分析。结果：T313M与一个保守磷酸化位点重叠，而L347P映射到保守活性位点残基，IDR分析提供了补充支持。对接分析显示，从野生型（–88.4 ± 8.2）到L347P（–81.9 ± 3.4）和T313M（–77.4 ± 4.9），结合亲和力逐渐降低，伴随静电稳定化（–357.6 → –260.6 → –235.5 kcal/mol）和埋藏表面积（1741.6 → 1624.1 → 1547.4 ?²）减少。分子动力学模拟表明，T313M产生了最大的结构和动态扰动，而L347P诱导了中等程度的去稳定化并增加溶剂暴露。尽管MM/PBSA分析表明所有系统的结合能量学大致相当，但PCA和DCCM分析显示突变复合物中构象灵活性增加和残基通讯改变，尤其是T313M。讨论：这些发现建立了突变诱导的结构动力学与Ser65处PINK1-泛素识别受损之间的机制联系，为理解ARPD早期线粒体自噬损伤提供了突变特异性框架，并支持未来的分子评估和靶向治疗开发。

帕金森病（Parkinson’s disease, PD）是全球第二大神经退行性疾病，其特征为黑质致密部多巴胺能神经元进行性退化，导致运动症状如运动迟缓、僵硬和静止性震颤。约8%的PD患者为遗传性，其中常染色体隐性帕金森病（autosomal recessive Parkinson’s disease, ARPD）主要由PINK1（PTEN-induced putative kinase 1）或Parkin基因功能丧失突变引起。PINK1是一种丝氨酸/苏氨酸激酶，定位于线粒体外膜，在受损线粒体上积累后，通过磷酸化泛素（ubiquitin, Ub）的Ser65位点以及磷酸化Parkin的泛素样结构域（ubiquitin-like domain, UBL），激活Parkin，启动线粒体自噬（mitophagy）。泛素在Ser65的磷酸化（生成phospho-Ser65-Ub）是线粒体自噬启动的关键事件，但实验方法难以解析疾病相关突变对内在无序区域（intrinsically disordered regions, IDRs）和翻译后修饰（post-translational modifications, PTMs）动态的影响。传统结构生物学技术受限于这些区域的灵活性和瞬态构象，无法精确表征突变如何扰动PINK1–Ub相互作用。因此，研究人员采用整合计算框架，模拟PINK1–Ub复合物在激活状态下的结构，系统筛选ARPD相关错义突变，研究其结构扰动机制。

研究人员通过整合计算流程，从公共数据库（Ensembl、ClinVar等）收集825个PINK1错义变异，多工具筛选后鉴定出T313M和L347P两个高置信度有害突变，位于激酶结构域。随后进行IDR和PTM分析、基于直系同源的序列保守性分析、活性位点预测、蛋白质-蛋白质对接（HADDOCK）、分子动力学模拟（GROMACS，200 ns×3重复）、MM/PBSA结合自由能计算、主成分分析（principal component analysis, PCA）和动态交叉相关矩阵（dynamic cross-correlation matrix, DCCM）分析，旨在揭示突变对PINK1–Ub复合物结构稳定性、动态行为和结合能量的影响。研究结果发表于《Frontiers in Molecular Neuroscience》。

主要关键技术方法包括：①多工具致病性筛选（SIFT、PolyPhen-2、EVE、SNP&GO、Missense3D、CUPSAT），从Ensembl和ClinVar数据库获取825个变异；②结构建模与序列分析（AlphaFold、PRALINE、PrankWeb），预测PINK1激酶结构域激活态（含pSer228、ATP和Mg²⁺）及Ub的pSer65状态；③蛋白质-蛋白质对接（HADDOCK），利用预测活性残基定义AIRs；④分子动力学模拟（GROMACS，CHARMM36力场，200 ns三重复）；⑤能量学分析（MM/PBSA，gmx_MMPBSA）；⑥构象动态分析（PCA、DCCM）。未使用患者样本，所有数据来源于公开数据库。

**3.1 模型优化与系统准备**：通过SWISS-MODEL优化AlphaFold预测的PINK1结构，Ramachandran图显示86.4%残基在最有利区域，结构可靠性良好。Ub结构（PDB:1UBQ）单体化，并将PINK1 Ser228和Ub Ser65磷酸化，加入ATP和Mg²⁺模拟催化活性环境。

**3.2 基于直系同源的保守性和催化口袋结构表征**：多序列比对（人、小鼠、赤拟谷盗）显示激酶结构域高度保守，HRDLKPEN催化环和DFG基序严格不变。ATP和Mg²⁺位于典型结合口袋，Asp384（DFG基序）与Mg²⁺配位（~2.0–2.3 ?），Lys364和Asp366稳定ATP磷酸基团。

**3.3 多工具致病性和稳定性评估**：825个变异经SIFT、PolyPhen-2、EVE、SNP&GO、Missense3D、CUPSAT筛选，仅T313M和L347P在所有分析中均表现出有害性，且位于保守激酶结构域。L347P热力学去稳定化最强（ΔΔG = -7.23 kcal/mol），T313M次之（ΔΔG = -2.79 kcal/mol）。

**3.4 内在无序区域与翻译后修饰组织结构**：IDR分析揭示两个主要无序区域（IDR-1: 1–136，IDR-2: 175–211），其中IDR-2具有高结合倾向（ANCHOR峰值0.95）。PTM分析识别16个位点，T313与磷酸化相关（NetPhos评分0.553），L347P位于保守催化基序附近。两突变均位于有序激酶结构域，但靠近IDR-2和PTM位点，可能通过远距离调控影响功能。

**3.5 进化保守性与活性位点特征**：PRALINE比对显示T313和L347在五个哺乳动物物种中完全保守，且位于PrankWeb预测的活性位点（Rank 1，评分9.28）内。基于此，将PINK1活性残基（T313、A321、A324、A327、L343、L347、T348、A352）与Ub活性残基（I44、L8、K49、S65、K68、R72、R74）用于对接。

**3.6 PINK1–泛素复合物的结合能量学**：HADDOCK对接显示，野生型复合物结合最强（HADDOCK评分-88.4 ± 8.2），L347P次之（-81.9 ± 3.4），T313M最弱（-77.4 ± 4.9）。突变导致静电稳定化显著降低（WT: -357.6 kcal/mol，L347P: -260.6，T313M: -235.5）和埋藏表面积减少（1741.6 → 1624.1 → 1547.4 ?²），表明磷酸化Ub识别受损。

**3.7 突变诱导的结构重排：柔性、紧凑性和溶剂暴露**：200 ns分子动力学模拟（三重复）显示，T313M的RMSD（1.0700 ± 0.2185 nm）和RMSF（1.0189 ± 0.4278 nm）最高，结构去稳定化最强；L347P RMSF（0.8416 ± 0.0593 nm）中等，但SASA最大（360.2930 ± 3.7158 nm²），表明溶剂暴露增加。Rg分析显示T313M松散（3.1760 ± 0.0733 nm），L347P与WT接近。氢键数量波动大，反映界面接触动态重排。

**3.7.1 MM/PBSA结合自由能分解**：MM/PBSA计算显示WT（-1551.84 ± 559.87 kcal/mol）、T313M（-1547.44 ± 556.01）和L347P（-1546.54 ± 558.16）结合自由能大致相当，标准偏差大，不支持突变显著改变结合亲和力。突变效应主要体现于结构动态而非热力学稳定性。

**3.8 构象景观分析：PCA衍生的自由能景观和动态交叉相关矩阵**：PCA显示WT占据紧凑低能盆地，T313M构象空间最分散，多亚稳态，L347P介于两者之间。DCCM显示WT具有协调的残基运动网络，T313M破坏相关/反相关模式，长程通讯受损，L347P扰动程度较轻。这些结果与RMSD、RMSF、Rg分析一致，表明T313M对动态架构破坏最大。

讨论部分总结：研究表明，T313M和L347P通过不同机制破坏PINK1–Ub相互作用。T313M位于保守磷酸化位点，削弱静电稳定化和构象一致性，导致整体结构和动态扰动；L347P位于保守活性位点，主要引起局部结构去稳定化和溶剂暴露增加。PCA和DCCM揭示突变增加构象灵活性并扰乱残基通讯，尤其是T313M。这些发现建立了突变诱导的结构动力学与Ser65磷酸化Ub识别受损之间的机制联系，为理解ARPD早期线粒体自噬损伤提供了突变特异性框架，支持未来分子评估和靶向治疗开发。研究完全基于计算分析，未进行实验验证，但为早期遗传筛查和机制研究提供了假设生成性见解。

研究结论部分翻译：本研究采用整合计算框架评估人PINK1突变在泛素结合界面上的潜在结构和功能影响。结构评估确认了建模PINK1的可接受立体化学特性。多步骤SNP排序始终将T313M和L347P鉴定为最有害变异，并经序列保守性、结构干扰和热力学不稳定性预测支持。IDR和PTM分析显示，尽管突变位于结构化的激酶结构域，但其邻近动态调控区域可能削弱构象柔性、磷酸化可能性和相互作用动力学。保守性分析和活性位点预测证实T313和L347占据进化保守的功能区域，并邻近泛素结合界面。L347P位于保守激酶结构域，可能直接损害结构完整性和催化残基排列；T313M对应预测的磷酸化位点，主要影响翻译后修饰调控。分子对接和动力学模拟表明，突变（尤其T313M）降低结合亲和力，削弱静电稳定化，并增加构象灵活性和残基通讯扰动。MM/PBSA显示结合自由能大致相当，但动态分析表明突变导致界面稳定性下降和构象异质性增加。这些发现共同表明，T313M和L347P诱导细微但功能显著的结构扰动，损害PINK1–Ub在Ser65处识别，为ARPD早期线粒体自噬损伤提供突变特异性机制框架，并支持未来分子评估和靶向治疗开发。