基于不完美二维晶体的三维重建：一种混合实空间/傅里叶空间方法

一、 研究团队、发表信息与学术背景

本研究由W.O. Saxton*、W. Baumeister以及M. Hahn合作完成。W.O. Saxton和W. Baumeister来自德国马丁斯里德的马克斯-普朗克生物化学研究所（Max-Planck-Institut für Biochemie），而M. Hahn则来自德国杜塞尔多夫大学（Universität Düsseldorf）生物物理与电子显微镜研究所。该研究成果以论文形式发表于1984年的《Ultramicroscopy》期刊第13卷，第57-70页。

该研究属于结构生物学与电子显微学交叉领域，特别是针对生物大分子二维晶体的三维结构解析。长期以来，利用电子显微镜（EM）和傅里叶变换方法解析二维晶体结构（如著名的紫膜结构）取得了巨大成功。然而，一个普遍存在的严重限制是，天然或合成的二维生物分子晶体极少能达到高度完美的周期性有序。晶体畸变（crystal distortions）会严重干扰基于传统傅里叶方法的分析，因为这些方法高度依赖于完美的晶体周期性。为了克服这一瓶颈，研究者们发展并应用了相关平均法（correlation averaging），该方法通过实空间图像处理来校正晶格畸变，从而从畸变的晶体图像中提取出平均化的、高质量的单元胞投影。本研究的目标，正是将这种实空间处理技术与传统的基于傅里叶空间的晶体学三维重建方法相结合，形成一种“混合实空间/傅里叶空间”工作流程，以实现对不完美二维晶体的可靠三维结构解析。文章以细菌细胞包膜的一种二维晶体（HPI层）为例，详细阐述了该方法的完整流程、技术挑战及初步应用结果。

二、 详细工作流程

本研究的工作流程是一个系统性工程，其核心思想是利用相关平均法处理每个倾斜角度的投影图像，获得平均化的单元胞，然后从这些平均单元胞中提取傅里叶变换数据（晶格线数据），最终通过三维逆傅里叶变换重建出密度图。图1的框图清晰地概括了这一过程，具体步骤如下：

1. 数据收集与样品制备： 研究使用负染（negative staining）技术制备样品，具体采用含0.05%葡萄糖的1%磷钨酸钠（PTA）溶液。这种组合在提供良好衬度的同时，能相对均匀地染色整个载网，并对HPI层结构提供了令人满意的保存（分辨率优于2纳米）。图像采集采用连续倾转系列（continuous tilt series）的方式，使用允许近乎无限倾转的样品杆，最高倾角可达约80°。这种策略便于后续精确确定实际倾转角，并确保所有投影图像具有共同的起源。为了最小化辐射损伤，采用了低剂量技术（总剂量约22000 e/nm²）。图像选择基于两个标准：正确的聚焦和消像散，以及样品被染色剂良好包埋。

2. 投影图像的平均化（相关平均法）： 对每个倾转角度下采集的原始投影图像，使用相关平均法进行处理。该方法的核心是：从同一张图像中选取一个小的“参考”区域，通过计算该参考区域与图像其他部分的互相关函数（cross-correlation function），定位所有单元胞的位置。然后，通过质心法（centre-of-mass criterion）精修这些位置，并迭代此过程（使用初次得到的平均图像作为新参考），最终将所有单元胞对齐并平均，得到一个代表该倾角下结构投影的高信噪比平均图像。本研究使用了标准化的、高度自动化的程序来完成此步骤。对于高倾角图像，由于焦深很浅，平均处理被限制在平行于倾转轴的狭窄条带区域内。

3. 倾转轴与倾转角度的确定： 显微镜测角仪的读数提供了名义倾转角，但需要精确确定样品平面相对于倾转轴的倾角（α）、倾转轴在图像中的方位角（δ）以及可能存在的角度偏移（φ₀）。本研究提出了一种基于投影图像中晶格基向量（base vectors）变化的几何方法。由于投影会导致图像在垂直于倾转轴的方向发生cosφ倍的收缩，通过拟合一系列不同倾角下观测到的晶格基向量变化，可以反推出上述所有参数。该方法利用包含高、低倾角的整个图像系列，显著提高了确定精度（例如，本研究中倾转轴方位角的标准误差为0.34°）。研究发现，对于高倾角（>60°）图像，仅依靠相邻图像的相关性或质心定位来对齐会引入累积误差，因此需要更全局的方法。

