SAR News No.17 (Oct. 2009)

構造活性相関部会・ニュースレター<1 October, 2009>

SAR News No.17「目次」///// Perspective/Retrospective /////
分子動力学による結合自由エネルギー計算

藤谷

秀章 ･･･２

浩明、広野修一 ･･･９///// Cutting Edge /////
MM-PBSA 法を用いた精密結合自由エネルギー計算
合田

QM/MM 分子動力学シミュレーションの生体高分子への応用

米澤康滋 ･･･１５///// Activities /////＜報告＞
構造活性フォーラム 2009「化合物と標的・非標的タンパク質との相互作用－創薬標的の同定と分子設計－」開催報告
SAR Promotion Award 2009 年度受賞者

竹田－志鷹

真由子 ･･･２１

米田

照代 ･･･２２＜会告＞第 37 回構造活性相関シンポジウム･･･２３SAR News No.17 (Oct. 2009)///// Perspective/Retrospective /////

分子動力学による結合自由エネルギー計算（株）富士通研究所

1.

藤谷

秀章

はじめに
生体機能はそれを構成する種々の生体高分子が有機的に機能することで実現されている。個々の生体高分子の機能を理解する事は生体機能を理解する出発点になるが、例
えば最も小さいタンパク質の一つである 36 個のアミノ酸残基からなる villin でも、その総原子数は 389 個にもなる。その性質や動的挙動を調べるには周りを取り囲んだ水分子も含めた取扱いが必要であり、計算化学の手法を用いて解析する事は極めて難しい問題になる。1998 年に Kollman は超大型スーパーコンピュータを用いて、villin が一本の鎖状態からタンパク質独特の立体構造を獲得するまでの畳み込み (protein folding) 現象を、
分子動力学法を用いて 1 マイクロ(10-6) 秒までをシミュレートした。しかし原子レベルの時間スケールがフェムト(10-15) 秒に対して villin が実際に畳み込まれる時間が 10 マイクロ(10-5) 秒ほどであるために、畳み込みが終了する時間には到達出来なかった[1]。
2002 年に Pande は世界中のパーソナルコンピュータにシミュレーションプログラム
をダウンロードしてもらう distributed computing の方法を用いて、villin の畳み込み計算を実行して実験値に近い畳み込みレートを求める事に世界で初めて成功した[2]。Pandeはスーパーコンピュータなどで用いられている従来の並列アルゴリズムではなく、世界中のインターネットで繋がっている数万台のコンピュータに分散させてシミュレートする独自の方法を考案して超難問のタンパク質の畳み込み計算を初めて実現した。スーパーコンピュータなどで用いられている message passing interface (MPI) などの並列化方法では CPU 台数が多くなるに従い CPU 間のデータ通信が占める割合が増加して効率的な計算を妨げる傾向があるので、例えば数万台の CPU で数万倍の効率的なシミュ
レーションを実行するためには、シミュレーションの根本アルゴリズムから検討する事が必要になる。
タンパク質とこれに特異的に結合する化合物(リガンド)の相互作用を調べる事は、生命現象の基礎的研究の大きな課題であると同時に応用の観点からも極めて重要である。
特にタンパク質と薬分子との結合自由エネルギーを精確に求める事は、創薬分野における永年の課題である。Docking simulation や molecular mechanics Poisson-Boltzmann surface
area (MM-PBSA) 法など様々な方法が試みられて来たが、実用レベルで有効な精度で結合自由エネルギーを予測する方法はなかった。我々は非平衡統計力学理論に基づいた計算方法を用いると共に、タンパク質などの生体高分子と有機化合物分子を統一的に扱う
分子力場を構築して、タンパク質とリガンドが離れた状態から結合するまでのマイクロ秒スケールの分子動力学計算を実行して結合自由エネルギーを求める方法を考案した。

-2-

SAR News No.17 (Oct. 2009)2.

薬とタンパク質
一般に薬分子は決まった標的タンパク質に結合してタンパク質本来の機能を阻害または補助して薬効を発する。図 1 は臓器移植やアトピーなどの病気治療に用いられる免疫抑制剤 FK506 がその標的とするタンパク質 FKBP (FK506-binding
protein) に結合している状態である。FKBP は 107
個のアミノ酸残基からなるタンパク質で原子数は 1666 個である。FK506 が FKBP に結合する事
で、T 細胞での脱リン酸化を阻害して免疫機能を抑制する。
タンパク質とリガンドは一般に化学結合を伴

図 1

FKBP と FK506 の結合状態

わない静電相互作用、van der Waals 相互作用、水素結合などによる極めて弱い結合力の集積で結合する。あるリガンドがタンパク質に強く結合していると、タンパク質の機能をリガンドがより良く阻害して薬効が良く現れる可能性が高くなる。この為、種々のリガンドがタンパク質に結合するエネルギーを精確に求める事が、計算機を使った薬設計(computer drug design) を実現する為の必須条件になる。この時問題になるのはタンパク質とリガンドが静的に引き合うエネルギーではなく、水に囲まれて全ての原子が熱運動をしている状態でタンパク質とリガンドを引き離す為に必要な熱力学的な自由エネル
ギーである。一般に薬分子とタンパク質の静的な引力が -80 kcal/mol から -100 kcal/mol
あっても、その結合自由エネルギーは -10 kcal/mol 程度であり、薬候補化合物間の優劣を判断するには 1 kcal/mol ほどの精度で結合自由エネルギーを求める必要がある。室温での熱揺らぎエネルギー kBT (kB: ボルツマン定数、T: 温度) が約 0.6 kcal/mol である事を考えると、1 kcal/mol の精度で結合自由エネルギーを求める事が如何に難しいかが納得出来る。
タンパク質とリガンドの結合
自由エネルギーを求める為に様々な計算方法が今迄に提案されているが、その多くが実験データに基づいた経験的な方法であり、実験データの無い未知の系に適用するのは難しい。タ
ンパク質とリガンドが水中 (細胞液) で結合した状態と、それが別々に水中に溶けている二つの異なった熱平衡状態の間の自由

図 2

-3-

FKBP リガンド

SAR News No.17 (Oct. 2009)

エネルギー差を求めるには、タンパク質の構造変化や分子の運動によるエントロピー変化を精確に取り入れる方法が必要になる。しかし薬は数十から三百くらいの原子からなる分子であり、タンパク質は数百から数十万原子からなる巨大分子である。また自由エ
ネルギーには周りを取り囲む水分子の運動が多大な影響を与えるので、これらを全て取り入れて結合エネルギーを 1 kcal/mol の精度で精確に計算するのは 20 世紀の最先端のスーパーコンピュータを用いても困難であった。
図 2 は今回結合自由エネルギー計算を行った FK506 とそのアナログである。L8、L9、
L13、FK506 の四種類のリガンドには FKBP との結合構造の X 線結晶データが Protein
Data Bank に登録されているが、他のリガンドについてはデータが無い。全てのリガン
ドはピペコリン酸エステルと α-ケトアミドを共通で持っており、X 線結晶データがあるリガンドでは、これらが図 1 で FKBP の真中に見られる疎水ポケットの底に入っている。
この事から他のリガンドも同様の結合形状を持っていると推測して結合構造の初期座標を作成した。

3.

非平衡等式と仕事量
タンパク質とリガンドの結合自由エネルギーを求めるには、リガンドが水分子やタンパク質など他の分子やイオン原子に及ぼす相互作用を消去する二重消去法を用いる。
∆Gcomplex: PL(SOL) → P(SOL) + L(GAS)(1)∆Gsol:

L(SOL) → L(GAS)(2)∆Gbind:PL(SOL) → P(SOL) + L(SOL)(3)

タンパク質とリガンドが水中で結合している状態 PL(SOL) から、リガンドが他分子と相互作用しないで空中に浮んだ状態 L(GAS) になる計算 (1) と、水中のリガンドが水分子と相互作用しなくなる計算 (2) を行い、二つの計算から得られた自由エネルギーの差を取って結合自由エネルギー ∆Gbind = ∆Gcomplex − ∆Gsol を得る (3)。もしリガンドがタン
パク質と全く相互作用しなければ ∆Gcomplex は ∆Gsol と同じ値になり結合エネルギーはゼロになるので、二重消去法を用いるとゼロ点から計った絶対的な結合自由エ
ネルギーが得られる。
1997 年に Jarzynski は熱平衡状態間の非平衡仕事量
W の指数関数平均が二つの状態間の自由エネルギー差になるという基礎的な等式∆G = −kBTln <e−W / kBT>av を証明した[3]。その後、幾つ

図 3

FKBP と L14 の相互作用ポテンシァル-4-SAR News No.17 (Oct. 2009)

