分子动力学揭示精神分裂症相关SLC6A1变异对GAT-1局部结构网络的重塑

《International Journal of Molecular Sciences》：Molecular Dynamics Reveals Structural Consequences of Schizophrenia-Associated SLC6A1 Variants

　　

SLC6A1基因编码GAT-1转运体，该转运体介导中枢神经系统（central nervous system, CNS）中γ-氨基丁酸（γ-aminobutyric acid, GABA）的再摄取。该基因的错义变异与神经精神疾病（包括精神分裂症）的发生相关。尽管体外实验表明转运活性丧失，但这些实验并未阐明特定氨基酸替换如何破坏蛋白质内部的局部相互作用。为描绘这些突变的结构后果，研究人员利用微秒级分子动力学（molecular dynamics, MD）模拟分析了人GAT-1的内向开放apo态（PDB: 7Y7V）。本研究将野生型转运体与突变形式进行比较，包括精神分裂症相关变异R211C和A305T，以及通过基于深度学习的致病性预测因子选出的A346S变异。结果表明，这些突变并未在整体上使蛋白质结构失稳；相反，它们诱导了局部接触网络内的特异性重排。R211C替换消除了细胞外环2（extracellular loop 2, EL2）与跨膜螺旋4（transmembrane helix 4, TM4）界面处的一个氢键。A305T改变了钠离子结合位点附近的TM1–TM6界面，而A346S增强了TM1/TM6–EL4相互作用。所识别的结构位移提示了致病突变如何重塑GAT-1相互作用网络并改变其构象动力学。这些发现为后续旨在药物重定位的分子对接以及完整GAT-1转运周期的全面MD建模奠定了基础。

**论文解读：精神分裂症相关SLC6A1变异的分子动力学结构后果**

精神分裂症是一种异质性的神经发育障碍，其GABA能假说认为皮质抑制不足导致兴奋-抑制平衡失调。SLC6A1基因编码GAT-1转运体，负责突触间隙中γ-氨基丁酸（γ-aminobutyric acid, GABA）的高效再摄取。SLC6A1错义变异已与癫痫、智力障碍、自闭症谱系障碍及精神分裂症等神经精神疾病相关联。然而，体外功能实验虽然能报告转运活性丧失或表面表达异常，却无法解释特定氨基酸替换如何在原子水平上破坏蛋白质内部的局部相互作用网络。GAT-1属于LeuT折叠家族，其功能依赖于跨膜螺旋的协同运动以及底物和离子的偶联，因此需要原子级分辨率的方法来阐明变异的结构后果。为此，研究人员以人GAT-1内向开放apo态结构（PDB: 7Y7V）为起点，利用微秒级分子动力学（molecular dynamics, MD）模拟，系统比较了野生型与三个错义变异体——精神分裂症相关的R211C和A305T，以及基于AlphaMissense致病性预测被选为可能良性的A346S——以揭示变异如何重塑局部接触网络并改变构象动力学，为后续药物重定位和完整转运周期模拟奠定基础。

在方法上，研究人员通过PyMOL进行计算机定点诱变，构建R211C、A305T和A346S突变模型；使用CHARMM-GUI在混合脂质双层中组装蛋白-膜系统，并由GROMACS完成每个体系四次独立的1微秒生产模拟。后续分析采用了几何判据定义堆积接触、氢键和盐桥，使用Cohen’s d与置换检验校准阈值，以硬过滤（|d|≥3且|Δocc|≥0.15）筛选可靠变化；HOLE2计算孔隙半径；基于广义互信息（generalized mutual information, gMI）构建每个重复的动态残基网络，分析通信路径和中心性变化。变异来源包括精神分裂症患者中de novo的R211C、SCHEMA数据集中的A305T以及AlphaMissense低风险评分的A346S。

研究结果分述如下：

**2.1 变异位点图谱（Variant Sites on the GAT-1 Map）**：R211位于细胞外环2（extracellular loop 2, EL2）紧邻跨膜螺旋4（transmembrane helix 4, TM4），A305位于TM6非螺旋片段（TM6a–TM6b断裂处）邻近底物位点，A346位于EL4。三者所处结构环境不同，提示它们可能通过不同机制影响转运体。

**2.2 孔道内衬残基选择（Selection of Pore-Lining Residues）**：通过HOLE2和3.5 ?距离判据定义渗透路径的孔道内衬残基，固定了30个缢缩残基，并筛选出133对用于堆积接触分析的配对，这些配对主要集中于TM1–TM6界面，包括位置305，但未覆盖R211或A346附近，因此堆积层无法直接报告这两个位点的接触变化。

**2.3 软效应量列表与置换偶然基线（Soft Effect-Size Lists Against the Permutation-Derived Chance Floor）**：在n=4的重复设计中，软阈值|d|≥0.8会产生约30%的偶然超过率；因此采用硬阈值|d|≥3且|Δocc|≥0.15。初步计数显示软过滤下的选择数接近随机期望，证明了硬过滤的必要性。

