利用分子建模研究错义突变对SYNGAP1双重Ras/Rap激活机制的影响

《Journal of Structural Biology: X》：Investigating the effects of missense mutations on the dual Ras/rap activation mechanism of SYNGAP1 using molecular modelling

　　

中文摘要翻译： SynGAP采用经典的RasGAP和Plexin样RapGAP机制。通过计算机模拟方法揭示了SynGAP–GTPase选择性的决定因素。体细胞癌症和SynGAP相关疾病（SRD）错义突变显示出显著的重叠。错义变异可影响SynGAP–Ras/Rap复合物的稳定性和催化活性。SynGAP对GTPase活性位点的预组织解释了其对Rap的活性比Ras高100倍的原因。

论文解读

一、研究背景与科学问题

SynGAP（突触GTPase激活蛋白）由SYNGAP1基因编码，在前脑神经元兴奋性突触的突触后致密区（PSD）中高丰度表达。SynGAP通过稳定GTPase的Switch I和Switch II区域使其进入催化活性构象，并向活性位点提供关键的“精氨酸手指”（Arg-finger）以稳定GTP水解的过渡态，从而加速Ras/Rap-GTPase的内在GTP水解，使其失活。SynGAP参与调控细胞生长、分化和细胞骨架动力学等信号通路，不仅在大脑中发挥作用，在甲状腺、支气管等其他人体组织中也有表达。

SynGAP蛋白约140 kDa，远大于其GAP结构域（约21 kDa），包含N端的PH（Pleckstrin同源）、C2和GAP结构域，以及C端的卷曲螺旋（CC）结构域。这种多结构域组织使其具有催化功能之外的多种功能：PH和C2结构域可能介导膜结合，CC结构域驱动三聚化。此外，SynGAP通过其α1亚型C端的PDZ结合基序，在长时程增强（LTP）过程中限制PSD95的PDZ结合位点对AMPAR–TARP复合物的可用性，从而调节突触可塑性。

目前存在的主要科学问题包括：SynGAP相关的错义突变诊断困难。SynGAP相关疾病（SRD）是一系列严重程度不同的神经发育疾病，包括智力障碍（ID）、癫痫和自闭症样行为，由SYNGAP1基因的新生突变引起。尽管SRD是一种罕见病（估计发病率约为每10万人6例），但它是儿童ID的主要遗传原因之一，约占所有病例的1%。截短变异容易被归类为致病性，但错义突变——可能占SRD病例的大多数——由于其微妙且往往难以捉摸的结构效应，诊断仍然具有挑战性。此外，虽然GAP功能丧失突变可能不是SRD的根本原因，但已有报道称它们会改变内在神经元兴奋性，因此在某种程度上可能是SRD相关症状和疾病表型的促成因素。

另一个关键问题是双重特异性GAP的非经典Rap激活机制尚不明确。在GTP水解过程中，RasGAP使用精氨酸手指稳定GTP的离去γ-磷酸，而Ras-Gln61定向亲核水分子。相比之下，Rap-GTPase依赖催化性“天冬酰胺拇指”（Asn-thumb）协调亲核水分子。双重特异性GAP（包括SynGAP）共享RasGAP折叠和精氨酸手指，但缺乏RapGAP的天冬酰胺拇指和同源性，因此其非经典Rap激活机制一直不清楚。虽然目前尚无实验解析的双重特异性GAP–Rap复合物结构，但来自Plexin–Rap X射线晶体结构的见解提示了一种“Plexin样”非经典机制，其中Rap-Gln63通过稳定开放的Switch II构象而被定位以进行催化。这种开放构象由Plexin特异的氢键网络稳定，而该网络在SynGAP等双重特异性GAP中缺失。相反，SynGAP被认为依赖其C2结构域来稳定Switch II，这与C2结构域对SynGAP和GAP1家族成员使Rap失活不可或缺的观察结果一致。

本研究旨在通过原子级分子动力学（MD）模拟、计算机模拟诱变和结合自由能计算，阐明SynGAP对双重Ras/Rap-GTPase激活的分子机制，以及体细胞癌症和SRD相关错义突变如何调节SynGAP–Ras/Rap的相互作用。

二、主要技术方法

研究人员主要采用了以下关键技术方法：基于AlphaFold数据库获取SynGAP-α2亚型的PH-C2-GAP区段结构（残基198–730），以人GAP–Ras复合物晶体结构（PDB: 1WQ1）为模板通过骨架叠合建模SynGAP–Ras复合物，以Plexin–Rap晶体结构（PDB: 4M8N）为模板建模SynGAP–Rap复合物；对野生型和突变型复合物进行原子级分子动力学（MD）模拟；使用FoldX-MD方法（结合FoldX能量计算与MD轨迹）和Rosetta Flex ddG方法进行计算机模拟诱变和结合自由能计算；通过MM/GBSA分析选择最低能量构象作为模板。变异数据来源于ClinVar数据库（N=322）和癌症数据库（N=199，包括COSMIC、Cancer Gene Census和dbSNP等）。