かの証明方法が得られているが、分子動力学でよく用いられる能勢-Hoover の温度制御を用いた場合にも成立する事が証明されている。しかし数値計算でこの等式が威力を発揮する為には注意するべき事が二点ある。一つは十分多くの非平衡仕事量をサンプルする事で、もう一つは計算に使用する初期状態を十分に熱平衡化する事である。タンパク質の様な自由度が大きな系の熱平衡状態を得るには、これまで行われていたピコ秒オー
ダーの分子動力学による平衡化計算では不十分で、遥かに長い平衡化計算が必要である。
たとえ結晶構造からスタートする場合でも自由エネルギー計算を始める前に 10 ナノ秒以上の平衡化計算を実行して正しい熱平衡状態を得る事が必要である。図 3 は平衡化計
算中の FKBP と L14 の相互作用ポテンシャルで、20 ナノ秒以降で変化が少なくなって
熱平衡状態に到達した事が読み取れる。
結合自由エネルギーを求める為に、リガンドの他分子に対する相互作用が完全に存在する状態 λ=0 から相互作用が無くなる状態 λ=1 までの異なった中間状態を考える。それぞれの λ 点で分子動力学計算を実行して、0.1 ピコ秒毎に隣の λ に遷移する為に必要な仕事量を両方向について計算する。この計算では 11 個の λC 点でリガンドの Coulomb
相互作用を消去してから 21 個の λLJ 点で van der Waals 相互作用を消去した。
図 4 は FKBP と L2 リガンドの系で、van der Waals の
λLJ が 0.775 と 0.8 の間の仕事量を 500 ピコ秒間収集した時の分布である。forward は
λLJ を 大 き く す る 方 向 でreverse は逆である。大きい
分布は、初期運動量分布が異なる 12 個の分子動力学計
算全てから得られた仕事量分布で、小さい山はその中の一個の分子動力学計算か
図 4 仕事量分布ら得られた分布である。点線の矢印は仕事量の指数関数平均から得た自由エネルギーで、実線は両方向の仕事量分
布を用いてベネット受容比 (Bennett acceptance ratio: BAR) 法で計算した自由エネル
ギーである。BAR 法による自由エネルギーの値が小さい山と大きな山で 0.2 kcal/mol ほど異なる。もし各 λ 点で一つの分子動力学計算で自由エネルギーを計算した場合には32 個の λ 点があるので、容易に 6 kcal/mol も違う値を得る可能性がある。十分多くの仕事量から自由エネルギーを求める事が、非平衡仕事量の等式 (Jarzynski 等式) を使う上では極めて重要である。

-5-

SAR News No.17 (Oct. 2009)

4.

MP-CAFEE
タンパク質とリガンドの結合自由エネルギーを求める為に、32 個の λ点でそれぞれ 12 個の初期運動量分布が異なる分子動力学計算を二重消去法の式 (1) と (2) の二種類について行うので、水を含んだ合計で 768 個の分子動力学計算を並列で自由エネル
ギーの値が収束するまで実行する。収束する為に必要な時間は扱う系によって異なるが、FKBP だと 3 ナノ秒
までには全てのリガンドで収束した。
図 5 は FKBP と L2 の系で仕事量の指数関数平均から求めた自由エネル
ギーの積算である。点が使用した 32
個の λ 点を表している。(a) は L2 の水和計算で Coulomb と van der Waalsの λ=1 の値を加えると L2 の水和エネ
ルギー ∆Gsol になる。(b) は FKBP とL2 の複合系で、同様に λ=1 の値から
∆Gcomplex を得て、この二つの値から
結合自由エネルギー ∆Gbind が求まる。
(b) の複合系では forward と reverse で

図 5 指数関数平均による自由エネルギー

僅かな違いが認められるが、両方向の仕事量分布を用いて BAR 法で計算した値は二本の線の間に入る。この計算方法を MP-CAFEE (massively parallel computation of absolutebinding free energy with well-equilibrated states) と名付けた。従来の計算法では両方向で値が大きく異なる履歴現象 (hysteresis) が生じて仕舞う為に、タンパク質とリガンドの間にバネなどを導入してリガンドの動きを制限して hysteresis を人為的に小さくしたりしていた。MP-CAFEE ではこの様な作為的な方法を用いずに結合自由エネルギーを精確に求める事が出来る[4]。
図 6 は MP-CAFEE で計算した縦軸の FKBP リガンドの結合自由エネルギーを横軸の実験値[5] と比較している。計算は二種類の力場パラメータを用いて行っている。リガ
ンドの力場パラメータは general Amber 力場 (GAFF) を用いているが[6]、三角印はタン
パク質 FKBP に Amber99 を用いてリガンドの電荷には AM1-BCC を用いている。L13と L14 を除く 8 つのリガンドの三角印の値は 2005 年に発表したが[7]、その当時から何故 3.2 kcal/mol も実験値からシフトするのかが疑問であった。その後、L13 と L14 の計-6-

SAR News No.17 (Oct. 2009)

算を行い Amber99 の力場パラメータの精度に疑問を持つようなった。任意の有機分子に対して GAFF 力場パラメータを統一 的 に 割 当 て る FF-FOM (force field
formulator for organic molecules) を作成して、タンパク質とリガンドに共通した力場パラメータを割当てるようにした。リ
ガンドの電荷はタンパク質の Amber 電荷と 同 様 に RESP (restrained electrostatic
potential) で決めた。この統一的な力場で計算した結果が四角印である。Amber99

図 6

FKBP リガンドの結合自由エネルギーを用いた計算では大きく外れていた L13 と L14 も実験値に良く合うようなった。L6、
L8、L9 で作られる三角形も Amber99 で計算したよりも小さくなり実験値によく一致している。L2 と L3 の値は Amber99 を用いた値よりも改善はしているが、他のリガンドよりも実験値との一致が悪い。これは L2 と L3 が FKBP の疎水ポケットよりも小さい分子で、結合位置から離れても、疎水ポケットの中で回転したりしてポケットの中に留まる為に、酵素活性の阻害を行っている事によると考えられる。図 6 の実験値は阻害定
数 (inhibition constant) から導出されたエネルギーで、もし ITC (isothermal titration
calorimetry) などで結合エネルギーを直接計れば、L2 と L3 もより計算値に近い値になると予想している。

5.

プログラムと計算時間
MP-CAFEE の分子動力学計算にはフリーソフトとして配布されている GROMACS を使用した。分子動力学では原子間力の計算に多大な CPU 時間を要するが、GROMACS
は計算が重い部分がアセンブラー言語レベルで personal computer チップ向けに最適化されたコードを持っており、他の分子動力学プログラムと比較すると極めて計算速度が速い。この速い GROMACS を用いても MP-CAFEE で FKBP と一つのリガンドの結合自由エネルギーを求めるのに、3 GHz の Pentium 4 を一台使用したとして約 1200 日の計算時間が必要である。しかし、MP-CAFEE のアルゴリズムはそのままで並列化できるので、600 台の Pentium 4 を使えば 2 日ほどで計算が終わる。
Obama 大統領の science adviser に選任された David E. Shaw が率いる研究所では、分子動力学専用計算機の開発を進めている[8]。既に 512 並列の Anton 試作機が 2008 年 9月から稼働しており、タンパク質の動的挙動を調べる研究に使用されている。このように計算機の発展によって、生体高分子の挙動を精確に調べる為に必要な速度を持った計算機が最近漸く登場し始めた。またタンパク質やリガンドを記述する力場パラメータの高精度化も進んでおり[9]、大規模計算機シミュレーションが生物の基礎研究や創薬分-7-

SAR News No.17 (Oct. 2009)

野で活躍する環境が整いつつある。

謝辞
Stanford 大学の Vijay S. Pande 教授と富士通 BioServer プロジェクトのメンバーに謝意を表します。GROMACS の改造では Stockholm 大学の Erik Lindahl 教授にお世話になりました。この研究は NEDO の『バイオ・IT 融合機器開発プロジェクト』の支援を受け
ました。
参考文献
[1] Y. Duan and P. A. Kollman: Pathways to a Protein Folding Intermediate Observed in a
1-Microsecond Simulation in Aqueous Solution, Science, 282, 740-744 (1998)[2] C. D. Snow, H. Nguyen, V. S. Pande, and M. Gruebele: Absolute Comparison of Simulated
and Experimental Protein-Folding Dynamics, Nature, 420, 102-106 (2002)[3] C. Jarzynski: Nonequilibrium Equality for Free Energy Differences, Phys. Rev. Lett., 78,
2690-2693 (1997)[4] H. Fujitani, Y. Tanida, and Azuma Matsuura: Massively Parallel Computation of Absolute
Binding Free Energy with Well-Equilibrated States, Phys. Rev. E, 79, 021914 (2009)[5] D. A. Holt, J. I. Luengo, D. S. Yamashita, H. J. Oh, A. L. Konialian, H. K. Yen, L. W.
Rozamus, M. Brandt, M. J. Bossard, M. A. Levy, D. S. Eggleston, J. Liang, L. W. Schultz, T. J.
Stout, and J. Clardy: Design, Synthesis, and Kinetic Evaluation of High-Affinity FKBP Ligands
and the X-ray Crystal Structures of Their Complexes with FKBP 12, J. Am. Chem. Soc., 115,
9925-9938 (1993)[6] J. Wang, R. M. Wolf, J. W. Caldwell, P. A. Kollman, and D. A. Case: Development andTesting of a General Amber Force Field, J. Comput. Chem., 25, 1157-1174 (2004)[7] H. Fujitani, Y. Tanida, M. Ito, G. Jayachandran, C. D. Snow, M. R. Shirts, E. J. Sorin, and V.
S. Pande: Direct Calculation of the Binding Free Energies of FKBP Ligands, J. Chem. Phys.,123, 084108 (2005)[8] D. E. Shaw, M. M. Deneroff, R. O. Dror, J. S. Kuskin, R. H. Larson, J. K. Salmon, C. Young,
B. Batson, K. J. Bowers, J. C. Chao, M. P. Eastwood, J. Gagliardo, J. P. Grossman, C. R. Ho, D.
J. Ierardi, I. Kolossvary, J. L. Klepeis, T. Layman, C. McLeavey, M. A. Moraes, R.Mueller, E. C.
Priest, Y. Shan, J. Spengler, M. Theobald, B. Towles, and S. C. Wang: Anton, a Special-Purpose
Machine for Molecular Dynamics Simulation, Proc. Symp. Comp. Arch., 1-12 (2007)[9] H. Fujitani, A. Matsuura, S. Sakai, H. Sato, and Y. Tanida: High-Level ab initio Calculations
to Improve Protein Backbone Dihedral Parameters, J. Chem. Theory Comput., 5, 1155-1165(2009)-8-

