[論文レビュー] Ab initio structure solution by charge flipping

G. Oszlányi, András Sütö|arXiv (Cornell University)|Aug 7, 2003

Enzyme Structure and Function被引用数 16

ひとこと要約

本論文は、高分解能X線回折データに対する、実空間と逆空間をフーリエ変換を用いて交互に繰り返す、シンプルなアバイト構造解決手法であるチャージフラッピングを紹介する。閾値δ未満の電荷密度の符号を反転させ、重みなしのFobsを用いて観測された構造因子の大きさを強制することで、事前の位相情報や対称性情報がなくても、小分子構造の解決に高い成功率を達成する。この手法は固定点ではなく、極限サイクルに収束する。

ABSTRACT

In this paper we present an extremely simple structure solution method termed charge flipping. It works ab initio on high resolution x-ray diffraction data in the manner of Fourier recycling. The real space modification simply changes the sign of charge density below a threshold, while in reciprocal space the modification is the Fobs map without any weighting. We test the algorithm using synthetic data for a wide range of structures, analyse the solution statistics and check the quality of reconstruction. Finally, we reconsider mathematical aspects of the algorithm in detail, showing that in this chaotic iteration process the solution is a limit cycle and not a fixed point.

研究の動機と目的

• 事前の化学的知識や位相情報が不要な、最小限で概念的にシンプルなアバイト構造解決手法の開発を目的とする。
• 理想化された高分解能の合成データを用いて、中心対称および非中心対称の多様な小分子構造に対して、チャージフラッピングアルゴリズムの実行可能性と頑健性を検証すること。
• 反復処理の数学的挙動、特に収束特性と解の性質を分析すること。
• 分離された精緻化プログラムを必要とせずに、構造内のすべての原子を回復できることを示すこと。

提案手法

• アルゴリズムはフリードル則を満たすランダムな位相セットから開始され、逆FFTを実行して初期の実空間における電荷密度を生成する。
• 各反復において、閾値δ未満の電荷密度ピクセルの符号を反転させ、δ以上のピクセルは変更しない。
• 変更された実空間密度をFFTにより逆空間に戻し、一時的な構造因子G(k)を計算する。
• G(k)の位相は保持されるが、その大きさは観測された構造因子振幅Fobs(k)に置き換えられ、F(0)は自由に変化可能である。
• 得られたF(k)を逆FFTで処理し、新しい実空間における電荷密度を生成することで、フーリエ再利用の1サイクルが完了する。
• 最適化は行わず、目的関数として用いるのではなくR因子のみをモニタリングする形で、繰り返し処理が行われる。

実験結果

リサーチクエスチョン

• RQ1 実空間におけるシンプルな符号反転ルールと、逆空間における重みなしFobs制約を組み合わせることで、アバイト結晶構造を効果的に解決できるか？
• RQ2 解の数学的性質は何か。具体的には、アルゴリズムが固定点に収束するのか、それとも極限サイクルに収束するのか？
• RQ3 閾値δの選択が収束、全電荷、および解の品質に与える影響は何か？
• RQ4 理想的なデータ条件下で、非中心対称構造を含む多様な小分子構造に対して、この手法がどの程度成功するか？
• RQ5 有限分解能およびノイズが加わった場合、アルゴリズムはどのように振る舞い、データの不完全性に対してどの程度頑健か？

主な発見

• チャージフラッピングアルゴリズムは、0.8 Å分解能の合成データを用いて、ケンブリッジ構造データベースに登録されたすべてのテスト小分子構造を、複数回の実行において高い成功率で解決した。
• 分離された精緻化ステップを必要とせず、テスト構造内のすべての原子が回復された。これは、本手法がアバイト構造解決において有効であることを示している。
• 全電荷F(0)は、位相のランダム化に起因して初期に過剰に増加した後、急速に減少し、長期間にわたり停滞した後、真の値のおおよそ30%低い値に収束した。
• 反復処理はカオス的であることが判明し、解は固定点ではなく極限サイクルに収束することが、位相の変化とR因子の挙動によって確認された。
• アルゴリズムの性能は、ポissonノイズや絶対スケール、イソトロープ的温度因子の誤差といったデータの不完全性に対しても頑健であったが、主な分析には理想データが用いられた。
• 収束後にδを小さくすることで、振動の抑制と全電荷の増加により解の品質が向上したが、δ = 0にするとゼロの海が消失し、過剰な電荷が生じた。