PowerPoint プレゼンテーション

Size: px

Start display at page:

Download "PowerPoint プレゼンテーション"

ふじきみさだい
6 years ago
Views:

1 SHELX による精密化の紹介量子科学技術研究開発機構 (QST) 量子ビーム科学研究部門東海量子ビーム応用研究センター平野優 1

2 Outline 1. SHELX の概要 2. SHELXL による精密化 3. 実例 (NADH シトクロム b 5 還元酵素を例にして ) 参考文献等 SHELX HP: 結晶構造精密化 SHELXLの使い方大場茂植草秀裕著三共出版 SHELXLの特徴を知ろう! 大場茂日本結晶学会誌 57, 87 (2015) 2

3 1. SHELX の概要 SHELX は X 線回折および中性子回折による低分子と高分子の結晶構造解析のためのプログラム群アカデミックフリー (Linux, Win, Mac OSX に対応 ) SHELX-2017 に含まれるプログラム SHELXT: 低分子の構造決定 SHELXS: 低分子の構造決定 ( 直接法 ) SHELXL: 低分子と高分子の精密化 PDB2INS:SHELXLのインプットファイル作成 CIFTAB, ShredCIF:SHELXLのCIFファイル編集 SHELXC, SHELXD, SHELXE: 高分子の位相決定 AnoDe: 異常散乱電子密度の作成等この後の実習 SHELX HP より : 3

4 Number of Entries タンパク質の精密化に用いられるプログラム SHELX による精密化の登録数は 2361 X 線結晶構造解析の約 1.8% Refinement software used PDB HP より : 4

5 SHELXL による精密化の特徴高分解能データの精密化に適している (2 Å より高分解能でできるだけ高分解能が望ましい ) ディスオーダー非等方性 ADP(atomic displacement parameter) 水素原子をモデルに含めた精密化が可能結合距離結合角度の esd 値 (estimated standard deviation) が得られる利用例 :Asp, Glu 側鎖のプロトン化状態 Charged O C O pka Protonated O C OH 1.25 Å 1.25 Å 1.21 Å 1.32 Å 3.8 (Asp) 4.3 (Glu) C O 距離同じ異なる (~0.1 Å) 5

6 完全行列最小二乗法による精密化モデル構造の全ての変数 p i = 10 N N = 原子数 ( 座標 3 異方性温度因子 6 占有率 1) R = hkl w hkl F o hkl 2 F c hkl 2 2 を最小にする条件 R p i = 0 Taylor 展開によって近似して連立一次方程式 (normal equation) を求める i ε i hkl w hkl F c hkl 2 p i F c hkl 2 p j = hkl w hkl F o hkl 2 F c hkl 2 F c 2 p j a ij b j i ε i a ij = b j (e i : p i のシフト量 ) 全ての行列要素について計算 p j の標準偏差 σ p j = a ij 1 R N o N p (N o, N p : データ数, パラメータ数 ) (A B 間距離 l については σ l = σ 2 A + σ 2 B ) 6

7 結合距離の esd に対する制約の影響 Concanavalin A(0.94 Å 分解能 ) における完全行列最小二乗法による精密化 Restrained Unrestrained 誤差の見積り過小評価 C C C N C O Restrained あり精密化の際の結合距離の s 2 1 σ res l = 2 1 σ diff 制約なし l σ geom l 結合距離の制約 International Tables for Crystallography (2012). Vol. F, p. 503, Figure を改変 7

8 2. SHELXL による精密化 Refmac, PHENIX などのプログラムによる精密化 (~1.5 Å 分解能等方性温度因子 ) SHELXL へ移行 (Rigid body refinement) 分解能の拡張 ( 共役勾配法による精密化 ) ( マルチコンフォメーションの追加非等方性温度因子制約の除外 ) 水素原子をモデルに付加して精密化完全行列最小二乗法による esd の計算 8

9 精密化において必要なファイル 1. インプットファイル (ins) TITL Rubredoxin in P1 (from 6RXN in PDB) CELL ZERR LATT -1 SFAC C H N O S FE UNIT DEFS RESI 1 MET C O N コマンド原子座標占有率温度因子 HKLF 4 END 2. 反射データファイル (hkl) (h k l I o s(i o ) flag) 9