SAR News No.17 (Oct. 2009)///// Cutting Edge /////

MM-PBSA 法を用いた精密結合自由エネルギー計算
北里大学薬学部

合田浩明、広野修一

はじめに

1.

近年、X 線結晶構造解析、NMR、電子顕微鏡などの構造解析技術および遺伝子工学技術の発展により、創薬の標的となるタンパク質の立体構造が比較的容易に決定されるようになった。それに伴い、標的タンパク質の立体構造に基づいた医薬分子設計 (Structure-Based Drug Design:
SBDD) 研究が理論的・計算化学的に大きく発展し、今や創薬戦略の中核をなしている。既に、上市された医薬品の開発において、SBDD が大きな役割を果たした例が数多く報告されている 1,2。
SBDD 研究における重要な要素の一つに、精密結合自由エネルギー計算がある。精密結合自由エネルギー計算は、主に SBDD 研究の最終段階（例えば、①de novo デザインなどにより考案された新規候補化合物の評価段階、および②スクリーニング計算などで発見したリード化合物の構
造最適化段階）において、①②

新規候補化合物はどの程度の結合親和性を有するのか？（絶対結合自由エネルギー、∆Gbind、の評価）
構造修飾により、リード化合物からどの程度結合親和性が改善されるか？（相対結合自由エネルギー、∆∆Gbind、の評価）を検討する際に、非常に重要となる。これまで、精密結合自由エネルギー計算を行うための手法
として、Free Energy Perturbation (FEP) 法と Thermodynamic Integration (TI) 法がよく使用されてきた 3,4。この二つの手法は、ある二つのリガンド間の∆∆Gbind 評価において、きわめて精度の良い結果を与えることができる。しかし、これらの手法には、・ ∆Gbind については評価できない。
・ 計算可能な系は、二つのリガンド間の化学構造的違いが非常に小さい場合(例えば、メチル基とエチル基の違いのような場合)だけである。
・ 計算中、あるリガンドからもう一つのリガンドに徐々に分子変換する必要があるため、非物理的な中間状態について長時間計算しなければならない（すなわち、計算時間がかかる）。
という短所がある。ところで、近年 Kollman 博士らにより、MM-PBSA (Molecular Mechanics
Poisson-Boltzmann Surface Area) 法を用いた精密結合自由エネルギー計算が提案された 5。FEP 法や TI 法と異なり、MM-PBSA 法を用いた計算には、・ ∆Gbind について計算できる。
・ そのため、化学構造が大きく異なるリガンド間の∆∆Gbind も評価できる。
・ FEP 法や TI 法よりも短時間で評価できる。
という長所がある。また FEP 法や TI 法ほどではないが、MM-PBSA 法を用いた計算も非常に精度が高い計算結果を与えることができる。それゆえ、現在では SBDD 研究において MM-PBSA法を用いた計算がさかんに行われるようになってきている 6,7。本稿では、先ず、MM-PBSA 法を用いた精密結合自由エネルギー計算の概略を示す。そして次に、適用例として、最近、我々の研究室で行われた天然由来のキチナーゼ阻害剤 Argifin および Argadin についての結合親和性解析を示す 8。

-9-SAR News No.17 (Oct. 2009)

2.

MM-PBSA 法を用いた精密結合自由エネルギー計算の概略

MM-PBSA 法を用いた精密結合自由エネルギー計算では、先ず、タンパク質－リガンド複合体、タンパク質単独、およびリガンド単独の溶液構造アンサンブルを用意する必要がある。したがって、それぞれについて水溶液中での分子動力学 (MD) シミュレーションを行い、それぞれの溶液構造アンサンブルを算出する必要がある。次に、複合体、タンパク質単独、およびリガンド単独の溶液構造アンサンブルを用いて、図１のような熱力学サイクルを考える。

このサイクルにおいて、求めるべき結合自由エネルギー ( ΔG bind (計算)) は、次のように表される。
ΔG bind (計算) = ΔG gas + G solv _ complex − G solv _ protein − G solv _ ligand

式 (1)

ここで、ΔG gas は気相中での結合エネルギーを表している。この項は、複合体、タンパク質単独、およびリガンド単独の溶液構造アンサンブルが持つ分子力学 (Molecular Mechanics) エネルギーを AMBER 等のパラメータ 9 を用いて計算し、差をとることで計算される。Over bar はアンサン
ブル平均をとることを示している。具体的に、 ΔG gas は次の項の和となる。
ΔG gas = ΔEint + ΔEVDW + ΔEelec − TΔSsolute

式 (2)

ここで、 ΔEint は結合時の構造変化に伴う内部エネルギー変化（結合長、結合角、二面角に関するエネルギー変化）の寄与、 ΔEVDW は van der Waals 相互作用エネルギーの寄与、 ΔEelec は静電相互作用エネルギーの寄与、−TΔSsolute は結合に伴う分子エントロピー変化の寄与、である。また、
式 (1) における、 G solv _ complex 、 G solv _ protein 、および G solv _ ligand は、それぞれ、複合体、タンパク質単独、およびリガンド単独の溶液構造アンサンブルについての水和自由エネルギーを表している。例えば、 G solv _ complex は、次のように極性項（ G PB _ complex ：水和自由エネルギーに対して電荷が寄与するエネルギー項）と非極性項（ G SA _ complex：水分子との van der Waals 相互作用エネル
ギーおよび空洞形成や水分子の再配置に必要なエネルギーを表す項）に分割されて、計算される。
G solv _ complex = G PB _ complex + G SA _ complex

式 (3)G PB _ complex は、Delphi10 等のプログラムを用いて Poisson-Boltzmann 方程式を数値的に解くことにより、G SA _ complex は表面積 (Surface Area) に依存した経験式により求められる。ところで、式 (1)中における、水和自由エネルギー項の寄与は、複合体の水和自由エネルギーから、タンパク質単独とリガンド単独の水和自由エネルギーを引いた形になっており、これはまさに結合に伴う水和
自由エネルギーの変化 ( ΔG solv ) を表している。
ΔG solv = G solv _ complex − G solv _ protein − G solv _ ligand-10-

SAR News No.17 (Oct. 2009)= G PB _ complex + G SA _ complex − (G PB _ protein + G SA _ protein ) − (G PB _ ligand + G SA _ ligand )= (G PB _ complex − G PB _ protein − G PB _ ligand ) + (G SA _ complex − G SA _ protein − G SA _ ligand )= ΔG PB + ΔG SA

式 (4)すなわち、MM-PBSA 法は熱力学サイクルを利用することにより、結合に伴う水和自由エネル
ギー変化をきちんと考慮している。それゆえ、精度の高い結合自由エネルギー値を与えることができる。最終的に得られる ΔG bind (計算) は次の項の和で計算される。
ΔG bind (計算) = ΔEint + ΔEVDW + ΔEelec − TΔSsolute + ΔG PB + ΔG SA

式 (5)以上のように、MM-PBSA 法を用いた結合自由エネルギー計算は、タンパク質－リガンド複合体、タンパク質単独、およびリガンド単独の溶液構造アンサンブルを用いて絶対結合自由エネル
ギー (∆Gbind) を直接的に評価する。したがって、化学構造が大きく異なるリガンド間の相対結合自由エネルギー (∆∆Gbind) 評価も可能となるため、SBDD において非常に有用な方法と言える。
次に、この MM-PBSA 法の特徴を利用した SBDD の実例を紹介する。

3.

キチナーゼ阻害剤 Argifin および Argadin の結合親和性解析

Argifin および Argadin は、共に、北里生命科学研究所において単離された天然由来の環状ペプ
チド性キチナーゼ阻害剤である (図２)11,12。キチナーゼ阻害剤には抗真菌薬、殺虫剤および喘息治療薬としての機能が期待できることが報告されており 13-18、現在、北里キャンパスにおいてこの二つの化合物に基づいた創薬研究が行われている。
ところで興味深いことに、
Argadin は Argifin
よりも非常に強いキチナーゼ阻害活性を示すことがわかっている (図２)。特に、霊菌のキチナー
ゼＢ (ChiB) に対する Argadin の阻害定数 (Ki = 20 nM) は、Argifin の阻害定数 (Ki = 33,000 nM)
よりも 1000 倍以上強い。
阻害定数は結合自由エネルギーに関連づけることができるので、Argadinの ChiB に対する結合自由エネルギー (∆Gbind(実験) = –10.92 kcal/mol) は、Argifin の値 (∆Gbind(実験) =–6.36 kcal/mol) より、4.56 kcal/mol も強いことになる。既に、Argifin–ChiB 複合体、およびArgadin–ChiB 複合体のＸ線結晶構造が報告されており (図３)、両者を比較することで Argifin とArgadin の結合様式の違いについては議論されていた 19。しかし結合様式の比較だけでは、結合親和性の違いを定量的に理解することはできない。そこで、我々は、精密結合自由エネルギー計
算を行い、エネルギー論の観点からこの違いを解析することにした。ところで、図２でわかるよ
うに Argifin と Argadin の化学構造は大きく異なっており、FEP 法や TI 法が適用できる系ではな
かった。したがって、MM-PBSA 法を用いた精密結合自由エネルギー計算を行うことにした。

