[論文レビュー] Molecular multiresolution surfaces

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

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

ひとこと要約

本稿では、van der Waals 表面における連続的溶媒の拡散を用いて、さまざまなプローブ半径で得られる表面を生成する統合的で多スケールなフレームワークとして、分子マルチリゾリューション表面を導入する。局所スペクトル的進化カーネル（LSEK）を用いた高速な1ステップ数値積分により、特異性の問題を回避しつつ、溶媒を排除した表面、溶媒にアクセス可能な表面、van der Waals 表面を高精度に捉えることが可能である。サイクロヘキサンおよび大規模なタンパク質を用いた検証で、既存の手法と高い一致を示した。

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.

研究の動機と目的

• 異なるプローブ半径におけるバイオ分子表面の統合的で多スケールなフレームワークの開発。
• 特に特異性や自己交差の問題を含む、従来の分子表面計算におけるトポロジカルおよび計算的課題の克服。
• van der Waals 表面、溶媒にアクセス可能な表面、溶媒を排除した表面の連続的で幾何学的に一貫した記述の提供。
• 新規な数値積分スキームを用いた、マルチリゾリューション表面の効率的かつ1ステップでの計算の実現。
• 複雑なバイオ分子を用いた検証を通じて、精度とスケーラビリティの両立を示すこと。

提案手法

• 曲率に依存する拡散係数を有する反応拡散方程式を用いて、van der Waals 表面における溶媒の拡散をモデル化する。
• 溶媒密度の不均一性およびフラックス相関を考慮するため、高次項（q > 1）を含む一般化されたフィックの法則を採用する。
• 初期条件を生成するため、溶媒にアクセス可能な表面（SAR）領域では rho(r,0) = 0、それ以外では rho0 を初期値として設定する。
• 局所スペクトル的進化カーネル（LSEK）を用いて、O(N) の計算量と制御可能な空間精度を達成しつつ、1ステップで拡散方程式を解く。
• LSEK はエルミート関数とウィンドウ付き畳み込みを用いて密度の進化を計算し、パラメータとして Mh=88、Mα=32、σα=3.05hα を使用する。
• マルチリゾリューション表面は、進化した溶媒密度の等高面として抽出され、van der Waals 表面はプローブ半径がゼロの極限として定義される。

実験結果

リサーチクエスチョン

• RQ1 1つの連続的フレームワークが、van der Waals 表面、溶媒にアクセス可能な表面、溶媒を排除した表面の表現を統合的に可能にするか。
• RQ2 分子表面の多スケール的性質が、物理的拡散過程によってどのように捉えられるか。
• RQ3 高速で1ステップの数値的手法が、反復的または複雑な幾何的アルゴリズムに代わって分子表面生成に適用可能か。
• RQ4 提案手法が、従来の表面計算に内在するトポロジカル特異性や自己交差を回避できるか。
• RQ5 大規模なバイオ分子へスケーリングする際、高い精度を維持しながら効率的であるか。

主な発見

• 本手法は、さまざまなプローブ半径において分子マルチリゾリューション表面を効果的に生成でき、van der Waals 表面はプローブ半径がゼロの極限として得られる。
• サイクロヘキサンの例では、計算された溶媒を排除した表面が、MSMSアルゴリズムと3つの断面で優れた一致を示した。
• 細胞分裂タンパク質の中間分子表面は、主要な構造的対称性を捉えつつ、細かな非本質的詳細を抑制した。
• LSEK 法は、固定された精度で O(N) の計算量を達成し、拡散方程式の1ステップ積分を効率的に実現した。
• 埋め込まれた原子表面と埋め込まれていない表面を区別する必要がなく、他の手法で一般的な特異性を回避した。
• 数値的検証により、小分子および数千原子を有する大規模タンパク質に対しても、本手法の精度と頑健性が確認された。