Cu2+与Zn2+离子非键相互作用对β-淀粉样蛋白1–16及25–35片段聚集的影响──分子动力学模拟研究

《ACS Chemical Neuroscience》：The Effect of Cu2+ and Zn2+ Ions’ Nonbonded Interactions on the Aggregation of β-Amyloid 1–16 and 25–35 Fragments─A Molecular Dynamics Simulation Study

　　

阿尔茨海默病与细胞外β-淀粉样蛋白（β?Amyloid，Aβ）聚集体的形成和积累有关，其毒性主要归因于可溶性寡聚物种，正如寡聚体假说所提出。与此同时，金属离子假说认为，过渡金属离子如Cu2+和Zn2+可直接调节Aβ片段的聚集过程及其结构稳定性。在本研究中，研究人员采用分子动力学（Molecular Dynamics，MD）模拟产生集体变量，首次在水溶液中研究这些离子及其浓度对β?淀粉样蛋白1–16和25–35片段聚集的影响。聚集的自由能分布表明，Cu2+离子的存在略微降低β?淀粉样蛋白1–16片段聚集相关的能量，提示其作为调节因子部分稳定寡聚物。相比之下，研究人员首次报道Zn2+离子在水溶液中并不降低β?淀粉样蛋白1–16片段聚集的能垒。Zn2+对Aβ1–16酸性残基表现出更高亲和力，与Cu2+相比建立更频繁但选择性较低的接触。主要发现是Cu2+和Zn2+离子调节β?淀粉样蛋白的早期聚集路径。特别是在高浓度下，Cu2+有利于形成富含反平行β?折叠、结构有序的小型原寡聚物，而非单纯增加聚集体尺寸，而在所研究条件下Zn2+未表现类似效应。这种金属诱导的富含β?折叠的低阶寡聚物稳定化（尤其在Cu2+中显著）将早期二级结构转变确定为淀粉样毒性的主要决定因素及阿尔茨海默病的关键分子事件。本研究提供的证据增强了对阿尔茨海默病寡聚体假说和金属离子假说的理解。

研究背景方面，阿尔茨海默病是一种进行性神经退行性疾病，主要损害记忆、认知功能和行为，最终导致失能和丧失自主性，其病理特征包括海马体和大脑皮层萎缩、细胞外β?淀粉样蛋白（Aβ）肽老年斑沉积以及由过度磷酸化tau蛋白组成的细胞内神经原纤维缠结。Aβ肽来源于淀粉样前体蛋白（Amyloid Precursor Protein，APP）经β?分泌酶和γ?分泌酶依次切割产生，常见变体为Aβ₄₀和Aβ₄₂，后者更易聚集并形成毒性沉积物，引起突触功能障碍、神经元死亡和脑功能进行性下降。目前获批治疗仅能缓解症状和延缓进展，无法治愈该病。为解释其机制，存在寡聚体假说和金属离子假说：寡聚体假说认为可溶性Aβ寡聚物引起快速突触损伤；金属离子假说指出内源性金属离子如Cu²⁺、Fe³⁺、Zn²⁺可调节Aβ寡聚物毒性，促进形成神经毒性寡聚物种并加剧突触功能障碍。Aβ_1–16包含N端第6、13、14位组氨酸（Histidine，His），对Cu²⁺和Zn²⁺配位至关重要，且与Cu²⁺相互作用可产生活性氧导致氧化应激；Aβ_25–35虽小但保持聚集、形成β?折叠及诱导与全长肽相似的神经毒性能力，适合计算模拟聚集及寡聚物初始形成。经典分子动力学（Molecular Dynamics，MD）模拟中过渡金属离子与肽片段相互作用的处理本质上是近似的，无法明确捕获部分共价性、可变配位几何和电子极化效应，本研究目的并非在量子力学水平模拟离子配位化学，而是探究非键金属–肽相互作用如何在相关时间和长度尺度上调节Aβ片段聚集行为，使用经典力场（离子表示为嵌入水模型中的带电粒子）通过采样集体变量系统探索聚集路径及相关能量分布，构成工作核心创新点，研究人员明确此方法提供金属–肽相互作用的简化描述而非详细配位说明，结论限于定性趋势，并讨论经典力场中离子参数化局限及对溶剂模型的依赖。