10 精密化に必要なファイル作成方法 1. ins ファイル PDB2INS を利用 pdb2ins Create.hkl file from structure factor file (cif) or PDB code? (y or n) [N]: N Enter name of PDB file to read (To download a PDB file enter '@<PDBCODE>'): start-model.pdb Enter HKLF code (3 for F, 4 for F-squared) [4]:4 Please enter Z (number of molecules per cell):4 No wavelength found in file. Enter wavelength in Angstroms [ ]: 0.7 Reset water occupancy to unity? (y or n) [Y]: Y Enter.ins filename to be created [start-model.ins]: ins 作成時 N 初期構造 4(default) 選択 Z 値測定波長水分子の占有率 insファイル名 2. hkl ファイル mtz2hkl を利用よりダウンロード可能 mtz2hkl data1.mtz hkl ファイル出力 (HKLF 4 フォーマット ) CCP4 の mtz2various を利用 Format: SHELX, Label: I, s(i), FreeR-flag 10

11 精密化における入出力の流れ 1. 入力 2. 実行 3. 出力 insファイル : 精密化のコマンド原子座標データ hklファイル : 反射データ ins ファイル hkl ファイルの拡張子前部分を揃える ( 例 : refine-1.ins, refine-1.hkl) shelxl refine-1 ( 並列計算利用の場合 t オプション -t4 など ) resファイル :insファイルと同じ形式精密化後の原子座標データ lstファイル : 精密化のログ結合距離角度 (esd) なども書き出される pdbファイル : 精密化後のpdb fcfファイル : 精密化後の反射リスト (Cootで電子密度図の描写に利用) ( 例 : refine-1.res, refine-1.lst, refine-1.pdb, refine-1.fcf) Coot で読み込み 11

12 Coot におけるファイルの読み込み File Open Coordinates res ファイル選択 File Open MTZ, mmcif, fcf or phs fcf ファイル選択モデル修正後 File Save Coordinates で新しい ins ファイルが保存される 12

13 ins ファイルの内容 1 コマンドの種類約 80 コマンド + 引数という形式になっている TITL converted from file start-model.pdb CELL ZERR TITL タイトル (76 文字以内 ) CELL l, a, b, c, a, b, g ZERR Z, esd(a), esd(b), esd(c), esd(a), esd(b), esd(g) REM Space group P LATT -1 SYMM x+0.5, -y+0.5, -z SYMM -x, y+0.5, -z+0.5 SYMM -x+0.5, -y, z+0.5 REM LATT SYMM コメント行空間格子の型単純格子は 1 対象心の無い場合( タンパク質 ) は負対象操作 13

14 ins ファイルの内容 2 SFAC C H N O P S DISP $C DISP $H UNIT SFAC モデルに含まれる原子種 DISP f, f, 吸収断面積 (Cu Ka, Mo Kaなど使用時は省略可 ) UNIT 結晶格子中の原子数 DEFS CGLS 20-1 SHEL ANIS FMAP 2 PLAN LIST 6 WPDB 2 DEFS 制約の標準偏差のdefault 値 (DFIX, CHIV, DELU, SIMU) CGLS 共役勾配法による精密化 ( サイクル数 -1の反射をR free の計算に使用 ) SHEL 分解能範囲 ANIS 非等方性温度因子の適用 FMAP フーリエ合成の種類 ( 2 はF o -F c ) PLAN フーリエ合成で探すピーク数 O, N 等の原子からの距離 (2.3 Å<d<4.0 Å) LIST 反射リスト (.fcf) の種類 ( 6 の時 Cootで読込み可 ) WPDB 1 ( 等方性 ) 2 ( 非等方性 ) 負の時水素をPDBファイルに書出し 14

15 ins ファイルの内容 3 RIGU DELU $C_* $N_* $O_* $S_* $P_* SIMU 0.1 $C_* $N_* $O_* $S_* $P_* XNPD SWAT MERG 4 ISOR_HOH 0.1 $O CONN_HOH 0 O RIGU,DELU,SIMU,ISOR 温度因子の制約 XNPD 温度因子が非正定値にならないように下限値を設定 SWAT バルク溶媒補正 MERG 反射データの平均化 4 の時フリーデル対を含め平均化 CONN 非水素原子間の結合を設定上記の場合水分子とタンパク質間の結合を回避 y x z DELU RIGU SIMU ISOR DU zz = 0 DU zz = 0, DU yz = 0, DU xz = 0 15