4. 投影图像的相对对齐（确定共同原点）： 在将不同倾角的投影数据合并到三维傅里叶空间之前，必须确定它们之间的相对位移，即建立一个“共同原点”。文章探讨了多种方法： * 质心法：对于HPI层这种大部分质量集中在中心“核心”、通过较轻“辐条”连接的结构，可以手动勾画出单个分子的边界，计算其质心。由于三维结构的质心必然投影到每个二维投影的质心上，因此可以通过对齐质心来实现投影对齐。 * 互相关法：计算相邻倾角投影图像之间的互相关函数，通过峰值位置确定相对位移。但该方法在倾角间隔大时可靠性下降，且可能产生误差累积。 * 基于晶格线数据的自举法：这是本研究采用并推荐的更复杂方法。其基本思路是：首先利用初步对齐的数据（例如，将所有图像与0°图像互相关对齐的结果）构建一个初始的、不完整的晶格线数据集。当要加入一张新的投影图像时，程序会根据已有的晶格线数据，通过沿晶格线短线段平均，生成一个对应于该新投影角度的估计傅里叶截面。然后将这张新投影图像与该估计截面进行互相关，找到最佳对齐位置后再将其数据加入晶格线数据集。这种方法利用了所有已积累数据的对称性信息，能进行大角度间隔图像间的有效比较，防止误差累积，并对垂直于倾转轴方向的对齐也提供了更可靠的基础。

5. 晶格线数据的提取： 二维晶体的三维傅里叶变换被限制在一组垂直于晶体平面的离散“晶格线”上。传统方法是从原始投影图像的傅里叶变换中提取这些线上的采样点。本研究则从相关平均法得到的平均化单元胞中提取等效数据。具体操作是：首先，将每个平均单元胞按照其投影后的晶格进行重新采样（re-sampled），使其正好被整数个采样点覆盖。然后，对这个重新采样的单元胞进行离散傅里叶变换（DFT），其变换值在特定频率位置上的采样，就对应于相应晶格线在特定高度（z*）上的一个数据点。z*的值由样品的空间几何取向（通过前述的旋转矩阵R_ij计算）决定。由于不同倾角下采样点在晶格线上的位置是不均匀的，且数据存在噪声，研究开发了专门的软件程序（在“SEMPER”系统上）来管理、显示和手工拟合这些晶格线数据。用户可以在屏幕上查看每条晶格线的模量和相位（或实部和虚部）数据点，并手动绘制平滑曲线穿过这些点，然后从曲线上均匀重采样，从而得到平滑、均匀采样的晶格线数据，用于后续重建。

6. 三维密度图的恢复与后处理： 从平滑插值后的晶格线数据恢复三维密度分布，是通过三维逆傅里叶变换实现的。本研究采用了逐线变换再组合的方法：首先对每条晶格线进行一维逆傅里叶变换，然后将所有线的对应采样点组合起来，构建出每一个垂直于晶体法线方向的“切片”的二维傅里叶变换，最后对每个切片进行二维逆傅里叶变换，得到实空间的密度切片。这些切片最初位于一个与晶体晶格对应的采样格点上，需要通过重采样恢复到常规的笛卡尔坐标网格。最终得到的三维密度体数据可以以灰度图或等高线图的形式逐层显示。

三、 主要结果

1. 完整流程的建立与验证：研究成功建立了一套从数据采集、图像处理到三维重建的完整混合方法工作流程。通过对细菌HPI层二维晶体的实际应用，证明了该流程的可行性。重建得到的三维密度图在0°和60°方向的计算投影，与原始图像的平均投影吻合良好，验证了流程的自洽性。
2. 技术难题的解决方案：
   • 倾转角优化：研究从采样理论出发，推导出为了在目标分辨率r下对厚度为t的样品进行均匀采样，最优的倾角间隔应与cosφ成正比，而非均匀间隔。这比之前提出的按正切均匀间隔的方案更高效。文章提供了针对不同t/r比值的优化倾角表。
   • 高精度倾转参数确定：通过拟合整个倾转系列中晶格基向量的变化，能够以高精度（<0.5°）确定倾转轴方位角、样品平面倾角以及测角仪偏移，即使对于未校准的高倾角（>80°）也适用。
   • 投影对齐的优化：比较了多种对齐方法，发现基于晶格线数据的“自举”互相关法能产生最平滑的晶格线数据，是解决高倾角投影对齐难题的有效方案。研究还测试了不同加权的互相关函数（如强制所有傅里叶分量模量为1），发现它们能产生更尖锐的互相关峰，但峰值位置与标准互相关基本一致。
   • 晶格线数据处理：开发了专门的软件工具用于晶格线数据的可视化、手工曲线拟合和重采样。研究发现，用实部和虚部形式展示数据比用模量和相位形式更便于手动拟合曲线，因为后者在过零点时会出现尖点和间断，容易引入误差。