研究人员开展的研究是通过分子动力学模拟考察Cu²⁺和Zn²⁺离子及其浓度对Aβ_1–16和Aβ_25–35片段在水溶液中聚集的影响，得出Cu²⁺轻微降低Aβ_1–16聚集相关能量、部分稳定寡聚物，Zn²⁺不降低其聚集能垒但对酸性残基亲和力更高且接触更频繁但选择性较低，高浓度Cu²⁺促进小而结构有序、富含反平行β?折叠的原寡聚物形成而Zn²⁺无类似效应，金属离子调节早期聚集路径且Cu²⁺诱导的β?折叠富集低阶寡聚物稳定化是阿尔茨海默病关键分子事件，该结果增强对寡聚体假说和金属离子假说的理解，重要意义在于从原子水平阐释金属离子如何定性重定向早期聚集路径及二级结构转变，为金属介导神经毒性提供分子基础并指向早期二级结构转变作为潜在治疗干预靶点，论文发表在《ACS Chemical Neuroscience》。

关键技术方法方面，研究人员使用GROMACS 2024软件开展MD模拟，采用TIP4P水模型，体系置于10×10×10 nm立方盒中并由约32000个TIP4P水分子溶剂化，添加0.15 mol L^–1盐浓度（Na⁺/Cl^–）及中和净电荷的反离子，先进行最大50000步或力达100 kJ mol^–1nm^–1的最陡下降能量最小化（范德华与静电截断半径1.2 nm，周期性边界条件），随后进行5 ns的NVT（粒子数、体积、温度恒定）系综平衡至310 K（leapfrog积分器，时间步2 fs，velocity rescaling恒温器τ=0.1 ps，LINCS约束重原子与氢键），再进行NPT（粒子数、压力、温度恒定）系综平衡（压力1 bar，τ=2.0 ps，指数弛豫压力耦合与随机细胞重新标度气压计c?rescale），生产步在同样温压条件下运行500 ns，长程与范德华作用用Particle Mesh Ewald法（实空间截断1.2 nm，傅里叶网格间距0.16 nm，四阶插值），非键作用含排斥、色散和库仑项，排斥与色散用Lennard?Jones (6–12)势，静电用库仑势；聚集分析采用GROMACS Colvar模块定义集体变量（Collective Variables，CVs），以0.6 nm截断代表残基间显著接触，每1000步采样，CV基于两片段原子对接近度计算，距离趋近截断时贡献增加；用Python脚本将CV值组织为45个区间直方图得到概率分布P(ξ)并归一化，按F(ξ)=?k_BT ln P(ξ)（k_B为玻尔兹曼常数，T=310 K）计算自由能分布并转为kJ mol^–1；径向分布函数（Radial Distribution Functions，RDF，g(r)）用GROMACS RDF模块计算金属离子与Aβ特定残基侧链原子在最后50 ns轨迹中的相互作用；Aβ_1–16取自PDB ID 2BEG的Aβ₄₂全长并用PyMOL处理，评估10个片段无金属、10/20/40个Cu²⁺、10个Zn²⁺及20个片段无金属与20个Cu²⁺体系，Cu²⁺用OPLS?AA/M力场，Zn²⁺用AMBER14SB_OL24力场并额外用同力场模拟Cu²⁺以比较；Aβ_25–35取自PDB ID 1QYT并用PyMOL处理，评估10/20/40个片段无金属及10个片段加10个Cu²⁺体系，无金属用OPLS?AA力场，含Cu²⁺用OPLS?AA/M力场。

研究结果部分保留原文小标题并依内容简述如下：