16 ins ファイルの内容 4 RTAB_* Omeg CA_+ N_+ C CA RTAB_* Phi C_- N CA C RTAB_* Psi N CA C N_+ REM HFIX_ALA 43 N REM HFIX_ALA 13 CA REM HFIX_ALA 33 CB CHIV_ALA C CHIV_ALA CA DFIX_ALA C O DANG_ALA C N RTAB 指定した原子間の距離角度等をlstファイルに出力 HFIX 水素原子の発生一桁目の3はriding model 使用時 CHIV キラル体積 ( 指定した原子に結合する4 原子が形成する四面体の体積 ) の制約 DFIX 結合距離の制約 DANG 結合角度の制約 HFIX 43 HFIX 13 HFIX 33 CHIV DANG H H H H d (Å) N C H C CA C N H CA 16

17 ins ファイルの内容 5 WGHT FVAR RESI A:20 ASP N = WGHT 反射の重み付け係数 AFIX w43 H0 2 h = 1 σ 2 I o + ap 2 + bp AFIX 0 (default f = 1/3) PART 1 CA = FVAR スケール因子 (1 番目 ) 自由変数(1 番目以外 ) PART 2 CA = PART 0 C = HKLF 4 END P = f max 0 or I o + 1 f I c (default a = 0.1, b = 0) 17

18 ins ファイルの内容 5 WGHT FVAR RESI A:20 ASP N = AFIX 43 H AFIX RESI 0 アミノ酸残基 ( リガンド水分子等 ) の指定 PART 1 A:20 Chain A, 20 番目の残基 (Chain IDはSHELX-2017で対応 ) CA Coot v で読み込み可能 = 以前のバージョンの SHELXでは (chain Aは +1000) などとしていた PART 原子名 2 SFAC 番号 x, y, z (fractional) 占有率 U CA (iso のときU iso ), U 22, U , U 23, U = 13, U 12 = は改行後の行も同一行とみなす占有率のは1に固定を示す PART 0 C = HKLF 4 END 18

19 ins ファイルの内容 5 WGHT FVAR RESI A:20 ASP N = AFIX 43 H AFIX 0 PART 1 CA = AFIX 43~AFIX 0 直前の原子に対する水素原子の取り扱いを指定 PART 原子名 2 SFAC 番号 x, y, z (fractional) 占有率 U iso CA 1 U iso のは直前の原子の U eq の1.2 倍となるような制約 = OH やCH 3 などでは倍とする PART 0 U C eq = 1 13 U 11 aa U 22 bb U 33 cc U 12 a b cos γ + 2U 13 a c cos β + 2U b c cos α = HKLF 4 END 19

20 ins ファイルの内容 5 WGHT FVAR RESI A:20 ASP N = AFIX 43 H AFIX 0 PART 1 CA = PART 2 CA = PART 0 C = PART 1,PART 2,PART 0 ディスオーダー ( マルチコンフォメーション ) 部分の指定 HKLF 4 END 原子の占有率はFVARの2 番目の値を参照することを示すは (1-FVAR(2)) を示し占有率の合計が1となるようにしている 20

21 ins ファイルの内容 5 WGHT FVAR RESI A:20 ASP N = AFIX 43 H AFIX 0 PART 1 CA = PART 2 CA = HKLF 反射データファイルフォーマットの指定 PART 0 4 の時は (h k l I o s(i o ) flag) C = HKLF 4 END 21

22 lst ファイルの内容読み込んだ ins ファイルの内容原子間結合のリスト最小二乗法の各サイクルのまとめ最終サイクル後の座標温度因子最終サイクル後の R 値,S 値等の結果結合距離角度等のリスト S 値 (Goodness of fit) S = hkl w hkl F o hkl 2 F c hkl 2 2 N o N p フーリエ合成で検出したピークのリスト res ファイルの末尾にも出力 22

3. 実例 (NADH シトクロム b 5 還元酵素を例にして ) NAD + NADH 酸化型酸化型二電子還元型 E-FAD E-FAD - NADH E-FADH - NAD + b5r NADH 結合ドメイン H - transfer FAD 結合ドメイン e - (Cyt b 5 ) E-FAD NAD + 再酸化型 E-FAD - NAD + 一電子還元型 e - (Cyt b H +

23 3. 実例 (NADH シトクロム b 5 還元酵素を例にして ) NAD + NADH 酸化型酸化型二電子還元型 E-FAD E-FAD - NADH E-FADH - NAD + b5r NADH 結合ドメイン H - transfer FAD 結合ドメイン e - (Cyt b 5 ) E-FAD NAD + 再酸化型 E-FAD - NAD + 一電子還元型 e - (Cyt b H + 5 ) (Red semiquinone) E-FADH NAD + 一電子還元型 (Blue semiquinone) アミノ酸残基数 :272 酸化型は FAD を結合非対称単位中に 1 分子 X-ray crystal structures of b5r from porcine liver M. Yamada et al., J. Mol. Biol., 425, (2013). Å, Å K. Takaba et al., Sci. Rep., 7, (2017). Å 23

b5r の X 線回折データ統計値 Resolution range (Å) 50-0.

87) Space group P2 1 2 1 2 1 Cell parameters

