[论文解读] Continual Model of Medium II: Universal adaptive algorithm for triangulation of smooth molecular surface

Oleg Kupervasser, N. E. Wanner|arXiv (Cornell University)|Jul 8, 2011

Various Chemistry Research Topics被引用 6

一句话总结

本文提出一种通用的自适应三角化算法，用于平滑分子表面——特别是通过滚动球体和圆环方法生成的溶剂非覆盖（SES）和溶剂可及（SAS）表面。该方法根据局部曲率及与活性位点的距离动态调整网格分辨率，从而在不产生网格伪影的情况下保留细小表面特征（如狭窄通道），实现溶剂化能计算与分子可视化应用的精确建模。

ABSTRACT

In the given paper the algorithm describing original and universal principles of a triangulation of a smooth molecular surface: solvent excluding solvent (SES), received by primary and secondary rolling, and solvent accessible surface (SAS) is presented. These surfaces are a boundary between molecule and solvent. Originality of the given paper consists in creation of the universal and adaptive algorithm of a triangulation. Universality of algorithm of a triangulation consists that it is suitable for not only for a surface, received by rolling and consisting of fragments of torus and sphere, but for any smooth surface, including any level surface. Adaptability of this algorithm consists in facts that the mesh size of a triangulation can vary depending on its location; reflecting even small, but smooth features of a surface; preventing "jump" to close, but not neighbor sites of the surface, excepting "cut off" of narrow necks and channels. It is reached by either decreasing triangulation lattice step to value smaller than two principal radiuses of curvature of the molecular surface or decreasing triangulation lattice step close to the active centre - closed, but not neighbor sites of the surface. The received triangulated surface can be used for the demonstration purposes in molecular editors (the algorithm is applicable for a triangulation of any smooth surface, for example, level surfaces) together with for calculation of solvation energy and its gradients for continual models of solvent.

研究动机与目标

• 开发一种适用于任意平滑分子表面的通用三角化方法，包括通过滚动探针生成的SES和SAS表面。
• 解决在表面离散化过程中保留拓扑特征（如狭窄颈道和通道）的挑战。
• 实现随局部表面曲率和与活性位点距离而变化的自适应网格分辨率。
• 防止网格伪影（如连接至非相邻表面点的“跳跃”或细小结构细节被“截断”）。
• 同时支持分子编辑器中的可视化以及连续溶剂模型中溶剂化能及其梯度的定量计算。

提出的方法

• 该算法使用滚动探针方法，通过组合球面和环面表面片段生成SES和SAS表面。
• 基于主曲率半径自适应调整网格分辨率，当网格步长低于最小曲率半径两倍时采用更细的分辨率。
• 通过在活性中心附近细化网格，防止连接至非相邻表面点，即使这些点在几何上并不相邻。
• 采用通用框架，可扩展至任意平滑等值面，确保广泛适用性。
• 三角化过程通过检测并保留狭窄特征（如通道和裂隙）来保证拓扑保真度。
• 该算法专为集成至分子建模软件而设计，支持可视化与能量计算工作流的双重需求。

实验结果

研究问题

• RQ1 如何使三角化算法普遍适用于任意平滑分子表面（包括复杂的SES和SAS表面）？
• RQ2 何种自适应准则可确保网格分辨率捕捉细微表面特征而不引入伪影？
• RQ3 如何防止算法在保留狭窄结构特征的同时出现连接至非相邻表面点的网格‘跳跃’现象？
• RQ4 基于曲率的细化如何提升溶剂化能计算的准确性？
• RQ5 该算法是否可应用于分子表面之外的一般平滑等值面？

主要发现

• 该算法成功对复杂分子表面（包括具有狭窄通道和高曲率区域的表面）进行三角化，且未引入网格伪影。
• 基于曲率的自适应网格细化确保了在离散化过程中即使微小的平滑表面特征也能被保留。
• 通过在这些特征附近减小网格步长，即使在拓扑上不相邻的情况下，也能防止狭窄颈道和通道被截断。
• 该算法具有通用性，适用于任意平滑表面（包括等值面），不限于分子表面。
• 三角化后的表面既适用于分子编辑器中的可视化，也适用于连续溶剂模型中溶剂化能及其梯度的计算。
• 通过动态调整分辨率以匹配局部表面几何特征，该方法实现了高保真的表面表示。