System Containing the Aβ_1–16Fragment：研究人员从PDB ID 2BEG获取Aβ₄₂中Aβ_1–16并用PyMOL制备，设置多条件体系（10片段无金属、10片段+10/20/40 Cu²⁺、10片段+10 Zn²⁺、20片段无金属、20片段+20 Cu²⁺），选用OPLS?AA/M（Cu²⁺）与AMBER14SB_OL24（Zn²⁺及对照Cu²⁺）力场。通过分析得出不同金属种类、浓度及片段数对聚集影响的体系基础设置及力场选择依据。

System Containing the Aβ_25–35Fragment：研究人员从PDB ID 1QYT获取Aβ_25–35并用PyMOL处理，设置10/20/40片段无金属及10片段+10 Cu²⁺体系，无金属用OPLS?AA、含Cu²⁺用OPLS?AA/M力场，指出所用片段与金属浓度高于生理值以提升采样并在可及模拟时间内观察聚集，因生理浓度下淀粉样聚集缓慢而无法在典型原子模拟时长内观察到。

Results and Discussions：研究人员通过MD模拟CV得到的聚集行为分析接触概率与自由能，结合RDF揭示金属离子与残基相互作用。对10个Aβ_1–16片段，Cu²⁺使P(ξ)重新分布，10个Cu²⁺时某些聚集构型采样频率更高，20和40个Cu²⁺因静电屏蔽效应减少高度紧凑构型群体；20个片段时有无Cu²⁺的P(ξ)在大部分ξ范围更相似，低浓度下ξ对金属调节更敏感而高密度下被肽–肽接触主导。Zn²⁺使P(ξ)移至不同ξ区域，几乎可忽略地微幅促进罕见ξ访问。Aβ_25–35有无Cu²⁺的P(ξ)相对群体不同但形状与范围相似，因该片段仅有Met35弱配位且水合主导使Cu²⁺影响有限。自由能分布F(ξ)显示所有体系低ξ对应低自由能、高ξ（3000–5000）对应低概率与能量增加；10个Aβ_1–16加10个Cu²⁺时低能态ξ区间改变，3000–5000 ξ区低能态数可能增加（瞬态状态），20和40个Cu²⁺时该区能垒升高；20个片段时含Cu²⁺系统能量增加延迟，部分稳定聚集体。Zn²⁺系统在高ξ区F(ξ)上升范围和斜率改变，不降低罕见ξ能垒且减少其热力学可达性，因观察中至多6个Zn²⁺直接参与聚集而Cu²⁺全参与且无分散离子；换用AMBER力场比较Cu²⁺与Zn²⁺仍保持Cu²⁺更易达高ξ聚集态、Zn²⁺高ξ能量陡峭上升趋势（高ξ为罕见事件需谨慎解释），RDF模式在不同力场一致。Aβ_25–35有无Cu²⁺的自由能曲线相似，因缺乏典型金属结合残基。RDF显示Cu²⁺与Aβ_1–16的Asp1、Asp7在～0.2 nm有明锐峰，Glu3、Glu11类似，Lys16末端羧基～0.25 nm有峰；His6、His13、His14在0.5–0.7 nm呈多峰宽分布，反映动态多样空间排列；Zn²⁺与Asp1、Asp7峰～0.4 nm较弱，His6/13/14呈0.3和2.0 nm双峰，Glu3/11～0.4 nm中等强度峰，显示Zn²⁺配置更多样且亲和力稍低；Aβ_25–35中Cu²⁺与Met35末端羧基在0.5 nm内有显著峰，反映频繁接近但非特异稳定配位，Asn27分布弥散。模拟体系演化为原寡聚组装（小型动态亚稳态聚集体，早期瞬态中间体），500 ns内未达全局平衡而采样亚稳态集合（连续重排与3–4片段小簇及约10片段大聚集体互变）；F(ξ)作为比较性有效热力学描述，用于评估金属如何调节紧凑寡聚构型可达性及重塑聚集态相对群体，不断言绝对定量或动力学机制。

