[论文解读] Molecular multiresolution surfaces

Gang Wei, Yuhui Sun+2 more|ArXiv.org|Nov 1, 2005

Supramolecular Chemistry and Complexes参考文献 32被引用 26

一句话总结

本文提出分子多分辨率曲面作为统一的多尺度框架，用于建模生物分子曲面，通过在范德华表面进行连续溶剂的扩散映射，生成不同探针半径下的曲面。一种快速的局部谱演化核（LSEK）实现扩散方程的单步数值积分，准确捕捉溶剂排斥曲面、可及曲面和范德华曲面，且无奇点问题，已在环己烷和大型蛋白质上验证，与现有方法高度一致。

ABSTRACT

The surface of a molecule determines much of its chemical and physical property, and is of great interest and importance. In this Letter, we introduce the concept of molecular multiresolution surfaces as a new paradigm of multiscale biological modeling. Molecular multiresolution surfaces contain not only a family of molecular surfaces, corresponding to different probe radii, but also the solvent accessible surface and van der Waals surface as limiting cases. All the proposed surfaces are generated by a novel approach, the diffusion map of continuum solvent over the van der Waals surface of a molecule. A new local spectral evolution kernel is introduced for the numerical integration of the diffusion equation in a single time step.

研究动机与目标

• 开发一个统一的多尺度框架，用于在不同探针半径下建模生物分子曲面。
• 克服传统分子曲面计算中的拓扑与计算挑战，特别是奇点与自相交问题。
• 提供范德华曲面、溶剂可及曲面和溶剂排斥曲面的连续、几何一致描述。
• 通过新型数值积分方案，实现高效、单步计算多分辨率曲面。
• 在复杂生物分子上验证该方法，展示其准确性和可扩展性。

提出的方法

• 该方法通过带有曲率依赖扩散系数的反应-扩散方程，对范德华表面的溶剂扩散进行建模。
• 采用广义Fick定律，引入高阶项（q > 1）以考虑溶剂密度的非均匀性及通量相关性。
• 在溶剂可及曲面（SAR）体积上初始化溶剂密度（rho(r,0) = 0），其余区域设为rho0，以生成初始条件。
• 使用快速局部谱演化核（LSEK）在单个时间步长内求解扩散方程，复杂度为O(N)，空间精度可控。
• LSEK利用埃尔米特函数与加窗卷积计算密度演化，参数为Mh=88, Mα=32, σα=3.05hα。
• 多分辨率曲面作为演化后溶剂密度的等值面提取，范德华曲面为探针半径趋近零的极限。

实验结果

研究问题

• RQ1 能否通过单一连续框架统一表示范德华曲面、溶剂可及曲面和溶剂排斥曲面？
• RQ2 如何通过物理扩散过程捕捉分子曲面的多尺度特性？
• RQ3 能否通过快速单步数值方法替代迭代或复杂几何算法，实现分子曲面生成？
• RQ4 所提方法是否能避免传统曲面计算中固有的拓扑奇点与自相交问题？
• RQ5 该方法在保持高精度的同时，能否高效扩展至大型生物分子？

主要发现

• 该方法成功在不同探针半径下生成分子多分辨率曲面，范德华曲面为探针半径趋近零的极限。
• 在环己烷上，计算得到的溶剂排斥曲面在三个截面处与MSMS算法结果高度一致。
• 细胞分裂蛋白的中间分子曲面准确捕捉了关键结构对称性，同时抑制了非必要的细节。
• LSEK方法在固定精度下实现O(N)复杂度，实现扩散方程的高效单步积分。
• 该方法无需区分埋藏原子表面与非埋藏原子表面，避免了其他方法中常见的奇点问题。
• 数值验证确认了该方法在小分子和含数千个原子的大蛋白质上的准确性与鲁棒性。