先ず、Argifin–ChiB (Argadin–ChiB) 複合体の X 線結晶構造、および結晶構造から取り出した
Argifin (Argadin) 単独構造を用いて、それぞれ 1700ps の MD シミュレーションを行った。この計算にはプログラム AMBER 720 を用いた。構造が平衡に達したと思われる後半 1000ps から 10ps

-11-

SAR News No.17 (Oct. 2009)

毎に全部で 100 個のスナップショットを取り出し、それぞれの系の複合体およびリガンド単独の溶液構造アンサンブルとした。ところで、Argifin–ChiB (Argadin–ChiB) 複合体中の ChiB の構造と ChiB 単独の X 線結晶構造 21 を比較したところ、座標に関する RMSD 値が約 1Åと非常に低いものであった。このことは、Argifin および Argadin が結合しても ChiB はほとんど構造変化を起こさないことを示唆している。そこで本研究では、ChiB 単独の溶液構造アンサンブルは、複合体の溶液構造アンサンブルからリガンドを取り除くことで近似的に用意した。また、これまでに非常に多くの系において、この近似法が有効であることが報告されている 22-25。このタンパク質
単独構造に関する近似法を利用できる場合、タンパク質単独の MD シミュレーションを省くことができるので大幅に計算時間を短縮できる。
計算結果を表１に示す。得られた ΔG bind (計算) は、Argifin に対して–6.98 kcal/mol、Argadin に対して–11.16 kcal/mol であった。∆Gbind(実験) は、それぞれ、–6.36 および–10.92 kcal/mol であるので、
MM-PBSA 法を用いた計算が実験値を非常によく再現できるものであることがわかる。また、式(5) の各項を詳細に比較することで、結合過程における物理化学的性質を議論することができる。
例えば、ΔEint の項は Argifin ではほとんどゼロであるが、Argadin では結合に対して 5.26 kcal/mol
不利になっている。本計算では先程記したタンパク質単独構造に関する近似法を用いているので、ここでの ΔEint は、結合時のリガンドの構造変化に伴う内部エネルギー変化を意味している。すなわち、ChiB に結合する際に、Argifin はほとんど構造変化を起こさないが、Argadin は 5 kcal/mol
程度のエネルギー損失に相当する構造変化を引き起こすことを示している。図４に MD シミュ
レーションで得られた各リガンドの複合体中における構造 (結合配座) と単独溶液構造の比較を示す。確かに、Argifin ではその結合配座と単独溶液構造がほぼ同じであること、対照的にArgadin の結合配座はその単独溶液構造と大きく異なっていることがわかる。このように、化学構造が大きく異なる二つのリガンド間の比較では、リガンド単独の溶液構造アンサンブルを算出
し、結合時のリガンドの構造変化に伴う ΔEint を見積もることが大切である。また、結合自由エ
ネルギーに対して電荷が寄与する項 ( ΔG elec + PB ) は ΔEelec と ΔG PB の和で表されるが、この値は、Argifin および Argadin 共に正の値になっている。これは、両者の複合体形成は、静電的には不利であることを示している。したがって、両者の複合体形成は共に van der Waals 相互作用 ( ΔEVDW )と水和自由エネルギーの非極性寄与 ( ΔG SA ) により安定化されていることがわかる。

本手法により計算された Argifin と Argadin の間の ΔΔG bind (計算) は 4.18 kcal/mol であり、これも実験値、4.56 kcal/mol、をよく再現していた。ところで表１より、両者の複合体形成に重要な役割を果たしている ΔEVDW と ΔG SA が共に Argadin において Argifin より有利になっており、これらが Argadin の高親和性を生じさせていることがわかる。特に、Argadin の ΔEVDW は Argifin の値より約 12 kcal/mol も有利であり、これが主な要因となっている。図５に Argifin と ChiB の各アミノ酸残基との間の van der Waals 相互作用エネルギー値から Argadin についての値を差し引いたものをプロットした。負の値を持つ残基は Argifin と、正の値を持つ残基は Argadin とより強く vander Waals 相互作用している残基である。これより、E144、M212、W220、Y292、I339、および-12-

SAR News No.17 (Oct. 2009)

W403 の６個の残基が Argadin とより有利な van der Waals 相互作用を形成していることがわかる。
特に、W220 と W403 はそれぞれ、4.60 および 4.53 kcal/mol と非常に大きく有利になっていた。
この二つの値の和は 9.13 kcal/mol となり、Argifin と Argadin の間の van der Waals 相互作用エネ
ルギー差のほとんどを占めていることがわかる。したがって、ChiB に対する Argadin の高親和性は、主に Argadin とこの二つの Trp 残基との間の非常に有利な van der Waals 相互作用に起因すると考えられる。ところで、W220 を Ala 残基に置換した ChiB 変異体 (W220A 変異体) に対して、Argifin および Argadin はそれぞれ、–4.18 および–7.52 kcal/mol の結合自由エネルギー値を示すことが実験的に報告されている 19。これらの値とワイルドタイプの ChiB に対する∆Gbind(実験) との差をとることにより、W220 を Ala 残基に置換したことによる結合自由エネルギー損失は、
Argifin および Argadin に対してそれぞれ 1.51 および 3.40 kcal/mol と評価される。すなわち、Argadin についての損失のほうがより大きいことがわかる。これは、Argadin と W220 の相互作用
が Argadin の高親和性に大きく寄与していることを支持する実験結果である。

4.

おわりに

MM-PBSA 法に基づく精密結合自由エネルギー計算により、化学構造が大きく異なる Argifinと Argadin の結合親和性の違いをエネルギー論の観点から定量的に解析できた。このことは、FEP
法や TI 法とは異なり、本手法が非常に幅広い系に対して適用可能であることを示している。したがって、SBDD 研究において、MM-PBSA 法に基づく計算が今後ますます重みを増すであろう。

5.

謝辞

キチナーゼ阻害剤 Argifin および Argadin の結合親和性解析に関する研究は、科学研究費補助金、財団法人武田科学振興財団、財団法人持田記念医学薬学振興財団などの助成を受けて行ったものであり、その資金援助に深く感謝します。

6.

参考文献

1. Klebe, G., Drug Discov. Today, 11, 580-594 (2006).
2. Ahmed, M. et al., J. Pharmacol. Sci., 102, 167-172 (2006).
3. Beveridge, D.L. et al., Annu. Rev. Biophys. Biophys. Chem., 18, 431-492 (1989).
4. Kollman P., Chem. Rev., 93, 2395-2417 (1993).
5. Kollman, P.A. et al., Acc. Chem. Res., 33, 889-897 (2000).
6. Yasuo, K. et al., J. Chem. Inf. Model., 49, 853-864 (2009).
7. Brown, S.P. et al., J. Med. Chem., 52, 3159-3165 (2009).
8. Gouda, H. et al., Bioorg. Med. Chem., 16, 3565-3579 (2008).
9. Cornell, W.D. et al., J. Am. Chem. Soc., 117, 5179-5197 (1995).

-13-

SAR News No.17 (Oct. 2009)

10. Honig, B. et al., Science, 268, 1144-1149 (1995).
11. Ōmura, S. et al., J. Antibiot., 53, 603-608 (2000).
12. Arai, N. et al., Chem. Pharm. Bull., 48, 1442-1446 (2000).
13. Kuranda, M. J. et al., J. Biol. Chem., 266, 19758-19767 (1991).
14. Merzendorfer, H. et al., J. Exp. Biol., 206, 4393-4412 (2003).
15. Takaya, N. et al., Biosci. Biotechnol. Biochem., 62, 60-65 (1998).
16. Cohen, E., Arch. Insect. Biochem. Physiol., 22, 245-261 (1993).
17. Zhu, Z. et al., Science, 304, 1678-1682 (2004).
18. Sutherland, T.E. et al., Clin. Exp. Allergy., 39, 943-955 (2009).
19. Houston, D.R. et al., Proc. Natl. Acad. Sci. USA, 99, 9127-9132 (2002).
20. Case, D.A. et al., AMBER7, University of California, San Francisco, (2002).
21. van Aalten, D.M. et al., Proc. Natl. Acad. Sci. USA, 97, 5842-5847 (2000).
22. Massova, I. et al., J. Am. Chem. Soc., 121, 8133-8143 (1999).
23. Chong, L.T.; et al., Proc. Natl. Acad. Sci. USA, 96, 14330-14335 (1999).
24. Masukawa, K.M. et al., J. Med. Chem. 2003, 46, 5628-5637 (2003).
25. Gouda, H. et al., Biopolymers, 68, 16-34 (2003).

-14-

SAR News No.17 (Oct. 2009)///// Cutting Edge /////

QM/MM 分子動力学シミュレーションの生体高分子への応用
大阪大学蛋白質研究所 蛋白質情報科学研究系
米澤 康滋

1.

