1. （1）フラグメント分子軌道（FMO）法^（注1）を用いて、タンパク質と固体表面の高精度相互作用モデリングに世界で初めて成功しました。
2. （2）従来は古典的力場に基づく手法でしか扱うことのできなかった、ナノ－バイオ複合系について、全電子（量子化学）計算を初めて実現しました。ペプチド中のアミノ酸残基やナノ結晶の部分構造との相互作用を電子状態レベルで高精度に解析することができます。
3. （3）本手法はナノ－バイオ境界の分子設計に広く拡張させることができ、生体親和性の高いインプラントの開発や、バイオセンサーやナノイメージングなどの医療工学、さらにはバイオミネラリゼーションによる機能性マテリアルの開発に役立つと期待できます。

1999年に北浦和夫氏（現 神戸大学 特命教授）らによって提案されたフラグメント分子軌道（FMO）法 ^[1]は、大規模分子系を量子論的に扱うことが可能であり、これまで生化学や生物物理化学、とりわけ創薬の分野で主に用いられてきました。これは、タンパク質のアミノ酸残基と薬剤分子との相互作用をフラグメント間相互作用エネルギー（IFIE）を指標にして詳細に解析することが直截に可能であり、薬剤分子とタンパク質の結合様態の理解、さらには薬剤の骨格・置換基の改変による最適化に好適であるためです^[2]。

一方で、固体と生体分子の相互作用は、インプラントの表面改質、微小粒子による薬剤投与（ドラッグデリバリ）などの医療工学、あるいは生体結晶析出（バイオミネラリゼーション）などの医療工学、生物工学などいわゆるナノ－バイオの境界領域の問題として近年重要視されています。特に、特定組成の固体面を認識する人造ペプチドは世界中のグループが研究開発を進めており、日本ではがん研究所の芝清隆氏らが先端を走っています。

こうした中、本研究グループでは、独自開発のABINIT-MPプログラムを拡張して4体までのフラグメント展開を行う「FMO4法」^（注2）を実装し、4体補正によって計算精度が飛躍的に向上することを実証しました^[3]。FMO4によって創薬分野においては官能基単位の相互作用解析が可能となり、高い空間解像度と定量性を両立させた次世代のIFIE解析として製薬企業から注目を集めています^[4]。今年度は、ABINIT-MPの高速化（日本電気株式会社と共同）や解析手法の洗練（神戸大学の田中成典教授、国立医薬品食品衛生研究所の中野達也室長らと共同）を進めてきました。そして今回、結晶系の高精度なモデリング手法を新たに開発し、ナノ－バイオ境界系へのアプローチが可能となりました。

FMO以外にも、大規模系の定番手法としては古典分子動力学（MD）もしくは古典／量子論混成（QM/MM）モデリングがありますが、古典MD法では電荷移動や化学反応などの精密な現象を扱うことができず、またQM/MMにおいては今回の系のようにペプチドと固体側の広範囲に渡る界面を量子論（QM）で扱わなければならない場合、実際の計算はきわめて難しくなります。また、固体物理で広く利用されている密度汎関数計算（CPMDなど）も、系の非周期性や分散力の欠如の点から適用は現時点では困難です。

ABINIT-MPおよびBioStation Viewerの詳細については、こちらをご覧ください。

今回、研究グループはFMO4法をシリカ（酸化ケイ素）表面への6残基のアミノ酸から成るペプチド（微小タンパク質）の結合状態のモデリングに応用することに世界で初めて成功しました。Si－O結合の3次元的ネットワークを持つ固体シリカの扱い^（注3）はFMO4法による高精度化とシリカ結晶の3次元的な分割手法の開発で可能になったものです。

対象とした6残基ペプチドは、芝清隆氏の研究グループが開発した、チタン、銀、およびケイ素の酸化物表面に対してのみ特異的に結合する12残基のペプチドの前半部分、表面固着能を決しているとされる6残基（Arg1-Lys2-Leu3-Pro4-Asp5-Ala6：RKLPDAと略記）^[5,6]です。このRKLPDA片をシリコン原子を257個含むシリカのクラスターモデルに組み合わせた複合系（水和も考慮）に対し、分散力を取り込める2次の摂動論（MP2）レベルで大規模な計算を行いました。実行にはノード当たり12個コア、それを168ノード持つFOCUSスパコンを用い、特に時間のかかる4体フラグメントの処理時間の短縮のためにコレスキー分解（CDAM-MP2法[7]）^（注4）を使っています。相互作用エネルギー解析では、フラグメント間の遮蔽効果を統計的に考慮した新しい方法：SCIFIE^{[8]（注5）}も採用しています。

図はIFIEによる相互作用を可視化したものですが、シリカ表面への結合ではRKLPDAの中で荷電したArg1、Lys2、Asp5の3つのアミノ酸残基がとりわけ重要であることが確認されました。こうした情報は、ペプチドの改変や最適化に対して計算の立場から有益な指針を提供するものと言えます。一方、シリカ側もペプチドの結合によって表面だけでなく奥も分極されており、固体側をより簡便な小型のクラスターモデルで近似することはできないことも分かりました。

