2009-2010年
- マルチコアCPU搭載マシンにおけるマルチスレッド化
- Half-Harmonic関数とMorse関数を使ったPseudo Force計算を導入
- Imai-Karplus式によるNMRの結合次数計算法（3JHHM）を導入
- 結晶計算の実体アンサンブルとなる単位格子を拡張するためのオプションを追加
2007-2008年
- 溶媒効果を利用したLogP計算
- 結晶構造計算の改良
- 力場のねじれ角パラメーターの自動生成機能
- ユーザーによる原子タイプの指定機能
2006-2007年
- 溶媒効果を取り入れた配座探索
- 結晶構造最適化
- アミノ酸残基置換機能の追加
- フロンティアモード追跡(FMF)法(配座変換遷移状態探索)
2005-2006年
- 偽ポテンシャル構造固定法
- ペプチド結合ねじれ角判定
- 動的基準振動解析法
- MMFF94s/GB/SA(溶媒効果)
2004-2005年
- MMFF94/MMFF94s
- PDBフォーマットRead/Write
- 配座クラスタリング
- UV/Vis/CDスペクトル計算
CONFLEX に含まれる力場の中で、 MMFF94 および MMFF94s では、分子に含まれるほとんどの原子上に、『結合電荷増分法』という手法を用いて電荷を置くようにしています。これらの電荷は、小さい分子での高精度 ab initio 分子軌道計算で求められた電荷分布を再現するようにパラメータが決められておりますが、計算の初期段階で電荷分布を決定した後は電荷の値を固定したまま最適化を行っています。
一方で MMFF/NQEq では、上記手法の代わりに『電荷平衡法(QEq 法)』という手法を用いて電荷および静電相互作用エネルギーを計算しております（QEq 法の概要は下記参照）。
QEq 法では、電荷分布が構造に依存するため構造変化するたびに電荷を計算し直さなければならないのですが、計算自体は簡便なため計算時間の増加は非常に少なく、かつパラメータは我々が独自に最適化したものを用いておりますので、高精度の静電相互作用計算を他の力場と同じくらいの計算時間で可能としております。ただし現バージョンの CONFLEX におきましては、この力場を用いた計算も可能ではございますが、他の相互作用を計算する際に用いるパラメータが MMFF で用いる値のままになっておりまして、最適なものにはなっておりません。現在必要なパラメータの最適化を行っている最中でありまして、決定したものから逐次導入しているところでございますので、ご使用の際はこの点をご理解ください。
電荷平衡法(QEq 法)：
静電相互作用エネルギーを、原子が電子を引き込む力(１中心項)と電荷を帯びた原子間に働くクーロン力(2中心項)の和として表し、各原子の電気陰性度を平衡化させることで、電荷分布およびエネルギーを求める方法。1991年に Rappe と Goddard により提唱された。
詳細はシーエムシー出版 機能材料１１月号(2005)
「溶媒効果を取り入れた新規電荷平衡法(NQEq)の開発と分子力場への応用」
中山尚史、長嶋雲兵、後藤仁志
をご覧下さい。
|