はじめに
生体高分子の分子シミュレーション研究は近年の計算機性能の著しい向上と相まって重要な研究手段となりつつある。現在、生体高分子シミュレーションに一般に用いられている分子力場（MM）は古典力学に基づいており、化学結合の切断や生成及び電子伝達等の量子効果を取り扱うことは出来ない。一方、現実の蛋白質の酵素反応や電子伝達機構などでは、この化学結合の切断や生成及び電子伝達等がごく普通に行われている。すなわちこのような現象を取り扱うためには
量子力学（QM）が必須であり、生体高分子シミュレーションで酵素反応及び電子伝達機構を解明するためには、QM の導入が不可欠である。分子の系に対しては QM 計算方法として分子軌道法が良く用いられているが、溶媒環境も含めた生体高分子の全系を分子軌道法で計算することは現在の計算機では不可能である。Warshel と Levitt は、分子力場だけでは実現が困難な酵素機能に関する分子シミュレーション研究を可能にする手法として、分子軌道法と分子力場のハイブリッド（QM/MM）シミュレーション法を初めて導入した[1]。彼らは、リゾチームの加水分解酵素反応
機構のシミュレーションで、リガンドの一部と活性部位の一部を半経験的分子軌道法で取り扱い、残りの部分を分子力場で取り扱って酵素反応を解析しその有用性を示した。このように QM/MM
シミュレーションは、研究対象分子の最も重要でかつ分子力場では記述できない領域を QM で取り扱い、残りの（相対的に）重要でない領域は MM で取り扱うことで、分子力場では得ることができない精密な分子シミュレーションを高速に実現する方法である。

2.

QM/MM ポテンシャルエネルギー
QM/MM シミュレーションのポテンシャルエネルギーは、以下の表式で表すことが出来る。（以後、分子軌道法で取り扱う領域を QM 領域、分子力場で取り扱う領域を MM 領域とする。）

Etot = EQM + EMM + EQM / MM
ここで Etot は QM/MM シミュレーション系の全エネルギー、EQM は QM 領域で分子軌道法を用いて計算したエネルギーで、EMM は MM 領域で分子力場を用いて計算したエネルギーである。最後の EQM/MM は QM 領域と MM 領域の相互作用エネルギーを表している。
QM/MM 境界に化学結合が存在しない場合には、QM/MM の相互作用項は大きく分けて２つの異なる取り扱い方法が存在する。ひとつは QM 領域と MM 領域が静電相互作用をしない方法で、もうひとつは QM 領域と MM 領域が静電相互作用する方法である。前者は初期の QM/MM シミュレーション
でしばしば用いられた方法で、プログラムへの実装が容易である特徴がある。後者は静電埋め込-15-

SAR News No.17 (Oct. 2009)

み法[2] と呼ばれており、現在の QM/MM シミュレーションではこの方法が広く用いられている。
静電埋め込み法では MM 領域にある原子の電荷は QM 領域の分子軌道計算に取り込まれて相互作用
エネルギーが計算される。この時 QM/MM 相互作用項は以下のように記述される。

EQM / MM = EvdW + Eele
ここで EvdW と Eele はそれぞれ QM 領域と MM 領域の原子間の van der Waals 相互作用項と静電相互作用項を示している。

3.

QM/MM 境界が化学結合に位置する場合の取り扱い
QM 領域と MM 領域の境界に化学結合が存在する場合は、QM/MM 相互作用の取り扱いは複雑になる。単純に QM/MM 境界で QM 領域を分割した場合 QM 領域末端には飽和されていない不対電子が生
じ QM 領域のエネルギーを大きく変化させるか、QM 領域の分子軌道計算自体を破綻させることもある。これを避けるために、QM 領域の末端を水素原子で終端するリンク原子法が最も簡単な処方として多用されている（図１）。

図１ リンク原子 L による QM/MM 境界の QM 領域終端例。Q と QB は QM 領域にある原子を示し、M と MBは MM 領域にある原子を示している。QB 原子と MB原子の間に QM/MM 境界が存在する場合に、リンク
原子 L は QM 領域の不対軌道を飽和させる目的で QB原子の近傍の適当な位置に配置される。

しかしこのリンク原子の導入は、(1) 系に本来存在しない余分な自由度を加える。(2) リンク
原子と MM 原子の相互作用によって QM 領域に過剰分極が誘起される、等の困難を生じることが報告されている。そのためリンク原子を用いる QM/MM シミュレーションでは、リンク原子を QB 原子と MB 原子の結合上に一定距離で束縛したり、MB 原子の電荷を消去するなどの処方が必要となる[7]。リンク原子に起因するこのような困難を避けるために、Gao 等による Generalized Hybrid
Orbital（GHO）法[3] や Rivail 等や Friesner 等による Localized Self Consistent Field（LSCF）
法[4]、Frozen Orbital 法[5] がそれぞれ提案されている。いずれの方法もリンク原子を使用せず QM/MM 境界に人工的な分子軌道を置くことでその困難を回避しているが、QM/MM 境界の問題を完全には解決しておらず今後の進展が期待される分野である。尚、これらの QM/MM 境界に対する方法は、リンク原子法も含めて、問題に応じて注意深くパラメータを設定することでほぼ等しい
精度を与えることが報告されている[6]。

-16-SAR News No.17 (Oct. 2009)

生体高分子への QM/MM シミュレーション応用例 – プロリンのシス－トランス異性化反応機構

4.

QM/MM シミュレーションの生体高分子への応用例として、我々が最近報告したプロリン残基のシス－トランス異性化反応機構の研究を紹介する[7]。X 線結晶構造解析や NMR 溶液構造解析か
ら、殆ど全てのアミノ酸が蛋白質の立体構造中でトランス状態を取るのに対してプロリンはおよ
そ 10 ％がシス状態を取っていることが知られている。このためプロリンのシス－トランス異性化は以前から蛋白質フォールディングの律速過程として重要視され研究が進められている。 最近、これに加えてプロリンのシス－トランス異性化反応は蛋白質機能の切り替えスイッチとして細胞内の信号伝達の調整機能を直接または間接的に司ることが明らかになり、新たに注目を集めている[8]。一般にアミノ酸残基の分子力場は、蛋白質中での安定状態を記述することに主眼が置かれており異性化反応などの遷移状態を記述することは考慮されておらずその精度は保証されていない。そこで我々は、このプロリンのシス－トランス異性化反応に QM/MM ポテンシャルを適用して分子動力学シミュレーション（QM/MM MD）することで、反応の自由エネルギーを精密に求めてこの異性化機構の詳細を研究した。
我々はプロリンのモデルとして N-acetyl-N'-methylprolineamide (Ace-Pro-NMe) プロリンジ
ペプチド分子を用いた。このペプチド全体を QM 領域とし、水の分子力場モデル TIP3P の MM 領域
でこのペプチドの周りを取り囲み、シミュレーション系を構築した（図 2）。ここで QM 領域の計
算は、制限ハートリーフォック法を採用し、QM/MM 間相互作用計算は静電埋め込み法を用いている。我々は自由エネルギーを求める反応座標として、
図 2（左）にしめす捻り角 ω’（∠Cδ-N-C-O）を使用した。この反応座標上では捻り角 0°がシス状態に、180°がトランス状態に対応している。

図 2（ 左 ） N-acetyl-N'-methylprolineamide(Ace-Pro-NMe) プロリンジペプチド、黄色の波線
で捻り角 ω’を示した。
（右）プロリンジペプチドの QM/MM シミュレーション系：TIP3P の数は 1355
分子で中心から球対象調和ポテンシャルで束縛されている。系の半径は 19Åである。

安定状態のペプチド結合の電子状態は sp2 混成を成して共鳴しているためにその C-N 結合長は
通常のσボンド結合長よりも 0.1Å程短縮している。これに加えてシス－トランス間にはこの共鳴による大きなエネルギー障壁が存在してペプチド結合を平面に保っている。すなわち室温環境
でこのエネルギー障壁を超える現象は非常にまれにしか実現せずこの異性化過程を現在の計算機でカノニカルシミュレーションすることは困難である。そこで我々は、反応座標に沿って人工的なポテンシャルを適用し、シミュレーションがまれにしか訪れない位相空間領域に局所化する方法（アンブレラポテンシャル法）を採用しこの困難を回避した。尚、この人工的なポテンシャ-17-

SAR News No.17 (Oct. 2009)

ルの影響は解析時に統計熱力学的に矛盾無く取り除くことができる。 我々は実際のアンブレラ
ポテンシャルとして、反応座標の捻り角 ω’をある一定値に束縛するような調和型ポテンシャ
ルを用いている。
我々は反応座標を 36 分割して、異なる束縛値 ω’を持つ 36 の独立したアンブレラサンプリ
ングシミュレーションを実行しこれらのシミュレーションから得られたデータを解析した。 以
下に示す結果は[7] の基底関数 4-31G を 6-31G に改良し、新たにシミュレーションして得たデー
タによる結果である。図 3（左）にシミュレーションから得られたプロリンジペプチドのシス－
トランス異性化自由エネルギー（Potential of Mean Force）の反応座標依存性を示す。

図 3（左）プロリンジペプチドのシス－トランス異性化自由エネルギー地形、反応座標 0°がシス状態に、180°がトランス状態に対応する。アンブレラサンプリングで分割された自由エネルギーはアンブレラ積分法[9] で結合して連続した自由エネルギー地形を得ている。（右）プロリンの C-N 結合長の反応座標依存性。

