生化学 SEIKAGAKU
Journal of Japanese Biochemical Society

タンパク質の機能を理解するには，その立体構造を可視化することが不可欠である．100万分の1ミリメートル程度しかないタンパク質の形を知るために，分子や原子レベルの極小の世界を分析できる技術，X線結晶構造解析が用いられてきた．この技術は，タンパク質の結晶を回転させながらX線を照射し，結晶が反射した光の模様（X線回折画像）を記録し，そこからコンピュータ上での解析を通じて，結晶を構成する原子の形（電子密度）を描き出す．高解像度で明瞭に構造を決定するには，従来，大型（＞200 µm）で良質の結晶を作ることが肝要であった．一方，より小さな結晶でも明瞭なX線回折画像を得るために，強力なX線を生み出す放射光施設が建設されてきた．

X線自由電子レーザー（XFEL）は，超高輝度・フェムト秒パルス・高空間コヒーレンスという特徴を有する新世代の光である．米国で世界初のXFEL施設LCLSが2009年に建設されたのに次いで，わが国では2011年，国家基幹技術としてSACLAが兵庫県播磨科学公園都市に建設された．SACLAが生み出す光はSPring-8など従来のシンクロトロン放射光施設が出す光の10億倍もの明るさを持つ．本稿では，XFELを利用して数十ミクロンもない多数の微結晶からタンパク質構造を解明する新しい実験手法「連続フェムト秒構造解析（serial femtosecond X-ray crystallography：SFX）」について実例を通して紹介する．

SFXは，インジェクターから噴出させた多数の微結晶に対し，XFELを連続的に照射してX線回折画像を収集して構造解析する手法である（図1）．SACLAでは，XFELが毎秒30発，光の弾丸として機関銃のように微結晶の流れに向けて発射される．利用できる光の波長範囲は0.62～2.76 Åである．1発のXFELはわずか10フェムト秒（100兆分の1秒）以下のパルス幅を持つ．微結晶にパルス光がヒットすると，その衝撃で結晶は一瞬で崩壊してしまう．しかし，結晶の崩壊過程よりも早くX線回折画像を記録できるため，構造解析に支障はない．数万から数十万個の微結晶についてのX線回折画像を数時間かけて測定し，常温下でダメージフリーのタンパク質の真の構造を決定できる．SFXの登場により，従来の放射光X線構造解析で問題となっていた，測定時の結晶凍結や放射線損傷が引き起こすタンパク質構造の変位や歪みを克服することが可能となった．

図1 連続フェムト秒構造解析

連続フェムト秒構造解析の概念図（左図）．わが国では兵庫県播磨科学公園都市にXFEL施設SACLAが大型放射光施設SPring-8に隣接して設置されている（右図）．

ここでは，地球上の窒素循環において重要な役割を担う銅含有亜硝酸還元酵素（copper-containing nitrite reductase：CuNiR）の解析事例を紹介する^{1, 2）}．本酵素は活性中心に銅原子を有し，亜硝酸イオンを一酸化窒素に一電子還元する反応（NO₂⁻＋2H^＋＋e⁻→NO＋H₂O）を触媒する．一般に，従来のシンクロトロン放射光を利用した金属タンパク質のX線結晶構造解析では，結晶に照射されたX線が水和電子を発生させ，活性中心の金属原子を還元してしまう現象（X線光還元）が避けられなかった．CuNiRも同様で，X線回折画像の収集時に銅原子のX線光還元が速やかに起こることで触媒反応が進行し，基質NO₂⁻から生成物NOや副生成物が生じてしまう．そのため，反応開始前の真の酸化型構造を観察することが不可能であり，反応機構の全貌は未解明のままであった．

そこで筆者らは，SACLAにてSFX実験を行い，光還元のない真の酸化型構造を常温（20°C）下で捕捉することを試みた．並行して，SPring-8の放射光X線を利用し，CuNiRの触媒反応をX線光還元により誘発させた構造を極低温（−170°C）と常温下で決定して，SACLAの構造との比較対照とした．従来から知られているように，X線照射前の緑色をした酸化型の結晶は，SPring-8でのデータ測定中のX線光還元に伴い無色に変化してゆく様子が観察された．興味深いことに，SPring-8で解いた構造は極低温と常温下のいずれの条件においても，基質NO₂⁻は銅原子に対してフェイスオン（face-on）型に近い配位様式をとっていたのに対し，SACLAの構造では直立（vertical）型で配位していることが明らかとなった（図2）．過去に報告された放射光で解かれたNO₂⁻複合体構造を調べると，すべてフェイスオン型に近い配位様式であった．このことから，反応前の酸化状態で直立型に配位したNO₂⁻は，反応中の銅原子の還元に依存して配位様式が倒れこむように構造変化すると考えられた．