三、研究结果

3.1 SynGAP对Ras的激活遵循经典RasGAP机制

通过结构建模和MD模拟，研究人员证实SynGAP–Ras复合物中，SynGAP的精氨酸手指直接与GTP的γ-磷酸相互作用，Ras-Gln61采用催化取向并与催化水分子结合。这一催化排列与经典RasGAP–Ras复合物（如HRasGAP–Ras，PDB: 1WQ1）一致，表明SynGAP对Ras的失活遵循涉及精氨酸手指和Ras-Gln61的经典RasGAP机制。

3.2 SynGAP对Rap的激活遵循Plexin样非经典机制

研究人员建立了首个SynGAP–Rap1b复合物的原子模型。在最终复合物模型中，SynGAP采用Plexin样排列，通过其GAP和C2结构域稳定Rap表面。关键相互作用包括Arg596和Ser604与Rap-Switch II的接触，以及SynGAP-Arg401与Rap-Asp69之间形成的盐桥。序列比对显示，SynGAP的α6螺旋基序中，Pro600和Ser604对应于Plexin的天冬酰胺残基。初始复合物模型显示SynGAP激活片段与Rap-Switch I之间存在空间冲突，而采用开放环构象后解决了这一冲突，实现了兼容的结合几何。在这一机制中，SynGAP稳定延伸的Rap-Switch II构象，将Rap-Gln63定位为类似于Ras-Gln61的作用，通过涉及C2（Arg401）和GAP（Arg596、Ser604）结构域的协调相互作用实现催化。这解释了SynGAP对Rap的失活效率比对Ras高约100倍的现象——GTPase活性位点的预组织是这一差异活性的基础。

3.3 错义突变对SynGAP–Ras/Rap复合物稳定性和催化活性的影响

研究人员分析了214个错义变异，主要来自癌症数据库，其中33个被预测会破坏复合物稳定性和/或GAP介导的GTP水解几何结构，包括8个ClinVar列出的可能与SRD相关的变异。通过FoldX-MD和Rosetta两种方法的结合自由能计算，识别出多个显著影响复合物稳定性的变异。例如，R272W、R401Q、R401W、R485P、K492E、R596H、K628N等变异在两个复合物中均表现出显著的结合亲和力降低（ΔΔG_bind ≥ 1.0 kcal/mol）；R485P和R485L-F484A组合变异导致显著的折叠稳定性降低（ΔΔG_fold ≥ 3 kcal/mol）。值得注意的是，R401Q和R401W在SynGAP–Rap复合物中引起显著的去稳定效应（ΔΔG_bind分别为3.9±0.6和3.0±0.6 kcal/mol），而在SynGAP–Ras复合物中也有明显影响（1.6±0.9和1.1±0.7 kcal/mol），这与Arg401在Rap结合中形成盐桥的关键作用一致。相反，K442R和K628R变异在两个复合物中均表现出稳定化效应（负ΔΔG值）。ClinVar注释的VUS变异（R401Q、E486K、K492E、S604L、N635S）和可能良性变异（R479Q、R596H、R621Q）中，多个被预测会破坏复合物稳定性，提示这些变异可能具有潜在的致病性影响。

四、讨论与结论

该研究通过分子建模和MD模拟，首次建立了SynGAP–Rap1b复合物的原子模型，揭示了SynGAP双重Ras/Rap特异性的结构基础。研究证实SynGAP对Ras的失活遵循经典RasGAP机制，而对Rap的失活遵循Plexin样非经典机制，其中C2结构域在稳定Rap-Switch II构象中发挥关键作用。SynGAP对GTPase活性位点的预组织解释了其对Rap的活性比Ras高约100倍的现象。

研究还发现，体细胞癌症和SRD相关错义突变存在显著重叠（60个变异在ClinVar和癌症数据库中均有报道），提示SYNGAP1的GAP功能失调可能在癌症发生和神经发育疾病中具有共同的病理机制。33个变异被预测会破坏复合物稳定性和/或GTP水解催化几何结构，其中多个ClinVar列出的VUS变异表现出显著的破坏效应，为这些变异的致病性评估提供了结构生物学依据。

该研究为理解SynGAP的双重GAP活性的分子机制提供了重要的结构基础，为SRD相关错义突变的致病性评估提供了计算预测工具，也为针对Ras/Rap信号通路异常的癌症治疗策略提供了潜在靶点。