この記事では、クライオ電子断層撮影法とそれに続く血球凝集素糖タンパク質のサブトモグラム平均を使用して画像化されたインフルエンザウイルスのデータ処理のためのプロトコルを紹介します。このプロトコルは、画像の前処理から最終的なモデルの改良まで、段階的なデータ処理をカバーします。
方法論記事
この記事では、クライオ電子断層撮影法とそれに続く血球凝集素糖タンパク質のサブトモグラム平均を使用して画像化されたインフルエンザウイルスのデータ処理のためのプロトコルを紹介します。このプロトコルは、画像の前処理から最終的なモデルの改良まで、段階的なデータ処理をカバーします。
クライオ電子断層撮影法は、不均一なサンプルを視覚化するための強力なツールであり、主な用途の 1 つは多形ウイルスの構造特性評価です。近年、これらの重要なタンパク質を無傷のビリオンの表面に直接可視化する方法として、ウイルス糖タンパク質のサブトモグラム平均化が登場しています。重要な標的の 1 つは、ウイルス エンベロープを密に覆い、インフルエンザ受容体の結合と膜融合に関与するインフルエンザ ウイルスの血球凝集素 (HA) 糖タンパク質です。インフルエンザHAのサブトモグラム平均が報告されていますが、クライオETに固有の信号対雑音比が低いことと、不均一なインフルエンザビリオンの分析に必要な手作業のために、その分解能は限られています。ここでは、インフルエンザビリオンの断層撮影データを効率的かつ堅牢に分析するために、いくつかのソフトウェアパッケージを統合したクライオET分析パイプラインを紹介します。このプロトコルは、最初の運動補正から最終的なモデル構築までのステップを通じて、インフルエンザビリオンからのHAの構造的決定を記述します。このパイプラインに続いて、A/プエルトリコ/8/34(PR8)インフルエンザ株から収集された2つのクライオETデータセットから、6.0 Åの分解能でのHA再構成が得られました。
クライオ電子断層撮影法(cryoET)は、タンパク質複合体、ウイルス、細胞、生物のスナップショットをキャプチャするために、過去数十年にわたって適用されてきました。クライオ電子顕微鏡(cryoEM)のモダリティであるクライオETは、生体サンプルを瞬間凍結し、傾斜1,2,3を介してさまざまな方向で画像化する構造生物学の方法です。次に、各方向で撮影された画像は、共通の傾斜軸に計算的に整列され、断層撮影に再構成されて3次元ビュー4を提供します。
X線結晶学や単粒子クライオEMでは、精製された構造的に均質な分子が必要ですが、クライオETでは、分子をそのネイティブコンテキスト内で直接画像化できます4。したがって、クライオETの主な利点の1つは、インフルエンザ5,6,7を含む膜ウイルスなどの多形サンプルを視覚化できることです。クライオETのもう一つの可能性は、スケールを超えて画像化する能力です。トモグラフィーは通常、5〜10 nmを超えて分解されませんが8、同じ粒子のコピーが識別され、整列され、平均化されるサブトモグラム平均化の統合により、リボソーム9,10などの一部の生体分子でほぼ原子分解能が得られる可能性があります。ただし、この分解能に到達できる分子の種類は限られています。サブトモグラムの平均は通常、10〜15 Åの解像度を超えません。対照的に、単粒子クライオEMは、分解能回転後に3〜4 Åの分解能を日常的に達成します11。クライオETデータ取得および解析ソフトウェアの両方の高スループットにおける最近の進歩により、ネイティブコンテキスト12、13、14、15、16、17、18内で追加の生体分子のサブナノメートル分解能の構造決定が可能になりました。
クライオETの一般的な用途の1つは、ウイルスの形態、組織、構造を視覚化することです。この技術では、単粒子クライオEMやX線結晶学に比べて分解能が低いにもかかわらず、クライオETとサブトモグラム平均化を組み合わせることで、ウイルスタンパク質が天然環境でどのように振る舞うかに関する情報を提供し、ビリオンの文脈におけるウイルスタンパク質の組織に関する重要な詳細を提供することができます。ウイルスのクライオETの一般的な標的は、治療薬やワクチンの主要な抗原や標的であることが多いため、宿主細胞の付着と融合に一般的に使用される表面糖タンパク質です。クライオET処理パッケージの最近の進歩により、これらの糖タンパク質19,20,21,22のサブナノメートル分解能平均を達成することがますます実現可能になってきています。そのような例の 1 つは、インフルエンザ ビリオンの表面にある主要なタンパク質である赤血球凝集素 (HA) です。このタンパク質は、受容体結合と膜融合の両方を行うだけでなく、単一のビリオン5に数百から数千のHAがあり、信じられないほど高密度にビリオンを覆っています。ここで紹介するプロトコル(図1)は、一般的に使用されるいくつかのパッケージを社内スクリプトと統合して、インフルエンザHAのサブトモグラム平均の前処理からモデル改良までの段階を描写します。
アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。
注:このプロトコルに使用されるサンプルデータセットは、このプロトコルに使用される2セットのチルトシリーズを含むEMPIAR-12864でアクセスできます。チルトシリーズは、2.09 Å/ピクセルの物理ピクセルサイズで精製されたインフルエンザAウイルスの手動で急落したグリッドから収集され、各チルトシリーズに複数のビリオンが含まれるように十分な視野を確保し、可能な限り最高解像度の再構成を実現します。ユーザー独自のデータセットの場合は、生のチルトムービーでワークフローを開始することをお勧めします。これらのデータセットは、高性能ワークステーションを使用して処理および視覚化されました。 資料表 には、このプロトコルに使用されるハードウェアとソフトウェアがリストされています。このプロトコルで使用されるすべてのソフトウェアパッケージはオープンソースであり、ダウンロードできます。インストールのリンクと手順 は、 材料表 に記載されています。クライオETデータセットを処理するために推奨されるワークステーションには、少なくとも8コアプロセッサ、6GBのVRAM、64GBのRAM、および2TBのローカルストレージを備えた専用GPUカードが装備されている必要があります。
1. Warp 23とIMOD 24におけるチルトムービーのデータ前処理とクライオ電子断層撮影の再構成
conda activate warp_environmentWarpTools create_settings --folder_data path/to/.tif --folder_processing warp_frameseries --output warp_frameseries.settings --extension “*.tif” --angpix 1.04 --gain_path gain_file.mrc --exposure 3.07WarpTools fs_motion_and_ctf --settings warp_frameseries.settings --m_grid 1x1x5 --c_grid 2x2x1 --c_range_max 7 --c_defocus_max 10 --c_defocus_min 4 --c_use_sum --out_averagesWarpTools ts_import --mdocs path/to/.mdoc --frameseries /path/to/frameseries --tilt_exposure 3.07 --min_intensity 0.3 --output tomostarWarpTools create_settings --folder_data tomostar --folder_processing warp_tiltseries --output warp_tiltseries.settings --extension “*.tomostar” --angpix 1.04 --gain_path gain_file.mrc --exposure 3.07 --tomo_dimensions NxNxNWarptools ts_stack --settings warp_tiltseries.settings --angpix 8.35WarpTools ts_import_alignments --settings warp_tiltseries.settings --alignments warp_tiltseries/tiltstack/ --alignment_angpix 8.35WarpTools ts_aretomo --settings warp_tiltseries.settings --angpix 8.35 --alignz 1000 --axis_iter 3 --exe AreTomo_executiveWarpTools ts_ctf --settings warp_tiltseries.settings --range_high 7 --defocus_min 2 --defocus_max 10 --auto_hand 4export WARP_FORCE_MRC_FLOAT32=1WarpTools ts_reconstruct --settings warp_tiltseries.settings --input_data input file names --angpix 8.35 --dont_invert2. 断層撮影の前処理と粒子ピッキング
conda activate isonet_envmkdir tomo_folder
mv tomograms*.mrc tomo_folder/isonet.py prepare_star tomo_folder --output_star tomograms.star --pixel_size 8.35isonet.py deconv tomograms.star --snrfalloff 0.7 --deconv_folder deconvolveconda activate eman_env
e2projectmanager.pycd path/to/tomogramse2spt_boxer_convnet.py --label label_namee2projectmanager.pyコマンドを使用して EMAN2 GUI を開きます。3. パーティクルキュレーション
4. 反復サブトモグラフィーの平均化と分類
WarpTools ts_export_particles --settings warp_tiltseries.setting --input_star pts2star.star --coords_angpix 8.35 --output_star bin4_export.star --output_angpix 8.35 --box 48 --diameter 140 --3drelion_convert_star --i bin4_export.star --o bin4_conv.starhead -n 30 bin4_conv.star >> subset.star & tail -n +31 bin4_conv.star | shuf -n 2000 >> subset.star
mpiexec -n 3 relion_refine_mpi --o init_ref/job001/run --auto_refine --split_random_halves --i subset.star --firstiter_cc --ini_high 20 --dont_combine_weights_via_disc --pool 3 --pad 2 --ctf --particle_diameter 300 --flatten_solvent --zero_mask --oversampling 1 --healpix_order 2 --auto_local_healpix_order 4 --offset_range 14 --offset_step 4 --sym C1 --low_resol_join_halves 40 --norm --scale --j 12 --gpu 0:1 --pipeline_control init_ref/job001mpiexec -n 3 relion_refine_mpi --o Refine3D/job001/run --auto_refine --split_random_halves --i bin4_conv.star --ref init_ref/job001/run_class001.mrc --firstiter_cc --ini_high 20 --dont_combine_weights_via_disc --pool 3 --pad 2 --ctf --particle_diameter 400 --flatten_solvent --zero_mask --oversampling 1 --healpix_order 2 --auto_local_healpix_order 4 --offset_range 16 --offset_step 4 --sym C1 --low_resol_join_halves 40 --norm --scale --j 12 --gpu 0:1 --pipeline_control Refine3D/job001relion_refine --o Class2D/job003/run --grad --class_inactivity_threshold 0.1 --grad_write_iter 200 --iter 200 --i Refine3D/job002/run_data.star --dont_combine_weights_via_disc --pool 3 --pad 2 --ctf --tau2_fudge 2 --particle_diameter 300 --K 20 --flatten_solvent --zero_mask -- strict_highres_exp 14 --center_classes --oversampling 1 --norm --scale --j 24 --skip_align --pipeline_control Class2D/job002。relion_star_handler --i input_file2.star --o output_good_class.star --select rlnClassNumber --minval goodclassnumber -maxval goodclassnumberrelion_star_handler --i "output_good_class1.star output_good_class2.star … output_good_classn.star" --o bin4_keep.star --combinerelion_star_handler --i input_file.star --o output_file.star --combinerelion_image_handler --i bin1_ref.mrc --o bin1_c3.mrc --sym c3MTools create_population --directory refine_m --name ha_finalMTools create_source --name source_1 --population refine_m/ha_final.population --processing_settings warp_tiltseries.settings
MTools create_species --population refine_m/ha_final.population --name ha_todaysdate --diameter 160 --sym c3 --temporal_samples 1 --half1 last_relion_refine/run_half1_class001_unfil.mrc --half2 last_relion_refine/run_half2_class001_unfil.mrc --particles_relion last_relion_refine/run_data.star --mask mask.mrcMCore --population refine_m/ha_final.population --refine_particlesMCore --population refine_m/ha_final.population --refine_particles --ctf_cs5. モデルの改良
アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。
この処理プロトコルの利用を実証するために(図1)、H1N1インフルエンザAウイルス株(A/プエルトリコ/8/1934)から得られた25枚の断層撮影を組み合わせた2つのデータセットに、前述のワークフローを適用しました。データ収集パラメーターの概要を 表 1 に示します。 図2 は、多形性インフルエンザビリオンの代表的な断層撮影とズームインビューを示しています。この断層撮影では、ビリオンが球形から楕円形/細長い形状までさまざまであるため、さまざまな形態が捉えられます。ほとんどのインフルエンザ粒子には、よく組織化されたM1およびvRNPアセンブリが含まれていますが、一部のビリオンはより無秩序で、重要な構造成分を欠いているようです。
このデータセットから、40,995のサブトモグラムの初期セットが、粒子のピッキングとキュレーション後のbin4(8.35 Å / pix)の...
アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。
重要なウイルスタンパク質の構造理解が深まることで、これらのウイルスに対する新しい治療法の発見を加速させることができます。過去10年間で、分解能革命により、単一粒子クライオEMを使用した高分解能ウイルス構造の決定が加速しましたが、この方法は、精製されたタンパク質または二十面体対称性を持つ非エンベロープウイルスのいずれかに限定されています。対照的に、クライオETは形態学的に多様な膜ビリオンを3Dで視覚化できますが、解像度は限られています。ここで紹介するパイプラインは、広く使用されているWarp-RELION-Mパイプライン9,32に基づいて構築され、一連のカスタムスクリプトをEMAN2の半自動畳み込みニューラルネットワークベースのパーティクルピッカーと統合して、エンドツーエンドのプロトコルを確立し、インフルエンザHAのサブナノメートル解像度の再構成を生成します。
クライオETにおけるより難しい問題の1つ...
アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。
著者らは開示するものは何もない。
著者らは、シファー・ラボとの有益な議論に感謝したいと思います。また、データ収集に協力し、サポートとアドバイスを提供してくれたマサチューセッツ大学チャンクライオEMコア施設にも感謝したいと思います。この研究は、国立総合医科学研究所の R01GM143773 から MS と R35GM151996 から CAS によって支援されました。
アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。
| 名前 | 会社 | カタログ番号 | コメント |
|---|---|---|---|
| AMD Ryzen Threadripper PRO 5965WX | AMD | https://www.amd.com/en/support/downloads/drivers.html/processors/ryzen-threadripper-pro/ryzen-threadripper-pro-5000wx-series/amd-ryzen-threadripper-pro-5965wx.html | |
| AreTomo 1.3.4 | UCサンフランシスコ | https://drive.google.com/drive/folders/1Z7pKVEdgMoNaUmd_cOFhlt-QCcfcwF3_ | |
| EMAN2 1952年2月99日 | ベイラー医科大学 | https://blake.bcm.edu/emanwiki/EMAN2 | |
| IMOD 4.12.27 | コロラド大学ボルダー校 | https://bio3d.colorado.edu/imod/ | |
| インフルエンザ分析スクリプト | UMass チャン医科大学 | https://github.com/jqyhuang/influenza-analysis | |
| IsoNet 0.3 | UCLA | https://github.com/IsoNet-cryoET/IsoNet | |
| M 2.0.0 | ジェネンテック | https://warpem.github.io/warp/home/m/ | |
| NVIDIA A4000 | NVIDIA | https://www.nvidia.com/en-us/products/workstations/rtx-a4000/ | |
| Open3D | インテル・ラボ | https://www.open3d.org/ | |
| フェニックス 1.21-5207 | ローレンス・バークレー国立研究所 | phenix-online.org | |
| RELION 4.0 | MRC分子生物学研究所 | https://relion.readthedocs.io/en/release-4.0/ | |
| Ubuntu 20.04 | Ubuntu(ウUbuntu | https://releases.ubuntu.com/focal/ | |
| UCSF キメラX 1.6.1 | UCサンフランシスコ | https://www.cgl.ucsf.edu/chimerax/ | |
| ワープ 2.0.0 | ジェネンテック | http://warpem.github.io/warp/ |
アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。
このJoVE記事のテキストまたは図の再利用許可をリクエスト
許可をリクエスト