図1　水和したシリカ－ペプチド複合系の相互作用の可視化

左側はペプチド側、右側はシリカ側に注目して合算している。赤は結合による安定化、青は不安定化を示す。

研究グループは、今回成功したFMO4計算によるペプチドとシリカ表面のモデル計算を発展させることを計画しており、別種のペプチドを組み合わせて相互作用のパターンの差異を見ることは興味深いと思われます。これは、表面に強く結合するアミノ酸残基の特定に役立つのはもちろん、むしろ吸着しない・させないペプチドを設計するのにも使えます。また、固体結晶の欠損やドープなどのモデリングやシリカ以外の固体表面でアルミナやゼオライトなどを対象にすることも考えています。さらに、酸化物に留まらず各種のエンジニアリングポリマーやダイヤモンド薄膜などを扱うことも視野に入っています。計算資源としてはHPCIを用いる予定です。こうした研究は、計算・シミュレーションを援用する「ものづくり」の文脈に合致すると言えます。

また、本計算手法は、FMO法の地球科学分野への応用に道を開くものです。シリカは、石英などの鉱物の基本構成要素であり、水和条件下での鉱物表面への各種イオンの吸着・脱着は重要な化学プロセスですが、電子状態計算による理解はこれまではほとんどなされていませんでした。このFMO4計算の技術を発展させることにより、熱力学的なパラメータを含めた詳細な描像が得られるようになると期待されます。

いずれにせよ、今回の成果はこれまで主に生化学・生物物理化学・薬学分野に限られてきたFMO法の適用範囲を拡大するブレークスルーと位置づけられます。

図2　シリカ表面とイオンの相互作用の可視化

左側はナトリウムイオン右側は塩素イオン。赤は結合による安定化、青は不安定化を示す。

Y. Okiyama, T. Tsukamoto, C. Watanabe, K. Fukuzawa, S. Tanaka, Y. Mochizuki, “Modeling of peptide-silica interaction based on four-body corrected fragment molecular orbital （FMO4） calculations”, Chemical Physics Letters 566 （2013） 25-31.

1. 注1フラグメント分子軌道（FMO：Fragment Molecular Orbital）法
   1999年に北浦和夫氏（現 神戸大学 特命教授）によって提案された生体分子系に対する効率的な計算手法[1]。巨大系を比較的小さなフラグメント（アミノ酸残基など）に分割し、各フラグメントのモノマーとダイマーの分子軌道計算を並列処理することにより、全系の電子状態をこれまでの手法よりはるかに短時間に高精度で解析することができる近似計算法。計算からは、フラグメント間の相互作用エネルギー（IFIE（Inter Fragment Interaction Energy）；別グループではPIE（Pair Interaction Energy）とも呼ぶ）が得られるため、創薬分野での応用計算が盛んに行われている^[2]。
2. 注2FMO4法
   フラグメントの展開を4体（テトラマー）まで拡張することにより、リガンドの複数分割とアミノ酸残基の主鎖・側鎖分割を伴う高解像度の解析を可能とするアプローチ^[3,4]。従来と異なり、FMO4計算ではダイヤモンドやシリコンなどの3次元の結合ネットワークを持つ固体をも扱うことが出来る。
3. 注3Si-O結合の3次元的ネットワークを持つ固体シリカの扱い
   3次元的に配置したSi-O結合をもつシリカ結晶の切断では、結晶のユニットセルを基本構造としたフラグメントに分割する。形式電荷やフラグメントサイズが均一になるように配置している。
4. 注4コレスキー分解（CDAM）による高速化
   系の4乗に比例する4中心の2電子積分をコレスキー分解によって近似的に3中心で評価する技法で、高速でありながら化学的な精度は担保される^[7]。FMO4のMP2計算では、コストは半分程度に抑えられる。
5. 注5SCIFIE（Statistically Corrected Inter-Fragment Interaction Energy）
   SCIFIEは、古典統計力学の多体論的な手法を用いてIFIEからフラグメント間の実効ポテンシャルを求める方法である。イオン性が強い系ではIFIEが過大評価される傾向にあるが、SCIFIEでは統計的な補正手法により遠方のフラグメントで期待される減衰効果を含めた値を得ることが出来る^[8]。

1. [1]K. Kitaura, E. Ikeo et al., Chem. Phys. Lett. 313 （1999） 701
2. [2]The Fragment Molecular Orbital Method: Practical Applications to Large Molecular Systems, CRC Press, Boca Raton, 2009.
3. [3]T. Nakano, Y. Mochizuki et al., Chem. Phys. Lett. 523 （2012） 128.
4. [4]http://www.iis.u-tokyo.ac.jp/publication/topics/2012/20120326press1.pdf （PDF/647KB）
5. [5]K. Sano, K. Shiba, J. Am. Chem. Soc. 125 （2003） 14234
6. [6]K. Sano, H. Sasaki et al., Langmuir 21 （2005） 3090
7. [7]Y. Okiyama, T. Nakano et al., Chem. Phys. Lett. 494 （2010） 84
8. [8]S. Tanaka, C. Watanabe, et al., Chem. Phys. Lett. 556 （2013） 272.