3. HPI层的三维结构模型：应用该流程，研究人员获得了HPI层的三维密度图，并据此构建了一个物理模型。模型显示了一个具有“粗糙”内表面（在完整细胞壁中朝向细菌外膜）的蛋白结构。然而，重建结果也揭示了当前方法的一些局限性。
4. 重建结果的局限性与挑战：
   • 密度分布不理想：理论上，负染样品应呈现蛋白（高密度）和染色剂（低密度）的二元分布。但实际重建的密度图显示，两者之间的过渡相对平滑，且密度值在晶体内部随远离中心切片而逐渐下降，导致结构看起来比实际更薄。文章分析认为，这主要不是由于高倾角数据缺失（“缺失锥”问题）造成的，因为本研究包含了超过60°的倾角数据，显著扩展了低阶晶格线的覆盖范围。根本原因在于从噪声数据中拟合出的晶格线数据质量有限，导致重建出的结构在垂直方向被一个高斯型扩展函数所模糊。
   • 阈值选择困难：由于缺失（0,0）晶格线（零频分量）的数据，每个密度切片的平均密度被强制设为零。因此，无法确定一个全局阈值来区分蛋白和染色剂，也无法直接计算出正确的分子体积。研究者不得不对每个切片主观地选择一个阈值，主要依据是等高线的密度梯度以及蛋白区域的连通性和切片间的连续性。文章对通过晶体褶皱图像获取（0,0）线数据的常用方法提出了质疑，因为HPI层的褶皱图像显示染色水平变化很大，且褶皱处的分子严重变形，无法提供纯净的侧面视图。

四、 结论与意义

本研究的主要结论是，通过结合相关平均法（实空间）和晶格线重建法（傅里叶空间），发展出了一种能够从结晶度不完美的二维晶体中获取可靠三维结构的方法。这项工作的重要意义在于，它极大地扩展了电子晶体学的应用范围，使得大量因晶格畸变而无法用传统傅里叶方法分析的生物样品（如许多天然存在的二维晶体）有望进行高分辨率三维结构解析。

研究的科学价值体现在提供了一套详细、可操作的技术方案，解决了从数据采集策略（倾转角优化、连续倾转系列）、图像处理（相关平均、高精度倾转参数确定、投影对齐）、到数据整合与重建（晶格线提取、手工拟合、三维逆变换）等一系列关键步骤中的实际问题。该方法不依赖于完美的晶体周期性，因此更具普适性。

五、 研究亮点

1. 方法学的创新与整合：首次系统性地阐述并实现了将“相关平均”这一实空间图像处理技术与传统的傅里叶空间三维重建流程相结合，形成完整的“混合”方法，专门针对不完美二维晶体。
2. 解决关键技术难题：
   • 从理论推导出最优倾角采样方案（cosφ比例），提高了数据采集效率。
   • 提出了基于整个倾转系列晶格基向量拟合的高精度倾转参数后验确定方法，精度高且通用性强。
   • 发展了基于已积累晶格线数据的“自举”投影对齐法，有效解决了高倾角图像对齐和误差累积问题。
   • 开发了用于晶格线数据可视化与交互式拟合的专用软件工具，提高了数据处理的灵活性和可靠性。
3. 流程的完整性与验证：提供了从样品制备到最终模型构建的详尽工作流程，并通过实际生物样品的重建和投影比对，验证了整个流程的自洽性和可行性。
4. 对局限性的坦诚分析与前瞻：不仅展示了成功，还深入分析了当前方法重建结果中密度分布不理想、阈值选择主观等局限性，将其归因于晶格线数据质量而非数据范围，并提出了可能的改进方向，如通过平均更大图像区域、进行传递函数修正、整合多个数据系列、或引入迭代精修、最大熵约束等更先进的算法来利用先验知识。

六、 其他有价值内容

文章在附录中提供了描述样品旋转的旋转矩阵的详细推导，为其他研究者实现相关计算提供了便利。此外，文章还讨论了高倾角成像中由于样品厚度增加导致的电子束衰减问题，指出在满足一定条件（厚度 t ≤ r²/4λ，r为分辨率，λ为电子波长）下，仍可近似用线性投影理论处理，且衰减效应可通过后续的归一化步骤进行补偿。最后，作者展望了该方法的潜力，指出未来可以结合更复杂的选择性平均程序，并利用阴影技术获得的表面轮廓信息等，来进一步优化和丰富三维结构模型。