24 b5r の X 線回折データ統計値 Resolution range (Å) ( ) Space group P Cell parameters (Å) a Å b 72.1 c 84.8 Completeness (%) 100 (100) I/s(I) 76.1 (4.4) R meas (%) 6.0 (45.5) PF BL-5 において回折データ収集 Wavelength 0.76 Å Detector ADSC Q315 Distance 100 mm 24

25 b5r の精密化の流れ PHENIX プログラムによる精密化 ( Å 分解能等方性温度因子 ) SHELXL へ移行 ( Å 分解能 Rigid body refinement) 剛体の範囲を AFIX 6 ~ AFIX 0 で囲む分解能の拡張マルチコンフォメーションの追加 ( Å 分解能等方性温度因子 ) 分解能の拡張非等方性温度因子の適用 ( Å 分解能 ) 制約の除外水素原子の導入 ( Å 分解能 ) 完全行列最小二乗法による esd の計算 ( Å 分解能 ) R work = 18.7%, R free = 23.0% R work = 17.9%, R free = 20.4% R work = 13.9%, R free = 16.9% R work = 13.1%, R free = 15.8% 25

マルチコンフォメーションと水素原子の取り扱い Asn198 RESI A:198 ASN Od1 PART 1 PART 1 PART 2 N 3 0.128026 0.442065-0.157955 251.00000 0.08620 0.06016 = Nd2 Cg Ca Cb 2Fo-Fc map contour 1.5s PART 2 N 3 0.130257 0.

362774 0.029087 11.00000 0.07826 0.07741 = 0.07281 0.01071 0.00458 0.00521 AFIX 13 HB 2 0.164677 0.362215 0.017819 11.00000-1.20000 AFIX 0 CG1 1 0.183985 0.380865 0.031791 11.00000 0.08496 0.

26 マルチコンフォメーションと水素原子の取り扱い Asn198 RESI A:198 ASN Od1 PART 1 PART 1 PART 2 N = Nd2 Cg Ca Cb 2Fo-Fc map contour 1.5s PART 2 N = PART 0 FVAR の 25 番目の値 PART 1 の占有率 :0.521 PART 2 の占有率 :0.479 Val68 Cg2 Cg1 Cb Ca Fo-Fc H omit map contour 2.5s RESI A:68 VAL CB = AFIX 13 HB AFIX 0 CG = AFIX 33 HG1A HG1B HG1C AFIX 0 26

27 完全行列最小二乗法による esd の計算最小二乗法のコマンド CGLS L.S. L.S. 1 ( サイクル数 1 とする ) 変数をシフトさせず標準偏差のみを見積もる DAMP 0 0 ( パラメータ数多い時は BLOC 1 も使用し原子座標の標準偏差のみ計算 ) 全ての制約を除外する ( 結合距離温度因子に関するものなど全て ) 結合距離 ( 角度 ) とesdの出力コマンド ( 水素原子を含む結合については BOND $H) SHELXLの実行 shelxl refine-ls -b t4 BOND -b オプション : パラメータ数 (default は 3000) Param thred < b5r 計算時間約 1 時間 20 分 (Core i7-4930k パラメータ数: thred 数 :4) 27

28 おわりに高分解能データ (d min ~ 1 Å) が得られた際には SHELXL を用いた精密化を検討してみてはいかがでしょうか謝辞データ収集でお世話になった PF ビームラインスタッフの皆さまに感謝いたします 28

ydkr_manu_v12.PDF

ydkr_manu_v12.PDF Yadokari-XG Ver. 1.2. 2002/02/22 4 Yadokari-XG c: ydkr elements.dat sg.dat sym.dat sym2.dat ydkr.exe ydkr.ini ysg.exe ortep3.exeortep-3 ntserver.exeortep-3 grfont.datortep-3 Yadokari-XG SHELX PATH c: windows