|
CONFLEX™フレキシブルな分子の配座空間を探索し、化学的に重要な配座異性体の最適化構造をもれなく見つけだします。今までの構造最適化プログラムでは、ユーザーが入力した初期構造に依存した、局所的な最適化構造しか求めることができませんでした。
CONFLEXはいくつかのアルゴリズムを用いて、最安定構造を明らかにするために、配座空間を余すところなく探索することができます。
- 貯水池注水アルゴリズム
- 可変サーチリミット
- ３モード変形操作 (corner flap, edge flip, stepwise rotation）
- 配座チェック
CONFLEXの計算結果は各種フォーマットに出力でき、ab-initio計算の初期座標として利用したり、BARISTA™のような表示ソフトにて使用できます。
創薬研究、材料開発などには、速くて、正確な配座解析が必要とされます。CONFLEXはこの問題を解決するために研究開発担当者の力強いツールとして威力を発揮します。
- 網羅的な配座探索
- 高速かつ高い精度
- 大きな分子に対応
- 並列計算オプション
- スタンドアローン及びBARISTAインタフェースによる計算
- 構造最適化
- 配座異性体
- エネルギー値リスト
CONFLEXは、官庁、大学、化学及び製薬会社の研究部門で使用されています。
CONFLEX Concepts
CONFLEX探索アルゴリズム
１．配座探索プログラム
ＣＯＮＦＬＥＸは、適当な分子力場のポテンシャル空間中に存在する、予め特定できないエネルギー極小点、すなわち配座異性体を発生するプログラムです。一般に、プログラムを用いて配座を発生する場合、出発構造の座標発生、構造最適化、そしてユニークな配座異性体であるかどうかのチェックといった、ある決まった手順（アルゴリズム）に従って配座空間を探索します。既往のプログラムの配座探索アルゴリズムは系統的手法、及びランダム法に大別できます。
ＣＯＮＦＬＥＸにおける配座空間の探索法は、貯水池（Ｒｅｓｅｒｖｏｉｒ）アルゴリズムと呼ぶ、独自の方法です。出発構造の発生方法として用いたＣｏｒｎｅｒ Ｆｌａｐに加え、Ｅｄｇｅ Ｆｌｉｐ及びＳｔｅｐｗｉｓｅ Ｒｏｔａｔｉｏｎを導入することによって、様々な分子構造（環式分子、非環式分子）に対応しました。また、新たに改良を加えたプレチェックや二面格差RMS比較法によって、効率的な配座空間の探索を実現しました。
２．貯水池(Reservoir)アルゴリズムのコンセプト
貯水池アルゴリズムは、空の貯水池に水が注ぐがごとく素早く最安定点に到達、満ちるがごとく重要な領域を調べ尽くしていきます。そして、上限の高さまで水位が到達すると、探索を終了します。例えば図に示したように、任意の安定配座(１) から探索を開始したとしましょう。この配座を基点とし、その周辺の空間にあるエネルギー極小配座（線で結ばれた配座）を調べることによって、水が流れ落ちていく方向を見きわめます。水の流れは最も低い方へ向かうはずであるから、見つかった配座の中から最もエネルギーの低いものを次の基点（２）とします。このように、基点周辺を綿密に調べ、低い方向へ移動する過程を繰り返すことによって、素早く最安定配座（３）に到達することができます。探索が最低地点に達すると貯水池の水位は徐々に上昇し始めると同時に、今度は基点となる配座をエネルギーの高いものへ移していきます。つまり、ここまでに得られた極小配座の中から、まだ基点になっていない最もエネルギーの低い配座（４）を選び、同様な操作を繰り返して行くことになります。そして、徐々に高エネルギーの配座を基点として用い（例えば７，８，．．．．．）、やがて予め設定しておいた上限（例えばLimit1の最安定配座から５kcal／mol）まで基点が到達した時点で、探索が完了したとみなします。
以上のようなアルゴリズムによる配座探索の最大の特徴は、（1）最安定配座を素早く、そして確実に発生することができ、(２)低エネルギー領域に存在する配座から集中的に発生するため、重要な配座異性体を見逃さず発生できます。また、（3）上限の高さ（基点として選択可能な配座のエネルギー最大値）を適切に調節することによって、探索する領域をある程度限定的できるため、既往の配座探索プログラムよりも計算時間が短くなります。
３．Edge FlipとStepwise Rotationの導入
他の配座探索プログラムに用いられているアルゴリズムの多くは、基点から配座を一つだけ発生するか、あるいは離れた配座空間に存在する配座、すなわち基点とは全く異なる構造を持つ配座を発生します。一方、貯水池アルゴリズムは、一つの、基点から複数の配座を発生させるが、基本的にこれらの配座は基点の配座と近い構造を持っています。つまり、空間的に近い（構造的がよく似ている）配座を発生することを意味します。基点を入れ換えながら配座空間を歩き回ることは、配座群の間のネットワークを考慮しているといえましょう。
基点周辺にある配座が発生する際、基点構造に対し部分的に小さな摂動を加え、適当な力場を用いて構造最適化を行うことにしました。局所的な摂動は、安定な配座が持っている安定な部分を乱さないため、素早く極小点に落ち、依然として安定な（空間的に近い）領域にとどまる確立が高い配座を発生すると考えられます。以前発表したＣＯＮＦＬＥＸは、この局所的な摂動として、Ｃｏｒｎｅｒ Ｆｌａｐを用いました。この操作は、我々が分子模型でしばしば行う、シクロヘキサンのいす形から舟形の配座交換をまねた操作です。その後研究が進むにつれ、この操作は次のような２つの問題点を持っていました。一つは、ｆｌａｐ原子を挟むねじれ角がＧａｕｃｈｅ+、-と続いた時に最も効果的な操作であって、同符号のＧａｕｃｈｅが続いた場合には原理的に全く変形を加えることができないということです。さらに、大きな環式分子に適用した場合、一方のねじれ角にＡｎｔｉが存在すると、予測していなかった無理な変形を行うことがあるということも問題でした。これは、本来環の外側に向けられるべき操作が、内側に向けて大きくＦｌａｐされることに起因していました。したがって、Ｃｏｒｎｅｒ Ｆｌａｐはより大きな環の配座空間を探索するには少し単純すぎて、配座発生能力に限界が生じてくることが判ってきました。
そこで、より大きな環式分子に対応した新たな摂動として、Ｅｄｇｅ Ｆｌｉｐを考案しました。この操作は、注目している環を構成する一本の結合のねじれ角がＧａｕｃｈｅであれば異符号のＧａｕｃｈｅとＡｎｔｉへ、また、中央の結合がＡｎｔｉであれば二つのＧａｕｃｈｅになるように局所的に摂動を加えます。実際にはＧａｕｈｅからＧａｕｃｈｅへ、およびＡｎｔｉからＧａｕｃｈｅへの摂動は、シクロヘキサンのいす型からねじれ船型への変換をまねた操作、いわばＤｏｕｂｌｅ Ｆｌａｐすることによって実現でき、そしてＧａｕｃｈｅからＡｎｔｉへは、そのＤｏｕｂｌｅ Ｆｌａｐを環の内側に向けて小さく行います。ここで必要なのは中央の結合がＡｎｔｉ（ＡｎｔｉＰｅｒｉｐｌａｎｅｒ）になれば良いのです。また、このとき、Ｆｌｉｐする結合とそれを挟む結合の３本のねじれ角パターンにしたがって，Ｆｌａｐする大きさや方向を調節することによって、大きく不安定化しない構造を発生できます。
以前発表したＣＯＮＦＬＥＸにおいて、側鎖に対する回転操作は、Ｆｌａｐの補助的な
意味しかもっていなかったため、非環式分子や長い側鎖を持つ分子に対して、十分に対応できませんでした。そこでＳｔｅｐｗｉｓｅ Ｒｏｔａｔｉｏｎを導入し、従来系統的探索法で用いる１２０度や６０度毎の結合強制回転を、Ｆｌａｐ／Ｆｌａｐと同じ意味、すなわち一つの局所的な摂動として用いることにしました。この操作は、結合の二面角を単に段階的に回転するだけであるが、確実に基点の配座とは異なる配座を発生し、また回転する結合の二面角以外の内部座標を変化させないので、最適化が速やかに完了します。
４．プレチェックと配座異性体のチェック
貯水池アルゴリズムと３種類の摂動の組み合わせによる出発構造の発生は、明らかに無駄な最適化ステップを含みます。つまり、基点配座Aのある部分に摂動を加え、構造最適化して得られた配座Ｂが基点となったとき，Ｂの同じ部分に同じ摂動を加えるとＡが再び発生する可能性が高くなります。このような配座変換の可逆性を確認すること自体は重要な意味を持つが、これによる時間的損失は非常に大きくなります。この無駄を取り除くため、構造最適化途中の構造をすでに保存した配座と比較し、もし同一であるとみなされればその時点で最適化を打ち切ることにしました。実際には、最適化繰り返し回数が、１０，２５，５０，１００、その後１００回毎に比較を行います。プレチェックとよぶこの手法を用いて、全計算時間を３０〜５０％ほども短縮することに成功しました。
プレチェックを含め、重複配座を比較するための方法として、二面角差ＲＭＳ値比較法も新たに開発しました。２つの配座間の対応するねじれ角の差を求め、これらの自乗平均を求めることにしました。ここで、Nは分子を構成するねじれ角の数であり、f A とｆ Bは配座ＡとＢのｉ番目のねじれ角を示します。この値が１０度以下のものを重複配座と認定することにしました。
プレチェックではこの比較のみを行うが、保存するユニークな配座異性体の認定基準としては、さらに高エネルギー配座の廃棄や、キラル中心や二重結合に対する絶対配置の検査も行います。
CONFLEXは豊橋技術科学大学の後藤仁志助教授により開発された配座創出プログラムです。
参考文献：J. Am. Chem. Soc., 1989, 111, 8950-8951.
J. Chem. Soc., Perkin Trans. 2, 1993, 187-198
CONFLEX5の一部は文部科学省「革新的技術開発研究費補助金」を受けて開発されています。 |
