Simulation of Protein Structure using a Coarse-Grained Potential incorporating the Backbone Dihedral Interactions

ハリケーンの中で、たった一粒の砂が踊る様子を見ようとしている自分を想像してみてください。次に、ハリケーンが吹き荒れる中で、砂粒のビーチ全体が渦巻き、衝突し、共に城を築き上げている様子を見ようとしている自分を想像してください。これが、私たちの体を動かしている微小で複雑な機械である「タンパク質」を研究する際に、科学者が直面する課題です。タンパク質は、アミノ酸の長い鎖からできており、それが複雑な3D形状へと折り畳まれます。タンパク質がどのように動き、相互作用するかを見るために、科学者は通常、すべての原子を個別のビリヤードの球のように扱う「全原子（all-atom）」シミュレーションを使用します。しかし、これは、砂粒のビーチ全体を追跡するために、一つ一つの砂粒を個別に記録しようとするようなものです。それには膨大なコンピュータパワーが必要となるため、ほんの一瞬の間、ごくわずかな砂の塊しか観察することができません。タンパク質が数千個集まって疾患の原因となったり、細胞壁を構築したりするような、大きな全体像を見るためには、科学者は視点を外側に広げる必要があります。彼らは「粗視化（coarse-grained）」モデルを使用します。これは、多くの原子を一つの「ビーズ」としてグループ化するもので、例えるなら、一掴みの砂を一つの大理石にまとめるようなものです。これにより速度は劇的に向上しますが、落とし穴があります。簡略化を急ぐあまり、モデルはしばしばタンパク質が折り畳まれるための「秘密のレシピ」を失ってしまうのです。それは、タンパク質の形を与える特定の角度やねじれを忘れ、マスターシェフを、単に材料を混ぜる方法だけを知っている人に変えてしまうようなものです。

ここから、この論文の物語が始まります。研究者の Kanika Kole、Abhik Ghosh Moulick、そして Jaydeb Chakrabarti は、「スマート」な粗視化モデルを構築したいと考えました。彼らは問いかけました。「ビーチ全体のダンスを観察できるほど速く、かつ、秘密の折り畳み角度を記憶できるほど賢い、簡略化されたタンパク質モデルを作れるだろうか？」と。彼らは、タンパク質鎖の各アミノ酸を単一のビーズとして扱う新しい手法を開発しましたが、それぞれのビーズに5つの特別な「自由度」を与えました。これらのビーズを、単なる静的な大理石としてではなく、揺れたりねじれたりできる柔軟な関節として考えてください。決定的なのは、彼らがこれらのビーズがどのように動くべきかを推測したのではなく、GB3と呼ばれるタンパク質の超詳細なスローモーション映画からルールを「盗んだ」ことです。その詳細な映画におけるあらゆる原子の動きを分析することで、ねじれや曲げに対する正確なエネルギーコストを抽出し、それらのルールを彼らの高速で簡略化されたモデルにプログラムしました。彼らはこの新しい「スマートビーズ」モデルを、小さな整ったグロビュールから、固定された形を持たない乱雑でフラフラとした鎖まで、さまざまなタンパク質を用いてテストしました。その結果、彼らのモデルは、実在の結晶構造や詳細なシミュレーションに見られる複雑な形状や折り畳みパターンを再現することに成功し、しかも遥かに高速で動作することを見出しました。さらに素晴らしいことに、あるタンパク質から学んだルールが、全く異なるタンパク質の形状を予測するために使用できることを示しました。これは、彼らのモデルが流暢に話すことができる、タンパク質折り畳みの普遍的な「言語」が存在することを示唆しています。

スマートビーズの物語

生物学の世界において、タンパク質は究極の変身者です。それらは（アミノ酸という）ビルディングブロックの長い鎖であり、その役割を果たすために特定の3次元構造へとねじれ、回転します。時には、折り紙の鶴のように硬く明確な形を持つこともあります。またある時は、「本質的に無秩序（intrinsically disordered）」、つまり、何かに捕まえた時に初めて形を成す、調理されたスパゲッティのボウルのように、フラフラとして混沌としていることもあります。これらの鎖がどのように動き、折り畳まれ、集まるのかを理解することは、生命や疾患を理解する上で不可欠ですが、それをシミュレートすることは非常に困難です。

この論文の著者たちは、コンピュータサイエンスにおける古典的な問題、すなわち「速度と詳細度のトレードオフ」に取り組みました。タンパク質の動きをリアルタイムで見たいのであれば、高速で簡略化されたモデルが必要です。しかし、簡略化しすぎると、「二面角（dihedral angles）」、つまりタンパク質の鎖が曲がる特定の角度を失ってしまいます。これらの角度こそが、タンパク質が頑丈なヘリックス（螺旋）になるか、平らなシートになるかを決定する「秘伝のソース」なのです。これらがなければ、モデルは単なる、本来あるべき姿の記憶を持たない、フラフラした紐になってしまいます。

チームの解決策は、「ハイブリッド」モデルを作成することでした。タンパク質の鎖を、一連のビーズの紐だと想像してください。彼らのモデルでは、各ビーズは一連のアミノ酸全体の中央を表します。しかし、これらは普通のビーズではありません。それらは、どのように曲がるのが好きかという記憶を備えた「スマートビーズ」なのです。研究者たちは、各ビーズに追跡すべき2つの追加変数、すなわち主鎖の二面角である $\phi$（ファイ）と $\psi$（プサイ）を与えました。これらはタンパク質の主鎖のねじれを定義する角度です。

