实验数据与密度泛函理论计算的对应关系阐明了变铅青链霉菌铜储存蛋白对亚铜离子的装载与卸载机制

《Inorganica Chimica Acta》：Correspondence between experimental data and density functional theory calculations elucidates mechanism of cuprous ion loading to and unloading from the copper storage protein from Streptomyces lividans

　　

亚铜（Cu(I)）在细菌细胞质中的储存通过Csp3（铜储存蛋白）蛋白家族的作用得以实现。理解Csp3成员用于安全结合并防止细胞内Cu(I)潜在毒性效应的配位化学，对于理解它们在细菌铜稳态中的更广泛作用非常重要。本研究报道了针对变铅青链霉菌Csp3（SlCsp3）金属结合位点模型的密度泛函理论（DFT）计算和从头算分子动力学研究。第一组计算针对四核Cu₄模型（这些模型被认为主导Cu(I)初始装载进入20个可用配位位点），阐明了配位几何受配体选择——硫醇盐（半胱氨酸）与咪唑（组氨酸）——控制的程度。相应模型的能量最小值与通过X射线晶体学在SlCsp3初始Cu(I)装载时实验确定的Cu₄簇良好匹配。发现咪唑配体强烈不利于四面体结构，且这种趋势在硫醇盐/咪唑杂配配合物中更为明显。这可能是决定SlCsp3外核（入口点）选择混合Cys/His Cu(I)配位的因素之一。对所有晶体学观察到的Cu(I)位点的进一步DFT计算支持两步序列的合理性，即外核先被Cu(I)离子填充，同时一个羧酸残基（Asp61）将外核与内核分隔，对进一步Cu(I)装载施加热力学屏障。最后但同样重要的是，对SlCsp3的脱辅基、Cu(I)装载和部分装载形式进行的分子动力学计算确认了预期的动态特征，尤其是在脱辅基和部分装载形式中，四螺旋束内的储存通道不对称地打开以暂时容纳Cu(I)离子的进入。

**研究背景与目的**

铜是生物体必需的过渡金属，但细胞内游离的亚铜离子（Cu(I)）具有潜在毒性，因此细菌进化出精密的铜稳态网络。铜储存蛋白Csp3（copper storage protein）家族可在细菌细胞质中安全储存大量Cu(I)，例如来自变铅青链霉菌（*Streptomyces lividans*）的SlCsp3，其四螺旋束单体可结合多达20个Cu(I)离子。前期实验通过X射线晶体学、定点突变和停流动力学已揭示SlCsp3的Cu(I)装载呈两步过程：Cu(I)首先进入由Cys/His组成的外核（位点15–20），形成杂配四核Cu(I)簇，随后跨越Asp61残基进入内核（位点1–14）。然而，控制这一过程的配位化学原理、高亲和力位点间的快速Cu(I)传递机制，以及蛋白质构象动态对装载的贡献仍不清楚。为此，研究人员开展了密度泛函理论（DFT）计算、从头算分子动力学（AIMD）和经典分子动力学（MD）模拟，旨在从电子结构和构象动态层面阐明SlCsp3的Cu(I)装载与卸载机制。相关成果发表在《Inorganica Chimica Acta》。

**主要技术方法**

研究采用DFT PBE0/def2-TZVP方法优化四核Cu(I)簇模型，配体以甲基硫醇盐（SMe?）和咪唑（Im）分别模拟半胱氨酸和组氨酸，考察同配与杂配体系的几何偏好。基于X射线晶体结构（PDB 6EK9）构建含20个Cu(I)位点的通道模型，通过单点能和冻结坐标优化计算各单Cu(I)位点的相对能量。采用AIMD（ADMP模型）验证局部配位稳定性。经典MD使用YASARA/AMBER14力场，对脱辅基（PDB 6EI0）、部分装载（PDB 6Q58）和完全装载（PDB 6EK9）的SlCsp3单体进行150 ns溶剂化模拟，分析骨架RMSD、残基波动、溶剂可及表面积等。