**2.4 TM1/TM6堆积接触（Packing Contacts at TM1/TM6）**：硬过滤后，WT–A305T仅保留Gly59–Leu306一个堆积接触；A305T–A346S对比保留Gly297–Ser302、Val58–Ala61和Gly59–Leu306三个接触。WT–R211C无堆积命中。A305T和A346S使同一局部基序发生相反方向的机械位移（r²=0.026，p=0.0038）。时间过程显示，Gly59–Leu306在A305T的四个重复中三个发生不可逆断裂，而野生型中基本维持，说明该接触丧失源于局部侧链重排而非整体螺旋束重排。

**2.5 氢键标记替换位点局部效应（Hydrogen Bonds Mark Substitutions Locally）**：硬过滤后共鉴定出21个氢键命中，显著高于置换偶然基线（P=4.6×10^?5），归结为12个独特残基对，其中三对包含替换侧链：Met197–Arg211（EL2）、Gly112–Thr305（TM6胞质面）、Ser346–Arg351（EL4）。

**2.6 盐桥：硬命中数为零（Salt Bridges: Hard Count Zero）**：34个候选盐桥中无一对通过硬阈值。虽然K278–E283表现出高|d|但Δocc极小；R211C可能的Arg211–Asp171盐桥在野生型中已占用为零（最小距离约6.6 ?），因此该变异不存在盐桥网络塌陷。

**2.7 细胞外门在该时间尺度上保持关闭（The Extracellular Gate Remains Closed on This Timescale）**：所有16条轨迹中，时间中位缢缩半径仅为0.34–0.74 ?，小于脱水钠离子的0.95 ?，且突变体与野生型无显著差异。因此观察到的TM1/TM6堆积改变和氢键变化均发生在稳定内向开放、细胞外门关闭的构象中。

**2.8 每重复动态网络揭示配体结合区的变异特异性重组（Per-Replicate Dynamic Networks Reveal Variant-Specific Reorganization of Ligand-Binding Regions）**：gMI网络显示整体网络差异处于重复间变化范围内，但中心性分析发现Tyr140、Phe294、Ser302在各变异中反复失去介数中心性，Asp451、Ser456失去节点强度。路径分析表明：R211C使从211位到S1底物结合区和噻加宾识别区的通信路径普遍变长，并削弱Leu218–Tyr140中继；A305T使到Tyr60/Ala61等路径变短，而到Tyr140/Ser396路径变长，形成不对称网络；A346S使到大多数S1靶点路径变短，形成Ser346–Phe339–Phe70–Phe294/Ser295通路。Gly297–Ser456强耦合边在所有三个突变体中均减弱，提示底物结合TM6段与细胞外门TM10侧的动态协调受损。

讨论部分指出，双阈值过滤有效降低了偶然选择；堆积接触和氢键信号均超过置换基线，而盐桥未显示可靠变化。三种变异呈现不同结构表型：R211C主要引起EL2局部极性重排（Met197–Arg211丢失），而非静电网络崩溃；A305T和A346S对TM1/TM6同一基序产生相反位移，其中A305T通过破坏Gly59–Leu306桥并产生不对称S1网络而干扰底物/离子偶联，与实验观察到的严重功能丧失一致；A346S则增强EL4-TM1/TM6耦合，可能导致内向开放状态过度稳定。结合功能数据和AlphaMissense评分，A305T严重性最高，R211C居中，A346S可能最轻。但该研究仅基于内向开放apo态，且缺少底物和离子，因此结论为机制性预测，需后续多状态模拟和实验验证。

结论翻译如下：分子动力学模拟表明，所研究的错义变异并未导致GAT-1整体失稳，内向开放构象和关闭的细胞外门仍然保持，三个替换均未破坏转运体的整体折叠。由于结构扰动是局部且变异特异性的，研究人员提出三种候选机制，这些机制仍是基于接触占有率推断的预测模型，有待实验验证。基于MD轨迹分析，研究人员预测：A305T可能重排底物-离子结合位点附近的TM1–TM6界面，破坏Gly59–Leu306接触并在中央口袋内产生不对称的相互作用网络；R211C可能间接改变蛋白功能，因为EL2中Met197–Arg211相互作用的丧失削弱了其与底物结合口袋和细胞外门的耦合；A346S可能引起最轻微的结构偏差，保留关键的Gly59–Leu306桥，但通过增强EL4与TM1/TM6的相互作用限制了结构可动性，可能过度稳定内向开放态。由于目前分析主要局限于单一结构状态，尚无法确定这些结构破坏是否在转运循环的其他阶段持续存在或如何影响状态转换。未来需系统研究GAT-1的多个功能状态：检验GABA和噻加宾结合下的野生型与突变体；模拟外向开放、闭合和内向开放之间的构象转变；以及在GABA、Na⁺和Cl^?存在下评估变异体。最终，将结构建模与遗传和功能数据整合，有助于从仅按序列位置分类变异转向识别离散的分子病理机制，为精准精神病学中的变异解读和药理策略提供结构基础。