得られた自由エネルギーはトランス状態がシス状態よりも 3-4 kcal/mol 程安定であることを示し実験と一致している。さらにエネルギー障壁の値（23 kcal/mol）も実験値と良く一致して
おり QM/MM シミュレーションの精度の良さを示している。 次にプロリンの平均 C-N 結合長の反応座標依存性を図 3（右）に示す。図から、シス状態とトランス状態では共鳴によって 1.34Å程の結合長であったプロリンの C-N 結合長がシス－トランス異性化の進行に依存して変化し、シス状態とトランス状態の中間状態では最大およそ 1.42Åに伸張していることが解る。この結果は
シス及びトランス状態では sp2 混成していた電子状態がシス－トランス異性化の中間状態では
sp3 混成電子状態に変換していることを明確に示している[10]。
ここで、プロリンの窒素原子を含むペプチド結合の平面性を測る尺度として pyramidality（以下、ピラミダリティ）= det(e1,e2,e3) を導入する。e1、e2 と e3 は、プロリン環の窒素原子 Nからそれぞれ Cδ、Cα、C 原子に向かう単位ベクトルとする。det はこれらのベクトルの行列式を意味する。 このピラミダリティは N 原子を頂点とするピラミッド様四面体の体積に比例しペプチド結合の平面性からのずれを示す良い指標となる。図 4 に QM/MM MD シミュレーションから得られたプロリンジペプチドの統計平均したピラミダィティ値の反応座標依存性を示す。
ピラミダィティは反応座標 0°付近（シス状態）では、ほぼゼロの値（平面構造）を持つ。そ-18-

SAR News No.17 (Oct. 2009)

して反応座標が正に増加するにつれて減少し 60°付近で一旦負の極値（負のピラミッド構造：
図 4(b)）を取る。そして反応座標がさらに増加すると急激に正の極値（正のピラミッド構造：
図 4(a)）に跳躍し、その後再び減少に転じて 180°（トランス状態）で再びゼロ（平面構造）に戻っている。反応座標が 180°から 360°に増加する際にも同様の現象をみることができる。

図 4 平均ピラミダリティの反応座標依存性と、(a) 正のピラミッド様コンフォメーション及び、(b) 負のピラミッド様コンフォメーションのスナップショット。

以前から分子軌道法によって正のピラミッド様コンフォメーションがシス－トランス異性化の中間状態として報告されているが[10]、我々の研究が明らかにしたプロリン異性化に伴う正と負の間のピラミッド様コンフォメーションの急激な跳躍は QM/MM MD シミュレーションで初めて予見された新規な現象である。さらに我々は詳細な解析の結果、プロリン窒素原子に関する溶媒
水分子の水素結合がこのピラミッド様コンフォメーションの安定化に寄与していることを解明した。正のピラミッド様コンフォメーションは、図４(a) に示されている（水色）矢印方向からの溶媒水分子による水素結合で安定化されており、負のピラミッド様コンフォメーションは図4(b) に示されている矢印方向からの水素結合で安定化されている。これらの水素結合は競合しておりその度合いはプロリンジペプチドの異性化に依存している。異性化の過程でこの競合関係は複雑に変化して、競合のバランスが反応座標の 60°付近と 250°付近で急激に逆転するために
正負のピラミッド様コンフォメーションの跳躍を引き起こす一因となると我々は考えている。

5.

まとめ
QM/MM シミュレーションの理論的基礎と応用例について概説した。QM/MM シミュレーションは分子力場の精度と限界に支配される分子シミュレーション（MM）と、その計算コスト故に計算対
象サイズが限定される分子軌道法（QM）の両者の欠点を補い合いかつ双方の利点を有効活用でき-19-

SAR News No.17 (Oct. 2009)

る優れた方法である。さらに QM/MM シミュレーションで注意を要する QM/MM 境界の取り扱いに言及しリング原子法に代表される化学結合を跨ぐ QM/MM 境界の様々な方法について紹介した。最後に溶媒中のプロリンジペプチドのシス－トランス異性化反応に関する QM/MM MD シミュレーショ
ンの結果を具体的研究例として紹介した。ここでは、分子力場を使用する分子動力学法と分子軌
道法の各々を独立に使用することでは解明し得なかった現象が QM/MM MD シミュレーションを適用することで初めて明らかになった例を挙げて QM/MM シミュレーションの有用性を示した。今後
QM/MM シミュレーションは、生体高分子やナノの研究領域でその重要性を増していくことが予想
される。

6.

謝辞
本研究に関しましてご協力及びご指導頂きました、阪大蛋白研 中村春木教授、理研 高田俊和

先生、そして共同研究者の皆様に深く感謝申し上げます。

7.

参考文献

[1] Theoretical Studies of Enzymatic Reactions: Dielectric Electrostatic and Steric Stabilization of the Carbonium Ion
in the Reaction of Lysozyme, A. Warshel and M. Levitt, J. Mol. Biol. 103, 227-249 (1976).
[2] QM/MM: What Have We Learned, Where are We, and Where Do We Go from Here?, H. Lin and D. G. Truhlar,
Theoret. Chem. Acc. 117, 185-199 (2007).
[3] A Generalized Hybrid Orbital (GHO) Method for the Treatment of Boundary Atoms in Combined QM/MM
Calculations, J. Gao, P. Amara, C. Alhambra and M. J. Field, J. Phys. Chem. A 102, 4714-4721 (1998).
[4] Quantum Chemical Computations on Parts of Large Molecules: The Ab Initio Local Self Consistent Field Method,
X. Assfeld and J.-L. Rivail, Chem. Phys. Lett. 263, 100-106 (1996).
[5] Ab Initio Quantum Chemical and Mixed Quantum Mechanics/Molecular Mechanics (QM/MM) Methods forStudying Enzymatic Catalysis, R.A. Friesner and V. Guallar, Annu. Rev. Phys. Chem. 56, 389-427 (2005).
[6] Frontier Bonds in QM/MM Methods: A Comparison of Different Approaches, N. Reuter, A. Dejaegere, B.
Maigret and M. Karplus, J. Phys. Chem. A 104, 1720-1735 (2000).
[7] Intra- and Intermolecular Interaction Inducing Pyramidalization on Both Sides of a Proline Dipeptide during
Isomerization: An Ab Initio QM/MM Molecular Dynamics Simulation Study in Explicit Water, Y. Yonezawa, K.
Nakata, K. Sakakura, T. Takada and H. Nakamura, J. Am. Chem. Soc. 131, 4535-4540 (2009).
[8] Prolyl Isomerase Pin1: New Findings on Post-translational Modifications and Physiological Substrates in Cancer,
Alzheimer's Disease and Asthma, K. Takahashi, C. Uchida, R. W. Shin, and T. Uchida, Cell Mol. Life Sci. 65,
359-375 (2008).
[9] Bridging the Gap between Thermodynamic Integration and Umbrella Sampling Provides a Novel Analysis
Method: Umbrella Integration, J. Kästner and W. Thiel, J. Chem. Phys. 123, 144104 (2005).
[10] A Computational Study of the Isomerization of Prolyl Amides As Catalyzed by Intramolecular Hydrogen
Bonding, K. N. Rankin and R. J. Boyd, J. Phys. Chem. A 106, 11168-11172 (2002).

-20-

SAR News No.17 (Oct. 2009)///// Activities /////
構造活性フォーラム 2009 開催報告
「化合物と標的・非標的タンパク質との相互作用－創薬標的の同定と分子設計－」
構造活性フォーラム 2009
実行委員長 竹田－志鷹真由子

平成 21 年 6 月 19 日（金）、北里大学薬学部コンベンションホール（港区白金 5－9－1）
において、構造活性フォーラム 2009 が標記の討論主題で開催されました。構造活性フォー
ラムは、1999 年に構造活性相関講習会として開始されて以来（2003 年からは構造活性フォー
ラムと名称を変更）、今回で 11 回目を迎えました。参加者は、総数 116 名（講師・実行委員・
学生アルバイトを含む）でした。当日のプログラムは以下の通りです。
1. 合田浩明（北里大学薬学部）
ヒト酸性キチナーゼを標的にした in-silico 創薬研究2. 石黒正路（新潟薬科大学応用生命科学部）
GPCR におけるリガンド受容の多様性
3. 奥野恭史（京都大学大学院薬学研究科）
多重標的創薬のためのインフォマティクス
4. 大田雅照（中外製薬㈱）Unwanted Interactions between Small Molecule and Multiple Proteins: Selectivity and Off-targets
5. 夏目徹（産業技術総合研究所）
タンパク質ネットワーク解析から展開するケミカルバイオロジー
6. 藤田稔夫（京都大学名誉教授）
多様な標的に対応する個々の構造活性関係の間の上位の関係の考察と解析 －SAR-omicsの提唱－
生体システムは生命現象を司る多種多様な分子の複雑なネットワークから成っており、創薬研究などにおいてもそれらの相互作用情報の活用が必須となっています。今回のフォーラムでは、化合物－タンパク質間およびタンパク質－タンパク質間の相互作用の多様性に焦点をあて、分子標的インシリコ創薬の実際、創薬標的としての GPCR、多重標的創薬のためのインフォマティクス、マルチターゲット化合物の選択性向上と副作用毒性回避、
タンパク質ネットワーク解析、構造活性相関研究との関連などについて、最先端でご活躍の 6 名の先生方にご講演いただきました。いずれのセッションも会場からの質疑も活発で大変有益な討論となり、盛会のうちに終了することができました。
末筆ながら、ご講演いただいた講師の先生方、ご参加いただいた皆様、開催にあたりご支援ご協力をいただきました皆様方に心より御礼申し上げます。

-21-

SAR News No.17 (Oct. 2009)///// Activities /////SAR Promotion Award 2009 年度受賞者（庶務幹事 新潟薬科大学 米田照代）