**研究结果**

*3.1 同配硫醇盐和咪唑配体的四核Cu(I)簇的内在几何偏好*

对[Cu₄(SMe)₄]⁰和[Cu₄(Im)₄]⁴⁺模型的DFT优化表明，两者的全局能量最低结构均为菱形（rhombic）排列。硫醇盐模型中，四面体Cu₄核心的能量比全局最低值高17.2 kcal/mol；咪唑模型中，四面体结构相对更不稳定，开放型3+1结构则更稳定。这说明咪唑配体强烈不利于紧凑的四面体几何，其原因与硫的软碱性质及18电子规则有关。

*3.2 杂配硫醇盐和咪唑配体的四核Cu(I)簇的内在几何偏好*

对[Cu₄(Im)₂(SMe)₂]²⁺和[Cu₄(Im)(SMe)₄]模型的计算显示，杂配体系的势能面与同配体系显著不同：菱形结构不再占据全局最低，四面体结构无法进入低能区域，而开放型3+1结构和近似线形结构获得额外稳定。该结果与SlCsp3外核初始形成的杂配Cys/His四核Cu(I)簇的晶体学观察一致。

*3.3 Cu(I)结合通道的热力学性质*

基于完全装载SlCsp3晶体结构（6EK9）的20位点模型，DFT单点能计算揭示了各Cu(I)位点的相对能量景观。全局最小值位于内核中心附近的位点11和13；外核位点18、19能量较低；而由Asp61羧酸根配位的位点14和15呈能量高点，构成外核通往内核的热力学屏障。该结果从理论上解释了为何Cu(I)先填满外核，再突破Asp61屏障进入内核，与实验提出的两步装载机制吻合。

*3.4 经典分子动力学研究Cu(I)结合通道的构象动态*

对脱辅基（6EI0）、部分装载（6Q58）和完全装载（6EK9）的150 ns MD模拟显示，三种状态下四螺旋束均保持折叠，但骨架具有明显动态性。脱辅基蛋白中，Asp61羧酸根溶剂暴露且活动性强；部分装载和完全装载后，该残基完全埋藏并与Cu(I)稳定配位，SASA降为零。半胱氨酸硫醇盐的溶剂可及性随Cu(I)装载逐步下降，完全装载后几乎不可及。此外，脱辅基形式中四螺旋束的铜装载端会不对称地张开，形成临时通道，而装载形式中该运动明显减弱。

**讨论与结论**

讨论部分指出，Cu₄簇在铜蛋白中并不常见，但类似结构见于N₂O还原酶的CuZ位点和铜抗性蛋白CopG，表明四核铜簇可能具有多种功能。SlCsp3外核的杂配Cu₄簇与DFT预测的开放菱形结构一致，组氨酸配体通过阻止过紧密的四面体几何来避免阻塞后续Cu(I)进入。Asp61作为热力学屏障，其在不同Csp3家族成员中可被Asn或Met替代，这解释了不同Csp3装载动力学差异。MD结果进一步支持Asp61作为负载驱动的“门”残基，且脱辅基蛋白的通道开放特性为Cu(I)进入提供了结构基础。

研究结论可总结为：第一，DFT计算表明Cu₄簇倾向于开放而非闭合的四面体结构，组氨酸配体强化了这一倾向，这为SlCsp3入口处混合Cys/His配位的选择提供了化学依据。第二，位点能量图识别出Asp61羧酸根是分隔外核与内核的热力学屏障，支持外核先装载的两步机制。第三，经典MD揭示了Asp61门控的动态特性，并证明脱辅基状态下四螺旋束的装载端在150 ns时间尺度内可逆地开放，为Cu(I)进入提供通道。这些结果综合展现了理论计算与实验数据之间的高度对应，深化了对细菌铜储存蛋白工作机制的理解。