1. 新しい cryoSPARC v3 プロジェクトの作成とデータのインポート
注:データは、ポートランドのオレゴン健康科学大学(OHSU)で、Falcon 3直接電子検出器を搭載した300kVタイタンクリオス電子顕微鏡を使用して取得されました。画像は、129フレームにわたって28.38e-/Å2の総線量と-0.5μm~-2.5μmのデフォーカス範囲、1.045ÅのピクセルサイズでEPUを使用して計数モードで収集された。AAV-DJのサンプルは大洲のスタッフから提供された。
- Web ブラウザーで cryoSPARC v3 を開き、「 プロジェクト」 ヘッダーをクリックします。[ + 追加 ] を選択して、新しいプロジェクトを作成します。それに応じてプロジェクトにタイトルを付け、ジョブとデータが保存される既存のディレクトリへのパスを指定します。
- プロジェクトのワークスペースを作成するには、プロジェクトを開き、[ + 追加] をクリックして、[ 新しいワークスペース] を選択します。ワークスペースにタイトルを付け、[ 作成] をクリックします。
- 新しいワークスペースに移動し、右側のパネルで ジョブビルダー を開きます。このタブには、クライオスパーク v3 で使用可能なすべての機能が表示されます。「ムービー のインポート」 をクリックし、ムービーパスを指定し、リファレンスファイルパスを取得し、集録パラメータを次のように設定します: 生ピクセルサイズ1.045 Å、 加速電圧300 kV、 球面収差 2.7 mm、 総露光量 28.38 e-/Å^2。
- 「キュー」をクリックし、ジョブを実行するレーンとワークスペースを選択して、「作成」をクリックします。
メモ:集録パラメータはサンプルと顕微鏡に依存します。
2. CryoSPARC v3 - ムービーアライメントとCTF推定
- パッチモーション補正(マルチ)を開きます。このジョブでは、ステップ 1.3 でインポートしたムービーを入力として必要とします。ワークスペースでムービーの読み込みジョブカードを開き、Imported_movies出力を新しいジョブのムービープレースホルダにドラッグします。ジョブをキューに入れます。
メモ: この記事で概説されている cryoSPARC メソッドの詳細については、cryoSPARC チュートリアル16 を参照してください。
- CTF 推定を実行するには、[ パッチ CTF 推定 (複数)] を開きます。手順 2.1 で生成された顕微鏡写真を入力し、ジョブを キュー に入れます。
- 平均顕微鏡写真とCTF補正顕微鏡写真を検査し、さらに処理するサブセットを選択するには、 Curate Exposures を開き、ステップ2.2で得られた露出を入力します。ジョブを キュー に入れます。
- ジョブが 待機 モードに入ったら、ジョブ・カードの 「対話」 タブをクリックし、パラメーターしきい値を調整し、さらに処理するために個々の顕微鏡写真を承認または拒否します。よく一致する推定CTFと実験的CTFを含む顕微鏡写真(図2)を受け入れ、高い乱視、CTF適合度が低く、厚い氷を持つものは廃棄してください。
- 現在のデータの処理中に、乱視の上限を 400 Å、CTF 適合分解能を 5 Å、相対氷厚を 2 に設定します。[ 完了 ] をクリックして、下流処理用の顕微鏡写真を選択します。
3. CryoSPARC v3 - 手動およびテンプレートベースのパーティクルピッキング
- 手動ピッカーを開き、ステップ 2.4 から 2.5 で受け入れられたエクスポージャーを入力し、ジョブをキューに入れます。[インタラクティブ]タブをクリックし、[ボックスサイズ(px)]を300に設定し、複数の顕微鏡写真で数百個の粒子をクリックし、重複する粒子を選択しないようにします。ここでは、29枚の顕微鏡写真にまたがる340個の粒子が選択された。完了したら、[ピッキング完了]をクリックします。パーティクルを抽出します。
メモ: このプロトコルは、手動パーティクルピッキングを使用して、自動選択用のテンプレートを生成します。ただし、他の方法も利用できます17。
- 自動パーティクル ピッキングのテンプレートを生成するには、[ 2D 分類] をクリックし、ステップ 3.1 で生成されたパーティクル ピックを入力します。 2D クラス の数を 10 に変更し、ジョブを キュー に入れます。
- [2D クラスの選択] を開きます。ステップ 3.2 で取得したパーティクルとクラスの平均を入力し、[インタラクティブ]タブをクリックします。SNRが良好な代表的な2Dクラスを選択し、[完了]をクリックします。
注: クラス平均は、さまざまなパーティクル ビューを反映します。各ビューを反映するクラス平均を選択します。目標は、自動ピッキングのために標本のさまざまなビューを表す明確に定義されたテンプレートを作成することです。
- テンプレートピッカーを開き、手順 3.3 で選択した 2D クラスと、手順 2.4 ~ 2.5 の顕微鏡写真を入力します。粒子径 (Å) を 220 Å に設定し、ジョブをキューに入れます。
- 自動ピックを検査するには、[ パーティクル ピックの選択] を開き、手順 3.4 で生成されたパーティクルと顕微鏡写真を入力して、ジョブを キュー に入れます。
- [パーティクル ピックの選択] ジョブ カードで、[インタラクティブ] タブをクリックし、[ボックス サイズ] (px) を 300 に設定します。個々の顕微鏡写真をクリックし、粒子がはっきりと見えるまでローパスフィルタを調整し、正規化相互相関(NCC)しきい値を0.41、電力しきい値を54000~227300に設定します。
- いくつかの顕微鏡写真を検査し、必要に応じて、ほとんどの粒子が誤検知を含めずに選択されるようにしきい値を調整します。完了したら、[ ピッキング完了] をクリックします。パーティクルを抽出します。
注: 真のパーティクルは通常、NCC と検出力スコアが高く、それぞれテンプレートに似ており、SNR が高いことを示します。
- 顕微鏡写真 から抽出を開き、 ステップ3.7の顕微鏡写真と粒子を入力します。抽出された ボックス サイズ (px) を 300 に設定し、ジョブを キュー に入れます。
4. クライオスパーク v3 - 2D 分類
- [2D 分類] をクリックし、手順 3.8 から抽出した粒子を入力します。[2D クラスの数] を 50 に設定し、ジョブをキューに入れます。
- さらに処理に最適な 2D クラスを選択するには、[ 2D クラスの選択] を開きます。ステップ 4.1 で取得したパーティクルとクラスの平均を入力します。[ インタラクティブ] タブをクリックし、クラス内の解像度とパーティクル数に基づいて2Dクラスを選択します(図3)。成果物を含むクラスは選択しないでください。選択したら、[ 完了]をクリックします。
注: 通常、パーティクルを除去するには 2D 分類を複数回行う必要がありますが、パーティクルは明確に定義された個別のクラスに収束しません。2D 分類のラウンドを必要な回数だけ実行して、このようなパーティクルをデータ セットから削除します (図 3)。
5. CryoSPARC v3 - アブ・イニティオ ・レコンストラクションと均質なリファインメント
- 初期3Dボリュームを生成するには、 Ab-initio Reconstruction を開き、ステップ4.2または最終的な2D分類から得られたパーティクルを入力します。 対称性 を二十面体に調整します。ジョブを キュー に入れます。
メモ: 対称性はサンプルに依存するため、それに応じて変更する必要があります。不明な場合は、C1 対称を使用します。
- 均質な絞り込みを開きます。ステップ 5.1 の体積と 4.2 のパーティクル、または最終的な 2D 分類を入力します。対称性を変更し、ジョブをキューに入れます。ジョブが終了したら、フーリエシェル相関(FSC)曲線を調べ、UCSF Chimera18で調べるボリュームをダウンロードします。
6. クライオスパーク v3 からパーティクル座標をエクスポートし、PyEM を使用して RELION-3 にインポートする
注: パーティクル座標は、各顕微鏡写真内の個々のパーティクルの位置に関する情報を運びます。パーティクルスタックの代わりに座標をRELION-3に転送すると、そうでなければ利用できないリファインメントステップを実行できます。たとえば、パーティクル研磨では、初期ムービーフレームにアクセスする必要があります。したがって、クライオスパーク v3 から RELION-3 にパーティクル座標をエクスポートする前に、ムービーをインポートし、RELION-3 でモーション補正と CTF 推定を実行します。詳細については、RELION-3 チュートリアル 19 を参照してください。
- RELION-3 プロジェクトディレクトリに移動し、RELION-3 を起動します。
- ジョブタイプブラウザから インポート を開き、手順1.3のようにムービーと集録パラメータへのパスを指定します。
- モーション補正を実行するには、RELION-3 GUIからUCSF MotionCor220を使用し、 モーション補正 を開き、UCSF MotionCor2マニュアル21のようにデフォルトのパラメータを設定します。手順6.2で読み込んだムービーへのパスを入力します。[ モーション] タブで、 motioncor2 実行可能ファイルへのパスを指定します。
メモ: MotionCor2 は、複数の GPU を使用して並列に実行できます。
- CTFFIND-4.122 を使用して RELION-3 GUI を使用して CTF 推定を実行します。CTF推定を開き、ステップ6.3で生成された顕微鏡写真.starを入力します。「CTFFIND-4.1」タブで、CTFFIND-4.1 実行可能ファイルへのパスを指定し、RELION-3.1 チュートリアル19 のようにパラメーターを設定します。
- パーティクルスタックを cryoSPARC v3 から RELION-3 にインポートするには、まず cryoSPARC v3 からパーティクルスタックをエクスポートする必要があります。cryoSPARC v3 で、ステップ 4.2 の 2D クラス選択 ジョブのジョブ・カードまたは最終的な 2D 分類を開きます。[ 詳細] タブで、[ ジョブのエクスポート]をクリックします。 エクスポートジョブ は、particles_exported.csファイルを出力します。
- cryoSPARC v3 から RELION-3 にパーティクル座標をインポートする前に、ステップ 6.5 のparticles_exported.csファイルを .star 形式に変換する必要があります。PyEM23 を使用して、次のコマンドを実行して、particles_exported.cs ファイルを .star 形式に変換します: csparc2star.py particles_exported.cs particles_exported.star
- RELION-3 で、[手動ピッキング] タブをクリックし、[I/O] タブで、手順 6.4 で説明した CTF 絞り込みの顕微鏡写真を入力します。[表示]タブで、次のパラメータを入力します:粒子径(A):220、ローパスフィルタ(A):-1、CTF画像のスケール:0.5。 ジョブを実行します。ManualPick というディレクトリが RELION-3 ホームフォルダに生成されます。
メモ: この手順は、RELION-3 で手動ピッキングフォルダ構造を作成するために実行されます。手動ピッキングの実行中に、RELION-3 GUI でのピッキングに使用される平均顕微鏡写真ごとに、ピッキングされたパーティクルの座標を含む単一の .star ファイルが作成されます。
- ステップ 6.6 の particles_exported.star ファイルを含むフォルダーに移動し、自家製のスクリプトを実行して、クライオ SPARC v3 粒子のピッキングに使用される平均顕微鏡写真ごとに 1 つの manualpick.star ファイルを作成し、これがステップ 6.5 でエクスポートされた最終的な 2D 分類に貢献しました。結果の座標ファイルは、ManualPick/Movies フォルダーに保存されます。
- RELION-3 に戻り、 手動ピッキング ジョブを再度開きます。[ 続行]をクリックします。これにより、以前に cryoSPARC v3 で選択したパーティクルが RELION-3 GUI に表示されます。いくつかの顕微鏡写真を調べて、粒子座標の転送が達成されたかどうか、および粒子が適切に選択されているかどうかを確認します。
7. RELION-3 - 粒子抽出と2D分類
- 「粒子抽出」をクリックします。「I/O」タブで、ステップ 6.4 の CTF 補正顕微鏡写真とステップ 6.9 の座標を入力します。[抽出]タブをクリックし、[パーティクルボックスサイズ(pix)]を300に変更します。ジョブを実行します。
- 2D 分類を実行して、cryoSPARC v3 で生成されたパーティクル セットをさらにクリーニングし、高解像度の再構成を実現します。 [2D 分類] をクリックし、[ I/O ] タブで、手順 7.1 で生成された particles.star ファイルを入力します。「 最適化 」タブで、「 クラス数 」を「50」に、「 マスク直径 (A)」を「 280」に設定します。ジョブ を実行します 。
メモ: マスクはパーティクル全体を網羅する必要があります。
- 最適な 2D クラスを選択するには、[ サブセットの選択] メソッドをクリックし、手順 7.2 の _model.star ファイルを入力して、ジョブ を実行します 。手順 4.2 の説明に従ってクラスを選択します。
- 手順 7.2 と 7.3 を繰り返して、収束していないパーティクルを削除します。
8. RELION-3 - 3Dリファインメント、マスク作成、後処理
- cryoSPARC v3 で生成されたマップ (ステップ 5.2) を RELION-3 の 3D リファインメントの初期モデルとして使用します。「インポート」方法を選択し、「I/O」タブで「生のムービー/顕微鏡写真のインポート: いいえ、生の入力ファイル: ムービー/*.mrc」のパラメータを設定します。
- MTF ファイルを指定し、手順 1.3 の説明に従ってムービー集録パラメータを入力します。[ その他 ] タブで、入力ファイルとして cryoSPARC v3 マップを選択し、[ ノードの種類] を 3D 参照 (.mrc) に変更して、ジョブ を実行します 。
- 「 3D 自動調整 」を選択し、「 I/O 」タブで、「 入力イメージ」 をステップ 7.3 または最後の選択ジョブの particles.star ファイルとして設定します。クライオスパーク v3 再構築を リファレンス マップとして指定します。「 参照 」タブをクリックし、「 初期ローパス・フィルター (Å) 」を「50」に変更し、「 対称」 を「二十面体」に変更します。「 最適化 」タブで、「 マスク直径 (Å)」を「280」に変更し、ジョブ を実行します 。
- ランニングが終わったら、UCSFキメラでrun_class001.mrcを開きます。
- UCSF キメラで、[ ツール ]をクリックし、[ボリューム データ]で [ボリュームビューア]を選択します。これにより、音量設定を調整するための新しいウィンドウが開きます。 [ステップ ] を 1 に変更し、マップにノイズがないレベル値に達するまでスライダーを調整します。この値は、次のステップでマスクの作成に使用されるため、記録します。
- 自動精錬から生成されたマップは、周囲の溶媒からのノイズが解像度を低下させるため、真のFSCを反映していません。後処理の前に、試料と溶媒領域を区別するマスクを作成します。
- 「マスクの作成」をクリックし、ステップ 8.3 から入力run_class001.mrc をクリックします。
- 「マスク」タブをクリックし、次のようにパラメーターを調整します: ローパス・フィルター・マップ (Å): 10、ピクセル・サイズ (Å): 1.045、初期二値化しきい値: ステップ 8.5 で取得したレベル値、バイナリー・マップの拡張 この数ピクセル: 3、およびこのピクセル数のソフト・エッジの追加: 3. ジョブを実行します。
- UCSFキメラのマスクを調べます。マスクがきつすぎる場合は、「 バイナリマップを拡張」をこの数 ピクセル増やすか 、このピクセル数のソフトエッジを追加します。鋭いマスクはオーバーフィットにつながる可能性があるため、柔らかいエッジを持つマスクを作成することが重要です。
- 「後処理」をクリックし、「I/O」タブで、ステップ 8.3 で作成したハーフマップを入力し、8.6 からマスクします。[キャリブレーション済みピクセル サイズ] を 1.045 Å に設定します。シャープタブで、次のように入力します:B因子を自動的に推定しますか?:はい、自動B適合の最低解像度(A):10、独自のB因子を使用しますか?:いいえ。[フィルター] タブで、[FSC 重み付けをスキップしますか?] を [いいえ。ジョブを実行します] に設定します。
9. RELION-3 - 研磨トレーニングと粒子研磨
- 粒子ビームによって誘発される動きを補正する前に、まずトレーニングモードを使用してデータセットに最適なモーショントラックを特定します。 ベイジアン研磨 を開き、[ I/O ]タブで、手順6.3の動き補正顕微鏡写真、手順8.3のパーティクル、および手順8.8の後処理.starファイルを入力します。[ トレーニング ]タブをクリックし、次のパラメータを設定します: トレーニング最適パラメータ :はい、 テスト用のフーリエピクセルの割合:0.5、 この多くのパーティクルを使用:5000。ジョブ を実行します 。
メモ: このスクリプトは、次の手順opt_params_all_groups.txt実行するために必要な最適化された研磨パラメータを含む RELION-3 Polish フォルダにファイルを生成します。
- トレーニングジョブが終了したら、[ベイジアンポリッシング]をクリックします。[トレーニング]タブをクリックし、[最適なパラメータのトレーニング?]を[いいえ]に設定します。「ポーランド語」タブを選択し、「最適化されたパラメーター・ファイル」で、ステップ 9.1 のopt_params_all_groups.txtファイルへのパスを指定します。[実行]をクリックします。
- 研磨された粒子のセットで3D精錬(ステップ8.3)と後処理(ステップ8.8)を繰り返します。
10. RELION-3 - CTFおよび粒子ごとの微細化
- 高次収差を推定するには、[ CTF リファインメント ]を開き、[ I/O ]タブの[パーティクル]で、最近のリファイン3Dジョブ(run_data.star)で研磨された パーティクルを含む.starファイルへのパスを選択します。
- 「後処理 スター・ファイル」で、最新の後処理ジョブからの出力へのパスを設定します (ステップ 9.3)。
- 「適合」(Fit) タブを選択し、次のパラメータを設定します: 推定値 (異方性倍率): いいえ、CTF パラメータ適合を実行しますか?いいえ、ビームチルトを推定:はい、トレフォイルも推定しますか?はい、4次収差を推定しますか?はい。ジョブを実行します。
- 前のジョブ(particles_ctf_refine.star)で生成されたパーティクル(Refine3Dから)を入力として使用して、手順10.1を繰り返します。「適合」タブで、「推定値 (異方性倍率)」を「はい」に変更し、ジョブを実行します。
- 前のジョブ(particles_ctf_refine.star)で生成された (Refine3Dから)パーティクル を入力として使用して、手順10.2を繰り返します。[ 適合] タブで、次のパラメータを設定します: 推定値(異方性倍率):いいえ、 CTFパラメータフィッティングを実行?:はい、 フォーカスを合わせる?:パーティクルごと、 乱視をフィット? 顕微鏡写真ごと、 B因子を適合?: いいえ、 位相シフトを適合:いいえ、 ビームチルトを推定?: いいえ、 4次 収差を推定?: いいえ、 実行 します。
注記: パーティクルのコントラストが十分であれば、[ 乱視に適合? ]タブを パーティクルごとに設定できます。このデータセットでは、 パーティクルごとの 乱視の絞り込みでは、マップの品質と解像度は向上しませんでした。
- 手順 10.3 のパーティクルで 3D リファインメントを繰り返し、[I/O] タブで [ 溶剤で平坦化された FSC を使用?] を [ はい] に設定します。実行が終了したら、後処理ジョブを実行し(ステップ8.8)、UCSFキメラでマップを調べる(ステップ5.2)。
11. RELION-3 パーティクル座標と 3D マップを Scipion 3 に転送する
- RELION-3 マップをさらに絞り込んで検証するには、まず最後の後処理ジョブ (ステップ 10.4) からボリュームとパーティクルを Scipion 3 にインポートします。Scipion 3 を起動し、新しいプロジェクトを作成します。
- 左側の[プロトコル]パネルで、[インポート]ドロップダウンを選択し、[パーティクルのインポート]をクリックします。 次のパラメータを変更します: インポート元: RELION-3、スターファイル: postprocess.star、およびステップ 1.3 のように取得パラメータを指定します。「実行」をクリックします。
- [インポート]ドロップダウンをクリックし、[ボリュームのインポート]を選択します。[インポート元] で、RELION-3 マップへのパスを指定します。ピクセルサイズ(サンプリングレート)Å/pxを1.045に変更して実行します。
12. Scipion 3 - 高解像度の改良
- まず、グローバルアライメントを実行します。[プロトコル] パネルの [絞り込み] ドロップダウンを選択し、[Xmipp3 - highres24] をクリックします。ステップ 11.2 および 11.3 からインポートしたパーティクルとボリュームをそれぞれフルサイズ イメージおよび初期ボリュームとして入力し、対称グループを 2 面体に設定します。[角度割り当て]の下の[画像の配置]タブで、[グローバル]を選択し、[最大ターゲット解像度]を3 Åに設定して、ジョブを実行します。
- ジョブが終了したら、[ 結果の分析] をクリックします。新しいウィンドウで、 UCSF キメラ アイコンをクリックして、洗練されたボリュームを調べます。さらに、[ 表示解像度プロット( FSC)]をクリックして絞り込み後にFSCがどのように変化したかを確認し、[角度 変更によるヒストグラムをプロット] をクリックしてオイラー角度の割り当てが変更されたかどうかを確認します。
メモ: 入力 RELION-3 構造体の解像度によっては、[角度割り当て] タブの [最大ターゲット解像度] に異なる値を設定して、この手順を複数回繰り返すことができます。詳細については、Scipion チュートリアル25 を参照してください。
- 手順 12.1 をローカル位置合わせで繰り返します。前のジョブをコピーし 、[前の実行の選択] を [Xmipp3 - 高解像度グローバル] に変更します。[ 角度の割り当て ] タブで、[ 画像の配置 ] を [ ローカル] に変更します。 最大ターゲット解像度 を 2.1 Å に設定します。
- UCSFキメラの洗練されたマップを調べ、FSCと角度割り当ての変化を分析します(ステップ12.2)。解像度が向上しず、角度割り当てが収束するまでローカル絞り込みを繰り返し、必要に応じて 最大ターゲット解像度 を調整します。
- Scipion 3 からの出力マップは、Phenix26 でさらに密度を変更し、シャープにすることができます。
13. Scipion 3 - マップ検証
- Xmipp3 - 高解像度で生成された最終的なマップのローカル解像度を調べます。Xmipp - ローカル MonoRes27 を開き、前のジョブの最終ボリュームとステップ 8.6 で生成されたマスクを入力します。解像度範囲を1~6Å、間隔を0.1Åに設定し、ジョブを実行します。
- 実行が終了したら、[ 結果の分析] をクリックし、解像度ヒストグラムと解像度で色分けされたボリュームスライスを確認します。
- パーティクルが適切にアライメントされているかどうかを確認するには、[ マルチリファレンス アライナビリティ]28を開き 、ステップ12.3のパーティクルと体積を入力します。 結果の分析 をクリックして、検証プロットを表示します。理想的には、すべてのポイントを周囲にクラスター化する必要があります (1.0,1.0)。
- Xmipp3 - オーバーフィットの検証を開きます。ステップ 12.3 のパーティクルとボリュームを入力します。実行が終了したら、結果を分析し、オーバーフィットプロットを調べます。整列ガウスノイズと整列パーティクル曲線の交差は、オーバーフィットを示します。