構造活性相関部会では，2005 年度より構造活性相関研究の発展を促進するための事業として当該制度を設け，部会員の国外での研究発表を奨励している．2009 年度は，6 月 19 日の常任幹事会において，受賞者を次の 1 名に決定した．

氏

名

吉田 達貞（よしだ たつさだ）

所

属

徳島大学大学院 ヘルスバイオサイエンス研究部 生体情報薬科学部門 分子情報薬学講座 創薬理論化学分野（助教）

参加学会名 Fifth International Symposium on Computational Methods in Toxicology andPharmacology Integrating Internet Resources (CMTPI2009)
開催期日

2009 年 7 月 4 日～8 日

開催場所

トルコ，イスタンブール

演題

Comparative QSAR Analysis of a Series of Benzene Sulfonamide Inhibitors Using Ab
Initio Fragment MO Calculation of Their Complex Structures with Carbonic Anhydrase（ポスター発表）

受賞者の報告は，次号の SAR News に掲載される予定である．

-22-

SAR News No.17 (Oct. 2009)///// Activities /////＜会告＞第 37 回構造活性相関シンポジウム
日時

平成 21 年 11 月 12 日（木）・13 日（金）

会場
主催
共催
懇親会
連絡先

北里大学薬学部コンベンションホール（東京都港区白金 5-9-1）
日本薬学会構造活性相関部会
日本化学会 日本分析化学会 日本農芸化学会 日本農薬学会
11 月 12 日（木）第 37 回構造活性相関シンポジウム事務局〒105-0014 東京都港区芝 3-17-15 クリエート三田 207
Tel 03-3798-5253
Fax 03-3798-5251 E-mail sar2009@event-convention.com

第 1 日目（11 月 12 日）
10:00 - 10:10 開会
10:10 - 11:00

横山祐作（東邦大・薬）一般講演（会場：薬学部コンベンションホール）
座長 ： 大田雅照

KO-01

KO-02

新型インフルエンザウィルス(H1N1)のモデリングデータベースの開発○寺師玄記、加納和彦、中村裕樹(北里大・薬)、高谷大輔、松本武久(理研)、竹田－志鷹真
由子、梅山秀明(北里大・薬、理研)
ドッキング条件検討システム(PALLAS)の開発と検証○幸瞳(理研・SSBC)、石田剛(理研・SSBC、日鉄日立)、佐藤朋広（理研・SSBC 、東大院・理）、津金沢恵子、田仲昭子、本間光貴（理研・SSBC）

11:00 - 11:40

招待講演（会場：薬学部コンベンションホール）
座長 ： 石黒正路

KI-01

医薬品に秘められたキラリティー
夏苅英昭（帝京大・薬）13:10 - 15:40 ポスターセッション（会場：薬学部 1 号館 1201 学生自習室）
13:10 - 14:25 奇数番発表
14:25 - 15:40 偶数番発表
15:40 - 16:45

一般講演（会場：薬学部コンベンションホール）
座長 ： 加藤博明

KO-03

KO-04
KO-05

サポートベクターマシンによる多種類の化学物質の発ガン性の予測○田辺和俊、栗田多喜夫(産総研・脳神経情報)、小野寺夏生(筑波大院・図情メディア)、貝原巳樹雄(一関工業高専)、鈴木孝弘(東洋大・経)有害性評価支援システム統合プラットフォームの開発○櫻谷祐企、佐藤佐和子、西川智、山田隼、前川昭彦、林真(NITE)
単純グラフ同型な NTG ネットワークを用いた薬物構造データマイニング○野村京平、大田黒空、高橋由雅(豊橋技科大)-23-

SAR News No.17 (Oct. 2009)

17:00 - 17:40

招待講演（会場：薬学部コンベンションホール）
座長 ： 粕谷敦

KI-02

新規物質の分子標的予測システムの開発と抗がん剤創薬への応用
矢守隆夫（癌研究会・癌化学療法センター・分子薬理）17:40 - 17:55

ポスター賞発表（会場：薬学部コンベンションホール）

18:30 - 20:30

懇親会（会場：恵比寿ガーデンプレイス ビヤステーション）
ポスター賞授賞式

ポスターセッション（会場：1201 学生自習室）
KP-01 HIV-1 膜融合阻害ペプチドにおける gp41/S138A 置換効果の MM と PBSA 解析○寺川幸宏、渡部毅、渡辺健太郎、大野浩章(京大・薬)、泉和樹(京大・ウイルス研)、大石真也、藤井信孝、北浦和夫(京大・薬)KP-02 取り下げ
KP-03 ドッキング条件検討システム(PALLAS)の開発と検証○幸瞳(理研・SSBC)、石田剛(理研・SSBC、日鉄日立)、佐藤朋広（理研・SSBC 、東大院・理）、津金沢恵子、田仲昭子、本間光貴（理研・SSBC）
KP-04 ファーマコフォアの空間配置に基づく低分子-タンパク質間相互作用記述子（Pharm-IF）の開発と機械学習を用いたインシリコスクリーニングへの応用○佐藤朋広(東大院・理、理研)、幸瞳、本間光貴(理研)、横山茂之(東大院・理、理研)
KP-05 タンパク質－リガンドドッキングにおけるリガンド初期配座の影響○小田彰史(東北薬大・薬、阪大・蛋白研)、小林佳奈、高橋央宜(東北薬大・薬)、山乙教之、
広野修一(北里大・薬)KP-06 L 糖含有糖鎖の合成とその線虫ガレクチンとの結合活性の検討
西山和沙、○山田篤、高橋美貴、武内智春、笠井献一、夏苅英昭、高橋秀依(帝京大・薬)KP-07 炭酸脱水酵素－ベンゼンスルホンアミド阻害剤複合体の非経験的フラグメント分子軌道法を用いた相互作用解析と QSAR
○宗井陽平、島本和典、岡田耕平、吉田達貞、中馬寛(徳島大院・薬)KP-08 フェノール水素原子のラジカル引き抜き反応の非経験的分子軌道法解析に基づくフラボノイ
ドの DPPH 活性の QSAR
○廣隅公治、吉田達貞(徳島大院・薬)、河合慶親、寺尾純二(徳島大院・栄養)、中馬寛(徳島大院・薬)KP-09 N1 ノイラミニダーゼ-シアル酸誘導体複合体の非経験的フラグメント分子軌道法計算に基づく相互作用解析
○比多岡清司、郡恵理、原田政隆、吉田達貞、中馬寛(徳島大院・薬)KP-10 分子動力学計算に基づく Arg344 のヒト Cathepsin A 活性への影響の解析○郡恵理、比多岡清司、原田政隆(徳島大院・薬)、門田佳人(徳島文理大・薬)、堀川靖、吉田達貞、伊藤孝司、中馬寛(徳島大院・薬)KP-11 In silico screening による非ペプチド性β-セクレターゼ阻害剤の探索○松下泰雄、中込泉、山乙教之、広野修一(北里大・薬)KP-12 カゼインキナーゼ 2 複合体に対する構造最適化と結合エネルギー評価○浅田直也、北浦和夫(京大院・薬)-24-

SAR News No.17 (Oct. 2009)

KP-13
KP-14
KP-15
KP-16

KP-17
KP-18

KP-19
KP-20
KP-21
KP-22
KP-23

KP-24

KP-25

KP-26
KP-27
KP-28
KP-29
KP-30

ファルマコフォアによる化合物多様性評価○石川誠(日産化学)
取り下げ
AutoGPA モデルによる受容体ポケットの特徴抽出○朝川直行、小林誠一、後藤純一(菱化システム)、平山令明(東海大・医)アルドース還元酵素阻害活性を有する新規三環系化合物の Structure-based 3D-QSAR○生川佳代、山乙教之、合田浩明(北里大・薬)、竹田大輔、峰平大輔(富山大病院・薬)、王旭(富山大院・薬)、加藤敦、足立伊左雄(富山大病院・薬)、松谷裕二、豊岡尚樹(富山大院・薬)、広野修一(北里大・薬)非経験的分子軌道法による分子間相互作用解析に基づく Hammett σの理論的考察
清水美帆、○吉田達貞、中馬寛(徳島大院・薬)
シクロフィリン A を標的とする抗 HIV 薬の開発を指向した構造活性相関研究○杉本裕昌、田雨時、岡本晃典(阪大院・薬)、川下理日人(阪大院・薬、阪大微研)、安永照雄(阪大微研)、高木達也(阪大院・薬、阪大微研)
パラオキシ安息香酸エステル類の構造とモル毒性の関係に関する研究
石原良美(東海大・理)、○片岡啓一(東海大院・理)、齋藤寛、髙野二郎(東海大・理)
化学物質の反復投与毒性試験データの解析とカテゴリーアプローチへの応用○西川智、櫻谷祐企、佐藤佐和子、山田隼、前川昭彦、林真(NITE)
カテゴリーアプローチによる化学物質の生物濃縮性予測○池永裕、櫻谷祐企、佐藤佐和子、山田隼(NITE)
ヒト発癌性物質の解析と Support Vector Machine による予測モデルの構築○三輪祐太朗、坂本久美子、山内あい子(徳島大院・薬)Identification of the 3D-Pharmacophore for Human ABCG2 Transporter Substrates○山乙教之(北里大・薬)、齊藤光(東工大院・生命理工)、石川智久(理研)、広野修一(北里大・薬)置換安息香酸アシルグルクロニドの分解速度定数に関する電子的及び立体的記述子を用いた構造活性相関○吉岡忠夫、馬場暁子(北海道薬大)N-アセチルトランスフェラーゼ 2 の構造と薬物動態○大倉一人、生城山勝巳(千葉科大・薬)、篠原康雄(徳大・疾患ゲノム研)、堀均(徳大・ソシオ
テクノサイエンス研)
新規経口急性毒性 Alert の探索○赤堀有美(化学物質評価研究機構(CERI))
ヒト酸性キチナーゼに対して特異性を有する Argifin 誘導体阻害剤の論理的分子設計○合田浩明、寺嶋真一、広野修一(北里大・薬)トリプレット表現に基づくヒト匂い受容体の構造特徴解析○森下千識、加藤博明(豊橋技科大)ドッキング計算と熱力学測定による血清アルブミンの構造センシング○加藤祐樹、田浦俊明(愛知県立大院・情報)ドッキングスタディを利用した抗インフルエンザ薬の薬剤耐性評価法○安田匡志(阪大院・薬)、川下理日人(阪大院・薬、阪大微研)、柏田理恵、田雨時、岡本晃典(阪大院・薬)、川瀬雅也(長浜バイオ大)、安永照雄(阪大微研)、高木達也(阪大院・薬、阪
大微研)

