[論文レビュー] Expectation-maximization for structure determination directly from cryo-EM micrographs

Shay Kreymer, Amit Singer+1 more|ArXiv.org|Feb 24, 2023

Advanced Electron Microscopy Techniques and Applications被引用数 4

ひとこと要約

本論文は、粒子ピッキングを経由せずに、低温電子顕微鏡（cryo-EM）の顕微鏡写真から直接3次元分子構造を再構築するための確率的近似期待最大化（SAEM）アルゴリズムを提案する。未知の粒子位置および方位を周辺化することで、低SNR（信号対雑音比）条件下でも高分解能の構造回復を達成し、SNR=0.75の31³ボクセルのBPTI変異体を、SNR=3.5のより小さな11³ボクセルの体積を高精度で再構築した。

ABSTRACT

A single-particle cryo-electron microscopy (cryo-EM) measurement, called a micrograph, consists of multiple two-dimensional tomographic projections of a three-dimensional (3-D) molecular structure at unknown locations, taken under unknown viewing directions. All existing cryo-EM algorithmic pipelines first locate and extract the projection images, and then reconstruct the structure from the extracted images. However, if the molecular structure is small, the signal-to-noise ratio (SNR) of the data is very low, making it challenging to accurately detect projection images within the micrograph. Consequently, all standard techniques fail in low-SNR regimes. To recover molecular structures from measurements of low SNR, and in particular small molecular structures, we devise an approximate expectation-maximization algorithm to estimate the 3-D structure directly from the micrograph, bypassing the need to locate the projection images. We corroborate our computational scheme with numerical experiments and present successful structure recoveries from simulated noisy measurements.

研究の動機と目的

• 特に小分子構造に対して顕著な低SNR領域におけるcryo-EMの粒子ピッキングの根本的限界を克服すること。
• 不要なパラメータ（位置および方位）を明示的に推定するのではなく、周辺化することで3次元構造決定の一貫性のある推定器を開発すること。
• 特に40 kDa未満の分子に対して、事前の粒子抽出を経ずに顕微鏡写真から直接高分解能の3次元再構築を可能とすること。
• さまざまなSNRおよび体積サイズのシミュレートされた顕微鏡写真を用いて手法を検証し、耐障害性および分解能の忠実性を示すこと。

提案手法

• 本手法は、未知の粒子位置および方位を統合しながら、3次元構造を反復的に推定するための確率的近似期待最大化（SAEM）フレームワークを採用する。
• 2次元投影が未知の回転および平行移動の下でどのように形成されるかを考慮した尤度モデルを用い、顕微鏡写真のパッチ上での確率的モデルとして定式化する。
• Eステップでは、有限個の回転角の集合に対してモンテカルロサンプリングを用いて、ポーズパラメータ（回転およびシフト）の後方分布を近似する。
• Mステップでは、球面調和係数を最大次数ℓ_maxまで用いて、期待対数尤度を最大化することで3次元構造推定値を更新する。
• アルゴリズムは反復的に適用され、3次元ボリューム推定値を改善し、FSC曲線を用いて収束をモニタリングする。
• 粒子ピッキングに類似した手順は、シフトに関する周辺尤度を計算することで副次的に得られ、推定された構造から位置推定が可能になる。

(b)

実験結果

リサーチクエスチョン

• RQ1 低SNR領域において、特に粒子ピッキングを経由せずに顕微鏡写真から直接3次元cryo-EM構造を再構築することは可能か？
• RQ2 未知の粒子位置および方位を周辺化することで、一貫性があり高分解能の3次元再構築が可能になるか？
• RQ3 SNRが極めて低い状況下で、直接再構築法と粒子ピッキングベースのパイプラインの性能を比較するとどうなるか？
• RQ4 提案されたSAEMアルゴリズムは、ノイズの多い顕微鏡写真から、細かな構造的詳細（例：ℓ_max=10まで）をどの程度回復できるか？
• RQ5 主たる粒子ピッカーとしての目的ではないが、尤度に基づくシフト推定を用いて間接的に粒子位置を同定できるか？

主な発見

• SNR=0.75の顕微鏡写真から、31³ボクセルのBPTI変異体構造を高分解能の忠実性を保ちながら再構築に成功した。
• より小さな11³ボクセルのBPTIボリュームでは、SNR=3.5で正確な再構築が達成され、視覚的およびFSCベースの検証で真の構造と強い一致を示した。
• 31³ボクセル再構築のFSC曲線は、ℓ_max=10の分解能に達しており、高周波数の構造的詳細の回復を示している。
• 従来の自己相関に基づく手法に比べ、分解能とノイズ耐性の両面で優れており、以前の研究でℓ_max=2であったのに対し、本手法ではℓ_max=10まで回復できた。
• 尤度最大化に基づく間接的な粒子ピッキング手順は、SNRが十分に高い場合には真の粒子シフトを的確に同定でき、手法の内部的一致性を示した。
• 本アルゴリズムは、さまざまな体積サイズおよびSNRレベルにおいても耐障害性を示し、厳しい条件下でも直接3次元再構築が可能であることを確認した。

(e)