周启涛生物技术工作室
结构与蛋白质组2025/08/20· 编译自 Nature

IF50.5|电子衍射实验测定原子部分电荷的新方法

iSFAC建模让每个原子的部分电荷可直接从晶体结构解析中获得,无需特殊软件或高分辨率数据。

iSFAC建模让每个原子的部分电荷可直接从晶体结构解析中获得,无需特殊软件或高分辨率数据。

这篇按生信解读号的读法拆:先看背景和问题,再看作者怎么做、做出了哪几条硬发现,最后把数据和代码收走。

研究背景

原子部分电荷是理解分子结构、相互作用和反应性的核心概念,但缺乏精确的量子力学定义,实验测定一直困难。传统X射线晶体学对电子细节不敏感,而电子衍射因电子与静电势直接作用,理论上可提供电荷信息。此前虽有汇聚束电子衍射等方法用于无机材料,但通用性有限。本文提出iSFAC建模,将离子散射因子与中性散射因子线性叠加,为每个原子引入一个额外参数,在标准电子晶体学流程中即可精修得到绝对标度的部分电荷。

科学问题

卡在哪里: 此前没有通用的实验方法能定量测定化合物中每个原子的部分电荷。X射线衍射对电子分布不敏感,而电子衍射虽有潜力,但现有方法仅适用于少数无机材料或需要极高分辨率数据,无法广泛应用于有机分子和生物分子。

本文要回答: 本文旨在建立一种基于三维电子衍射的通用方法,在常规晶体结构解析中同时获得所有原子的部分电荷,并验证其可靠性和适用范围。

技术路线

作者主要用了:三维电子衍射、iSFAC建模、SHELXL精修、量子化学计算

作者在标准电子晶体学流程中引入iSFAC建模:对每个非氢原子,其散射因子表示为中性原子和离子散射因子的线性组合,比例系数作为额外参数精修;氢原子则直接用H+散射因子。通过SHELXL的FVAR机制实现,无需专门软件。该方法在环丙沙星、酪氨酸、组氨酸和ZSM-5沸石上验证,并与多种量子化学方法计算的电荷进行Pearson相关分析。

Fig. 1 · 原文图,按文末许可署名使用
Fig. 1 · 原文图,按文末许可署名使用

图 1 Fig. 1。图1展示了iSFAC建模得到的原子部分电荷示意图,正负号表示电荷符号,直观呈现分子中电荷分布的不均匀性,为后续具体化合物分析奠定概念基础。

核心亮点

亮点 1|iSFAC建模实现原子级部分电荷测定

iSFAC建模为每个原子引入一个额外参数,表示离子散射因子的比例,从而得到绝对标度的部分电荷。该方法在标准电子衍射数据上即可运行,无需特殊软件或超高分辨率。对环丙沙星、酪氨酸、组氨酸和ZSM-5的测试显示,有机分子的实验电荷与量子化学计算(如ADCH、CM5、NPA)的Pearson相关系数达0.8以上,验证了方法的可靠性。

Fig. 2 · 原文图,按文末许可署名使用
Fig. 2 · 原文图,按文末许可署名使用

图 2 Fig. 2。图2对比了传统精修(每个元素类型一个散射因子)与iSFAC建模(每个原子独立散射因子)的差异。iSFAC通过线性组合中性原子和离子散射因子,为每个原子精修出独特的散射曲线,从而提取部分电荷。

亮点 2|氢原子处理是iSFAC建模的关键

传统晶体学中氢原子位置和热参数常被约束。iSFAC建模要求氢原子使用H+散射因子并自由精修,否则非氢原子的部分电荷变得无意义。对组氨酸的分析显示,包含氢原子后,实验电荷与计算电荷的Pearson相关系数从低于0.2提升至0.8以上。这一处理还揭示了氢键中质子位置的细微差异,如组氨酸中N1–H1⋯O2氢键的两种构型,其能量差异达59 kcal mol−1。

Fig. 3 · 原文图,按文末许可署名使用
Fig. 3 · 原文图,按文末许可署名使用

图 3 Fig. 3。图3给出了环丙沙星分子的实验部分电荷值,标注了标准不确定度。图中可见羧基碳C18带正电荷(+0.11e),而大多数非氢原子带负电荷,氢原子带正电荷,符合化学直觉。

亮点 3|部分电荷反映化学环境敏感性

组氨酸的两种多晶型(空间群P212121和P21)具有相似的分子结构但不同的晶体环境,其部分电荷存在可检测的差异。例如,羧酸根氧原子O2的电荷在不同数据集中从−0.24e到+0.10e变化,与氢键角度相关。这表明iSFAC建模不仅能测定分子内禀电荷,还能捕捉分子间相互作用对电荷分布的影响,为理解晶体场效应提供实验依据。

Fig. 4 · 原文图,按文末许可署名使用
Fig. 4 · 原文图,按文末许可署名使用

图 4 Fig. 4。图4展示了酪氨酸和组氨酸的实验部分电荷,按元素类型着色。两种氨基酸的羧酸根碳均带负电荷(酪氨酸C9为−0.19e,组氨酸C6为−0.25e),反映了羧酸根中电子的离域效应。

亮点 4|iSFAC建模对X射线数据不适用

作者测试了iSFAC建模对高质量X射线数据的适用性,包括Ca2+酒石酸盐和已发表的组氨酸结构。结果显示,X射线数据下精修不收敛,部分电荷超过1e甚至出现负电荷氢原子,均不合理。这归因于X射线对电子分布不敏感,而电子衍射直接探测静电势,因此iSFAC建模是电子衍射特有的优势。

Fig. 5 · 原文图,按文末许可署名使用
Fig. 5 · 原文图,按文末许可署名使用

图 5 Fig. 5。图5显示了组氨酸二聚体中N1–H1⋯O2氢键的四种构型及其DFT计算能量。DS5构型(氢键角175°)比DS1(155°)稳定59 kcal mol−1,而传统riding atom模型能量更高,说明iSFAC建模能揭示氢键几何与电荷的关联。

亮点 5|方法稳健性及分辨率依赖性

通过分辨率截断测试,发现当分辨率优于1 Å时部分电荷稳定,低于此阈值偏差增大,与Sheldrick规则一致。对氢原子施加温和约束可恢复可靠性。此外,Rcomplete交叉验证表明,单个反射的缺失对部分电荷的影响与常规精修中原子位移参数的影响相当,说明方法稳健。

生信可带走

这一块是给做分析的人用的,不是科普点缀。

  • 方法栈: 三维电子衍射、iSFAC建模、SHELXL精修、量子化学计算
  • 公开数据: 原文未给出公开组学登录号
  • 代码: 原文未给出公开 GitHub/GitLab 仓库
  • 谁该点开原文: 从事电子晶体学、计算化学或分子模拟的研究人员,以及对原子电荷实验测定感兴趣的化学和材料科学工作者。

带走一句

iSFAC建模为晶体学分析增加了电荷维度,可在常规电子衍射数据上获得原子级部分电荷,为分子模拟和反应机理研究提供实验基准。

本页供学习交流,不替代原文,也不代表原作者、期刊或资助方立场。引用科学结论请以正式发表版本为准。 图表按 CC BY(可转载,需署名) 使用。 版权声明

微信二维码

需要同类分析?扫码加微信,一对一沟通需求,交付后终身售后。

添加微信

需要同类分析?扫码加微信,一对一沟通需求,交付后终身售后。