Implications：研究人员指出聚集呈非线性且高度依赖浓度，金属–肽化学计量与局部配位环境关键重塑聚集景观；低浓度Aβ_1–16与Cu²⁺形成大型原寡聚组装（约9–10肽，伴相近数Cu²⁺），结构无序或微α?螺旋，无β?折叠；Cu²⁺增至20和40时原寡聚群体共存不同尺寸簇，仍无序且无明显α?螺旋；加倍肽与Cu²⁺浓度时系统常居2–4片段小型原寡聚物种并伴有较大组装，重定向至低阶原寡聚体并伴随Aβ_1–16凝聚（反平行β?折叠）出现，Cu²⁺作为构象催化剂选择性稳定浓缩Aβ_1–16排列，小型富含β?折叠聚集体最神经毒性且与尺寸无关而与二级结构和次序有关。Zn²⁺在低浓度也促大型原寡聚组装但参与金属离子较少，作用较弱或少价，调节稳定性与结构动态不同。Aβ_25–35与Cu²⁺作用较弱但自身易采取反平行β?折叠，说明序列本身编码β?折叠倾向。原寡聚形成、二级结构获得及Cu²⁺/Zn²⁺相互作用紧密耦合，Cu²⁺促进小型β?折叠丰富原寡聚体为金属诱导神经毒性提供分子基础，早期金属介导构象转变是关键治疗靶点；无金属Aβ_25–35体系中低ξ最常被采样，高ξ概率连续下降，浓度调节聚集态相对群体，无金属时罕见ξ仍难访问且能量不利。

Limitations：研究人员指出常规MD表征罕见事件困难，因其低频且被高自由能垒分隔，单一轨迹有限采样无法稳健估计平衡概率或区分随机波动，对初始条件敏感且独立模拟可能得不同路径，系统常受限在亚稳盆地内屏障跨越不足而无法可靠估计自由能差、转变率或态寿命，时间尺度不匹配（MD为ns–μs，罕见事件可能μs–ms以上）使观察可能为偶然采样或瞬态波动；因此研究中罕见事件不作为机制结论基础，而强调基于平衡轨迹部分的统计稳健可观测量（概率分布与沿CV的自由能分布），罕见构型仅作系统可访问结构极端与能量成本的定性示意，结论基于集合行为（时间平均与分布分析）而非单一转变叙述，承认严格表征需多独立轨迹或增强采样方法。

结论部分翻译总结：本研究结果与大量证明金属离子调节Aβ_1–16和Aβ_25–35片段聚集的实验证据一致。研究人员发现Cu²⁺和Zn²⁺对构象稳定性和动态行为施加不同影响，取决于特定片段和环境；RDF分析揭示Cu²⁺和Zn²⁺与带负电残基及组氨酸反复接触，在相对频率和构象上存在差异，未观察到持久配位位点，因为采用的经典力场描述倾向于非键相互作用且未明确解释稳定金属–配体配位复合物所需的共价特征；在此框架下，Cu²⁺可能表现出更复杂的调节作用，影响聚集动态但不一定促进直接纤维形成，而Zn²⁺可能主要表现为间接效应，保持大部分溶剂化并降低罕见高度聚集态的热力学可达性。尽管取得进展，完全阐明Cu²⁺和Zn²⁺离子与各种Aβ寡聚态的相互作用仍是挑战，此处观察提供初步见解以支持未来研究整合实验与理论方法进一步阐明金属离子在Aβ聚集中的作用。总体而言，最重要发现是Cu²⁺和Zn²⁺离子不仅以定量方式增强Aβ聚集，而是定性重定向早期聚集路径：提高Cu²⁺和肽浓度有利于出现小型原寡聚物而非大型组装，同时诱导反平行β?折叠形成（神经毒性淀粉样物种的定义性结构标志）。与广泛实验证据一致（证明小型Aβ寡聚物致病潜力），结果不直接涉及毒性而是提供早期聚集事件的结构见解；特别是模拟提示金属离子（尤其Cu²⁺）可调节原寡聚物内二级结构，促进独立于聚集体尺寸的浓缩构象，这些观察提供分子水平视角，有助于解释实验观察到的金属结合、早期结构组织与毒性之间的关系。在未来研究中，也将处理这些金属–肽相互作用的共价组分，以获得金属介导β?淀粉样蛋白聚集早期阶段的更全面图像。