図2 異なる条件で測定したCuNiRの亜硝酸イオン（NO₂⁻）複合体の構造比較

活性中心の銅原子への亜硝酸イオンの配位様式が測定条件ごとに大きく異なっていた．この違いは，銅原子がSPring-8のX線による光還元を受けることで酵素反応が進行することにより起こる．図は文献1）より改変して転載．

また，本稿では詳細な記載を省くが，銅原子の酸化還元状態に依存し，銅原子に配位するHis残基の側鎖が大きく回転する現象も見いだされた．これらの新しい知見から筆者らは，CuNiRの触媒反応機構におけるプロトン共役型電子移動の新しいモデルを提唱している．本研究成果は，新時代の光XFELと従来の放射光X線による構造解析手法を併用することにより，未解明であったタンパク質の作用機作を初めて発見できた事例として重要である．

筆者らがSFXに着手した当時，ほとんどのSFX構造は既知構造をモデルとする分子置換法によって解析されていた．SFXを従来の放射光施設では解けない新規標的に適用し，立体構造を常温かつダメージフリーの状態で明らかにするには，SFXでの，分子置換法によらない構造決定，すなわち実験的位相決定法の確立が不可欠であるが，その特殊性のため近年まで実証例は少なかった．

実験的位相決定には，同型置換法と異常散乱法の2種類と，双方を組み合わせた方法がある．同型置換法の場合は，ネイティブ結晶と重原子標識した誘導体結晶の反射強度の違いを利用する．異常散乱法の場合は，バイフット（Bijvoet）対と呼ばれる反射どうしの強度の差を利用する．いずれにせよ，微小な違いを正確に測定することが必要である．しかし，SFXでは数万個以上の独立した結晶から，パルス状のX線によってデータを収集するので，観測値のばらつきは非常に大きい．その原因としては，構造の非同型性・重原子占有率の違い・結晶形態やサイズの違いといった結晶に由来するものや，ビーム強度やスペクトルの違いといった光源に由来するものがあげられる．また，XFELならではの問題点として，すべての反射が部分反射であることによる部分度（partiality）の影響もある．モンテ・カルロ法によるSFXデータセットのマージとは，反射を何度も測定して多重度を高めることによって，これらの違いを平均化することに他ならない．

理論上，モンテ・カルロ法のデータ精度は多重度の平方根に比例して改善する．つまり，バイフット対の相対振幅差であるバイフット比が1％の系で位相決定をするには，バイフット比10％の系の100倍のデータが必要である．ただし，収束のオーダーが理論どおりだとしても，比例係数は実験データの分散に依存するため，現実的な枚数のX線回折画像では位相決定できない懸念がある．また，これは正確性ではなく精度の議論なので，データ処理に系統誤差がある場合，多重度をいくら増やしても位相決定できない可能性すら存在する．

2014年までにSFXで実験的に位相が決定された例は1件のみであった^3）．ここでは，ニワトリ卵白リゾチームのガドリニウム誘導体から単波長異常散乱（SAD）法で位相が決定された．この系は10％を超える大きいバイフット比を持つため，異常散乱シグナルがより小さい現実的な系で位相決定可能かは未知であった．特に硫黄によるSAD（S-SAD）は，重原子誘導体化を必要とせず，システインやメチオニンに含まれる硫黄原子から位相決定できる手法であるが，異常散乱シグナルが非常に小さいため困難が予想された．

我々はリゾチームのS-SADに挑戦し，指数がついた13万枚分のX線回折画像から位相決定に成功した^4）．この系は，非対称単位あたり10個の硫黄原子と1個の溶媒由来の塩化物イオンを含み，理論上のバイフット比は約1.6％であった．部分構造探索にはSHELXD，位相計算・位相改善にはSHELXEという，通常のX線結晶構造解析でも頻用されるプログラムを用いたが，分解能上限や溶媒比率などのパラメータを網羅的に振った大量の計算を実行することが成功に不可欠であった．また，人間にはほとんど解釈不能な電子密度マップでも，SHELXEでの主鎖自動トレースと位相改善を25回反復することで，自動モデリング・プログラムBuccaneerによる側鎖のあてはめが可能なマップに改良できることを驚きとともに知った．

上述したCuNiRの結晶は，非対称単位に六つの銅原子を含み，理論上のバイフット比は約1.7％となる．リゾチームのS-SADと同様にパラメータの網羅的探索を実行したところ，Cu-SADによる位相決定に成功した．SFXのために開発されたデータ処理プログラムCrystFELバージョン0.5.4を用いて，指数がついた15.6万枚の回折画像からの位相決定を論文報告した^1）．論文では最低必要枚数の検討を行わなかったが，実際には指数がついた回折画像10万枚で十分であった．