これらのビーズにどのように振る舞うべきかを教えるために、著者たちは単にルールを捏造したわけではありません。彼らは、GB3と呼ばれる小さなタンパク質の高精細な「全原子」シミュレーションを観察しました。彼らは、原子が時間の経過とともにどのように動くかを観察し、あらゆる可能なねじれや回転における「自由エネルギー」を計算しました。これは、谷がタンパク質の快適で安定した位置であり、丘が不快で高エネルギーな位置であるような、風景を描くようなものです。そして、彼らはこのマップを、一連のルールへと翻訳しました。

彼らはまた、ビーズが水とどのように相互作用するかについても確実にしました。彼らはビーズを2つのタイプ、「親水性（solvophilic）」と「疎水性（solvophobic）」に分類しました。水を好むビーズは水分子を引き付けるようにプログラムされ、水を嫌うビーズは水分子を遠ざけるようにプログラムされました。これは、実際のタンパク質がどのように折り畳まれるかを模倣しています。つまり、水を嫌う部分はコアの中に隠れ、水を好む部分は外側に留まるということです。

大いなるテスト：小さな鶴からフラフラのスパゲッティまで

モデルを構築した後、チームはそれをテストにかけました。まず、彼らはGB3から学んだルールを使用して、GB3自体をシミュレートしました。結果はどうだったでしょうか？ モデルはタンパク質の形状をうまく再現しました。彼らが「スマートビーズ」シミュレーションを元の高精細な映画や実際の結晶構造と比較したところ、形状は驚くほど一致しました。タンパク質の構造要素（ヘリックスやシートなど）の約80%が、簡略化されたモデルによって正しく特定されました。

しかし、本当の魔法は、彼らがGB3のルールを「他の」タンパク質に適用しようとした時に起こりました。彼らは同じルールセットを取り、以下を含むさまざまなタンパク質に適用しました：

• Homeodomain: 遺伝子を制御する小さなタンパク質。
• Ubiquitin: 他のタンパク質をリサイクル用にタグ付けするタンパク質。
• Chignolin, BBL, BBA: 折り畳みのためのテストケースとしてよく使われる小さなタンパク質。
• Di-ubiquitin および Adenylate Kinase: 複雑なタスクを実行する、2つのパーツからなる大きなタンパク質。

驚くべきことに、モデルはそれらのほとんどに対して機能しました。これらのタンパク質はGB3とは異なりますが、「スマートビーズ」はそれらを実生活の対応物と密接に類似した形状へと折り畳むことができました。例えば、Ubiquitinの場合、モデルは結晶構造に対して約67%の精度で正しい二次構造を捉え、他の詳細なシミュレーションと比較した場合にはさらに優れた92%を達成しました。これは、タンパク質折り畳みの基本的な「文法」――主鎖がどのように曲がり、水とどのように相互作用するか――が、驚くほど普遍的であることを示唆しています。

チームはまた、「本質的に無秩序なタンパク質（IDP）」、具体的には $\alpha$-シヌクレインと $\lambda$N についてもモデルをテストしました。これらは、固定された形を持たない「フラフラのスパゲッティ」のようなタンパク質です。シミュレーションを実行したところ、モデルはこれらのタンパク質が現実世界と同じように、主に展開され、混沌とした状態のままであることを正しく予測しました。興味深いことに、これらの無秩序なタンパク質における「弾性定数（主鎖がいかに硬いか、あるいは柔軟か）」は、構造を持つタンパク質よりもはるかに低いことがわかりました。これは、「フラフラであること」がこれらのタンパク質の特徴であり、彼らのモデルがそれを完璧に捉えたことを示しています。

なぜこれが重要なのか（そして、何ではないのか）

著者たちは、彼らの新しいモデルを、Martini や SIRAH といった他の一般的な簡略化モデルと比較しました。彼らは、Martini がタンパク質を硬すぎたり秩序立てすぎたりする傾向があり、SIRAH がタンパク質をフラフラさせたり無秩序にしたりする傾向がある一方で、彼らの新しいモデルが完璧なバランスを実現していることを見出しました。彼らのモデルは、実験的な結晶構造と高速な全原子シミュレーションの両方に高精度で一致し、タンパク質の詳細な構造的特徴を再現することができました。

しかし、論文は、このモデルが「できないこと」についても注意深く述べています。現在のモデルはタンパク質の主鎖に限定されており、側鎖（主鎖から突き出ている「腕」の部分）や水の双極子モーメントを明示的に含んでいません。これは、タンパク質がどのように折り畳まれ、集まるかを理解するには優れていますが、薬物分子がタンパク質の特定のポケットにどのように結合するかを正確に予測するには、詳細さが足りない可能性があることを意味します。また、ルールは特定のタンパク質のシミュレーションから導出されたものであり、それらが他の多くのタンパク質に対しても機能するように見えますが、存在する「あらゆる」タンパク質に対してこの「転移性」を証明することは、この単一の研究の範囲を超えた、極めて困難な課題であることを著者らは認めています。

結局のところ、この論文は強力な新しいツールを提供しています。原子の複雑なダンスを記憶するように簡略化されたビーズを教えることで、研究者たちは、大量のタンパク質が長期間にわたって動き、相互作用することをシミュレートする方法を作り出しました。これは、パーキンソン病のような疾患で見られる塊がどのようにタンパク質が集積するか、あるいは生物学的材料がどのように自己組み立てを行うかといった、以前は観察することが不可能であった現象を研究するための扉を開くものです。これは、単一の砂粒だけでなく、ビーチ全体のダンスを見ようとする試みなのです。