-25-

SAR News No.17 (Oct. 2009)

KP-31

KP-32

KP-33
KP-34
KP-35

KP-36

ホモロジーモデリングとドッキング手法による thiolactomycin と D-aspartate oxidase の結合様式の推定
○中込泉、山乙教之、合田浩明、片根真澄、本間浩、広野修一(北里大・薬)GPCR structure modelling and ligand docking: GPCRDock2008○加納和彦(北里大・薬)、高谷大輔(理研)、寺師玄記、竹田-志鷹真由子、梅山秀明(北里大・薬、理研)知識ベースによる合成経路を包括した化合物バーチャルライブラリの構築○西村拓朗、船津公人(東大院・工)単純グラフ同型な NTG ネットワークを用いた薬物構造データマイニング○野村京平、大田黒空、高橋由雅(豊橋技科大)バーチャルスクリーニングを用いた新規脱皮ホルモン受容体リガンドの探索○原田俊幸、中川好秋(京大院・農)、山田豊(サイバネット)、大江武弘(オープンアイ)、宮川恒(京大院・農)事例データにもとづく薬物分子グラフの NTG ネットワークの生成と応用○寺本岳史、高橋由雅(豊橋技科大)第 2 日目（11 月 13 日）
9:00 - 10:20
一般講演（会場：薬学部コンベンションホール）
座長 ： 藤原巌
KO-06
KO-07

プロテインキナーゼ ATR と阻害剤の複合体モデリング○米田照代、古寺弘昭、田宮実、石黒正路(新潟薬大・応用生命)
カゼインキナーゼ 2 複合体に対する構造最適化と結合エネルギー評価○浅田直也、北浦和夫(京大院・薬)座長 ： 清水良

KO-08

KO-09

N1 ノイラミニダーゼ-シアル酸誘導体複合体の非経験的フラグメント分子軌道法計算に基づく相互作用解析
○比多岡清司、郡恵理、原田政隆、吉田達貞、中馬寛(徳島大院・薬)非経験的フラグメント分子軌道法による薬物－受容体相互作用解析と QSAR への応用
○吉田達貞、宗井陽平、中馬寛(徳島大院・薬)10:40 - 11:40

特別講演（会場：薬学部コンベンションホール）
座長 ： 横山祐作

KS-01

ユニークな光化学原理に基づくバイオイメージングプローブの開発研究
大学発の創薬研究の新たな動き
長野哲雄（東大院・薬、東大・生物機能制御化合物ライブラリー機構）

13:10 - 14:15
KO-10
KO-11

一般講演（会場：薬学部コンベンションホール）

座長 ： 黒木保久
AutoGPA:グリッドとファーマコフォアモデルの連携による 3D-QSAR モデルの自動構築○小林誠一、朝川直行、後藤純一(菱化システム)、平山令明(東海大・医)シクロフィリン A を標的とする抗 HIV 薬の開発を指向した構造活性相関研究○杉本裕昌、田雨時、岡本晃典(阪大院・薬)、川下理日人(阪大院・薬、阪大微研)、安永照雄(阪大微研)、高木達也(阪大院・薬、阪大微研)-26-

SAR News No.17 (Oct. 2009)

KO-12

多変量解析による化学物質の加水分解反応予測○栗花落昇平(阪大院・薬)、尾形直紀、藤本貴男(阪大・薬)、山崎広之、日高伸之介、岡本晃典(阪大院・薬)、川下理日人(阪大院・薬、阪大微研、阪大微研・感染症共同研)、、高木達也(阪大院・薬、阪大微研、阪大微研・感染症共同研)

14:30 - 16:20

一般講演（会場：薬学部コンベンションホール）
座長 ： 中山章

KO-13

KO-14

バーチャルスクリーニングを用いた新規脱皮ホルモン受容体リガンドの探索○原田俊幸、中川好秋(京大院・農)、山田豊(サイバネット)、大江武弘(オープンアイ)、宮川恒(京大院・農)
薬物排出に関わる P-糖タンパク質の基質認識機構 －さまざまな化学物質の構造と ATPase
活性との関係－○赤松美紀、金岡怜志、木村泰久、中川好秋、植田和光（京大院・農）
座長 ： 岡島伸之

KO-15
KO-16

KO-17

KO-18

オーファン GPCR のリガンド探索手法の開発○光山倫央、荒川正幹、船津公人(東大院・工)Organic anion transporting polypeptide (OATP) 1B1 と OATP1B3 の基質を用いたファーマコフ
ォア解析、3D-QSAR 解析○渡辺悦郎(東大院・薬)、山乙教之(北里大・薬)、楠原洋之(東大院・薬)、広野修一(北里大・薬)、杉山雄一(東大院・薬)座長 ： 久保寺英夫
Identification of the 3D-Pharmacophore for Human ABCG2 Transporter Substrates○山乙教之(北里大・薬)、齊藤光(東工大院・生命理工)、石川智久(理研)、広野修一(北里大・薬)GPCR structure modelling and ligand docking: GPCRDock2008○加納和彦(北里大・薬)、高谷大輔(理研)、寺師玄記、竹田-志鷹真由子、梅山秀明(北里大・薬、理研)16:20 - 16:30

閉会（会場：薬学部コンベンションホール）-27-

SAR News No.17 (Oct. 2009)///// Activities /////
構造活性相関部会の沿革と趣旨
1970 年代の前半、医農薬を含む生理活性物質の活性発現の分子機構、立体構造・電子構造の計算や活性データ処理に対するコンピュータの活用など、関連分野のめざましい発展にともなって、構造活性相関と分子設計に対する新しい方法論が世界的に台頭してきた。このような情勢に呼応するとともに、研究者の交流と情報交換、研究発表と方法論の普及の場を提供することを目的に設立されたのが本部会の前身の構造活性相関懇話会である。1975 年 5 月京都において第１回の ｢懇話会｣（シンポジウム）が旗揚げされ、1980 年からは年１回の ｢構造活性相関シンポジウム｣ が関係諸学会の共催の下で定期的に開催されるようになった。
1993 年より同シンポジウムは日本薬学会医薬化学部会の主催の下、関係学会の共催を得て行なわれることとなった。構造活性相関懇話会は 1994 年にその名称を同研究会に改め、シンポジウム開催の実務担当グ
ループとしての役割を果すこととなった。2002 年 4 月からは、日本薬学会の傘下組織の構造活性相関部会
として再出発し、関連諸学会と密接な連携を保ちつつ、生理活性物質の構造活性相関に関する学術・研究の振興と推進に向けて活動している。現在それぞれ年一回のシンポジウムとフォーラムを開催するとともに、部会誌の SAR News を年二回発行し、関係領域の最新の情勢に関する啓蒙と広報活動を行っている。
本部会の沿革と趣旨および最新の動向などの詳細に関してはホームページを参照頂きたい。
（http://bukai.pharm.or.jp/bukai_kozo/index.html）

編集後記
日本薬学会構造活性相関部会誌 SAR News 第 17 号をお届けいたします。ご多忙の中、ご執筆いただきました諸先生方に心よりお礼申し上げます。今回は、薬物分子がその活性を現す最重要局面である、標的生体高分子との結合の強度や構造を解析する上で不可欠な方法論である分子動力学に焦点を当てました。Perspective/Retrospective では、藤谷秀章先生（富士通研究所）に Jarzynski 等式に基づくΔG 計算法について、Cutting Edge では、合田浩明先生・広野修一先生（北里大学薬学部）に MM-PBSA 法に基づくΔG 計算法をご解説いただきました。Innovative な理論に裏打ちされたこのような分子計算によって、絶対値として化学的精度が達成されつつあることは刮目に値します。また、米澤康滋先生（大阪大学蛋白質研究所）には、古典力場を越えて量子力学との融合をなす方法論である QM/MM 分子動力学のご紹介をいた
だきました。いずれも目の離せない重要なご研究です。この SAR News が、構造活性相関研究の先端情報と展望を会員の皆様にご提供できることを編集委員一同願っております。今後ともよろしくお願い申し上げます。（編集委員会）
SAR News No.17 平成 21 年 10 月 1 日
発行：日本薬学会 構造活性相関部会長 石黒 正路
SAR News 編集委員会（委員長）久保寺 英夫
藤原 巌
黒木 保久
福島 千晶
粕谷 敦
＊本誌の全ての記事、図表等の無断複写・転載を禁じます。

-28-