S-SADについては，のちにソーマチン^5）や膜タンパク質であるA2a受容体^6）でも成功が報告された．それぞれの理論上のバイフット比は約2.1％と1.9％であり，指数がついた回折画像がそれぞれ12.5万枚，50万枚必要であった．これらの報告により，バイフット比が2％未満の系でもSFXでSAD法による位相決定ができることが確立された．しかし，いずれの場合でも膨大な枚数のデータが必要であり，サンプルやビームタイムの消耗が著しい．必要枚数を減らすには，シグナルを強めるのと，データ処理の質を高めるという二つの方向で取り組める．

シグナルを強めるアプローチとして，我々は膜タンパク質の重原子誘導体の新しい調製法を開発した（図3）^7）．膜タンパク質の結晶構造では，分子表面に界面活性剤や脂質分子を認めることが多々ある．そこで，ヨウ素を三つ含む重原子試薬I3C，通称magic triangleをカプリル酸とアミド結合させた化合物HAD13aを合成した．これをバイセル中で得られたバクテリオロドプシン微結晶と混合するだけで結晶の重原子標識が可能であり，7 keVでデータ収集した2.3万枚の指数がついた回折画像から，SAD法で位相決定することに成功した．非対称単位あたり2分子のHAD13a，つまり六つのヨウ素が結合しており，10.5％という高いバイフット比を得られたのが，少ない枚数での成功につながった．仮にHAD13aを使わず，S-SADをするとすると，バイフット比はわずか0.9％である．次に，HAD13aを含まないネイティブ結晶のデータと組み合わせて，異常散乱効果を伴う単一重原子同型置換（SIRAS）法を行ったところ，ネイティブ3000枚・誘導体4000枚，合計わずか7000枚の指数がついた回折画像から位相決定ができた．さらに，測定したすべての回折画像を用いてSIRASを行うと，分解能3.3 ÅというSFXでの位相決定における最低分解能記録で構造解析することもできた．これは同型置換シグナルと異常散乱シグナルが相補的な位相情報を提供するためと考えられ，効率のよい位相決定には，ネイティブ結晶の併用が効果的であることを確認できた．なお，HAD13aは，人工的な脂質二重膜である脂質キュービック相（LCP）の中で結晶化したGタンパク質共役受容体にも結合でき，異常散乱シグナルを与えることを確認しており，今後さまざまな膜タンパク質に応用可能と期待している．このほか我々は，水溶性タンパク質で水銀^8）やプラセオジム^9）などの重原子を用いた位相決定にも成功している．

図3 膜タンパク質を効率的に重原子標識してSFXで構造決定する技術を開発

下段の（A）と（B）は，それぞれ，膜タンパク質の全体構造とHAD13a結合部位の拡大図を示している．図は文献7）より改変して転載．

SFXのデータ処理法には，まだ改善の余地がある．CrystFEL付属のpartialatorというソフトでスケーリングできるようになりつつあるが，温度因子（B factor）が発散するデータセットも散見され，万全とはいいがたい．結晶品質にバラつきがある場合の処理も問題である．一定の分解能に達しない結晶を棄却したり（min-resオプションを使用），結晶ごとにマージに使う分解能を変えたり（push-resオプション）するとデータの質がよくなることがある．一番の課題は部分度である．post-refinementによる部分度の補正はいくつか報告されており^{10, 11）}，BinABタンパク質の，複数の誘導体と異常散乱効果を用いたMIRAS法による位相決定^12）でも用いられたが，筆者らの経験では，データセットによって効果は異なり，かえって悪化することも多い．我々のCuNiRの系では，CrystFELの（執筆時点での）最新バージョン0.6.2を使うことで，必要な枚数を10万枚から8万枚に減らすことができている．また，post-refinementによってさらに枚数を削減することに成功しており，論文報告を準備中である．

SFXを応用すると，タンパク質が働く際の構造変化をムービーのように捉えることができる（時分割解析）．これまでに我々は，バクテリオロドプシンや光化学系などの膜タンパク質でその技術を実証することに成功している^{13, 14）}．時分割解析には，上述の実験的位相決定において経験したSFXで微小な違いを正確に検出する手法の開発が活かされている．構造生物学は，タンパク質の静止構造から動的構造を解明するダイナミクス研究の新時代に入ったといえよう．

なお，本稿で紹介した研究は，文部科学省X線自由電子レーザー重点戦略研究課題「創薬ターゲット蛋白質の迅速構造解析法の開発」，JST-ERATO「村田脂質活性構造プロジェクト」，および，JSPS科研費15K18487の支援を受けて実施した．計算解析はSACLA HPCシステムおよびミニ京スーパーコンピュータシステムを利用して実行した．関係各位に深く感謝したい．