CHARMM
CHARMM(チャーム、Chemistry at Harvard Macromolecular Mechanics)は、多粒子系で広く応用されている分子動力学法のための力場の名称であり、それらと関連した分子動力学シミュレーションならびに解析パッケージの名称である[1][2][3]。CHARMM開発プロジェクトにはマーティン・カープラスと彼のハーバード大学の研究グループと共に世界中の開発者のネットワークがCHARMMプログラムの開発およびメンテナンスに関わっている。このソフトウェアのライセンスは学術界で研究している人物およびグループに対して有償で提供されている。ソースコード (charmm) は学界 、政府関係、非営利の研究室には無償で提供されている。企業向けにはCHARMmと呼ばれる商用版がある。
開発元 | マーティン・カープラス、Accelrys |
---|---|
初版 | 1983年 |
最新版 |
c40b1
/ 2015年 |
最新評価版 |
c41a1
/ 2015年 |
プログラミング 言語 | FORTRAN 77/95 |
対応OS | Unix-like |
種別 | 分子動力学法 |
ライセンス | プロプライエタリ |
公式サイト | charmm.org |
CHARMM力場
編集元々、CHARMM力場は生体分子の研究にフォーカスして開発された。ここでいう生体分子とは、ペプチド、タンパク質、補欠分子族、小分子配位子、核酸、脂質、および炭水化物などの溶液、結晶、そして膜状態にあるものをいう。
タンパク質のためのCHARMM力場には以下のものがある: 融合原子 (united-atom) CHARMM19[4]、全原子 (all-atom) CHARMM22[5]、その2面角ポテンシャル補正版CHARMM22/CMAP[6]。CHARMM22タンパク質力場では、原子の部分電荷はモデル化合物と水との間の相互作用の量子化学的計算から導かれている。さらに、CHARMM22は明示的なTIP3P水模型(陽溶媒)についてパラメータ化されている。にもかかわらず、CHARMM22力場は陰溶媒を用いて頻繁に使用されている。2006年、CHARMM22/CMAPの特別版が陰溶媒GBSWを用いた矛盾のない使用のために再パラメータ化された[7]。
DNA、RNA、脂質については、CHARMM27[8]が使用される。一部の力場は組み合わせることができる。例えばタンパク質-DNA結合のシミュレーションのためにはCHARMM22とCHARMM27を組み合わせて使用する。加えて、NAD+、糖、フッ素化化合物等についてのパラメータもダウンロードできる。これらの力場のバージョン数はそれらが最初に現われたCHARMMのバージョンを示しているが、当然CHARMM実行プログラムの後続のバージョンと共に使用することができる。同様に、CHARMM力場群はそれらをサポートするその他の分子動力学プログラム内で使用することができる。
2009年、薬らしい(ドラッグライクな)分子のための一般力場(CGenFF)が発表された。CGenFFは、「多数の複素環骨格を含む生体分子ならびにドラッグライクな分子中の存在する幅広い化学基を扱う」[9]。CGenFFは化学基のいかなる組合せも扱うことができるように設計されている。これは不可避的に特定の種類の分子を表現する際の精度の低下を伴う。開発者のMacKerellのウェブサイトにおいて、使用者が特化した力場が既に存在する分子(タンパク質、核酸等)についてCGenffパラメータを使わないよう繰り返し警告されている。
CHARMMは2つの手法を用いた分極可能力場も含む。1つは揺らぎ電荷(fluctuating charge、FQ)モデルに基づいている。このモデルは電荷平衡(Charge Equilibration、CHEQ)とも呼ばれる[10][11]。もう1つはドルーデ殻(分散振動)モデルに基づいている[12][13]。
これらの力場全てについてのパラメータは無償でMacKerellのウェブサイト[14]からダウンロードできる。
CHARMM分子動力学プログラム
編集CHARMMプログラムでは、幅広い分子シミュレーションの生成と解析が可能である。最も基本的な種類のシミュレーションは与えられた構造のエネルギー最小化と分子動力学トラジェクトリのプロダクションランである。
より高度な機能としては、コンフォメーションおよびパスサンプリング法、自由エネルギー摂動法(FEP)、擬調和エントロピーの見積もり、相関分析、QM/MM法がある。
CHARMMは最も古い分子動力学プログラムの一つである。CHARMMは数多くの機能を集積しており、それらの一部はわずかな違いがある複数のキーワード下で重複している。これは、世界中で多くのグループがCHARMMに関わっていることの不可避的な結果である。更新履歴とCHARMMのソースコードに主要な開発者らの名前と所属を見ることができる。ミシガン大学のチャールズ・L・ブルックス3世のグループの貢献が顕著である。 CHARMMは非常に多くの計算ツール群を提供している。例えば、コンフォメーションおよびパスサンプリング法、自由エネルギーの見積もり、分子エネルギーの最小化、動力学、分析テクニック、そして分子モデル構築などができるようになっている。
CHARMMプログラムはシリアル構成および並列構成ともに非常に多くのプラットホームに移植されている。開発者のマーティン・カープラスは、2013年ノーベル化学賞を受賞した。
プログラムの歴史
編集1969年頃、小分子のためのポテンシャルエネルギー関数の開発に大きな関心が持たれていた。CHARMMはハーバード大学のマーティン・カープラスのグループに起源がある。カープラスと彼の当時の大学院生Bruce Gelinは任意のアミノ酸配列と一連の座標(例えばX線構造から)を扱うことができ、この情報を原子の位置の関数としての系のエネルギーを計算するために使うことを可能とするプログラムの開発の機が熟したと判断した。カープラスは
- ワイツマン科学研究所のSchneior Lifsonのグループ、特にハーバード大学にポスドクとして赴き、自身のconsistent force field(CFF)プログラムをもたらしたアリー・ウォーシェル
- コーネル大学のHarold Scheragaのグループ
- タンパク質についてのマイケル・レヴィットの先駆的なエネルギー計算
を含むCHARMMプログラムの開発における主要なアイデアの提供の重要性を認めている。
1980年代、最終的に論文が発表され、CHARMMが公開された。その時までにはGelinのプログラムは大幅に再構築されていた。発表のため、Bob BruccoleriはHARMM(HARvard Macromolecular Mechanics)という名称を思い付いたが、これは適切ではないように見えた。そこで、彼らは先頭に化学(Chemistry)のCを加えた[15]。
2016年現在の最新版はc40b2(2015年)である[16]。
脚注
編集- ^ Brooks BR, Bruccoleri RE, Olafson BD, States DJ, Swaminathan S, Karplus M (1983). “CHARMM: A program for macromolecular energy, minimization, and dynamics calculations”. J Comp Chem 4 (2): 187–217. doi:10.1002/jcc.540040211.
- ^ MacKerell, A.D., Jr.; Brooks, B.; Brooks, C. L., III; Nilsson, L.; Roux, B.; Won, Y.; Karplus, M. (1998). "CHARMM: The Energy Function and Its Parameterization with an Overview of the Program". In Schleyer, P.v.R.; et al. (eds.). The Encyclopedia of Computational Chemistry. Vol. 1. Chichester: John Wiley & Sons. pp. 271–277.
- ^ Brooks BR, Brooks CL 3rd, Mackerell AD Jr, Nilsson L, Petrella RJ, Roux B, Won Y, Archontis G, Bartels C, Boresch S, Caflisch A, Caves L, Cui Q, Dinner AR, Feig M, Fischer S, Gao J, Hodoscek M, Im W, Kuczera K, Lazaridis T, Ma J, Ovchinnikov V, Paci E, Pastor RW, Post CB, Pu JZ, Schaefer M, Tidor B, Venable RM, Woodcock HL, Wu X, Yang W, York DM, Karplus M (29 July 2009). “CHARMM: The biomolecular simulation program”. Journal of Computational Chemistry 30 (10): 1545–1614. doi:10.1002/jcc.21287. PMC 2810661. PMID 19444816 .
- ^ Reiher, III WH (1985). “Theoretical studies of hydrogen bonding”. PhD Thesis at Harvard University.
- ^ MacKerell, Jr. AD (1998). “All-atom empirical potential for molecular modeling and dynamics studies of proteins”. J Phys Chem B 102 (18): 3586–3616. doi:10.1021/jp973084f.
- ^ MacKerell, Jr. AD, Feig M, Brooks, III CL (2004). “Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations”. J Comput Chem 25 (11): 1400–1415. doi:10.1002/jcc.20065. PMID 15185334.
- ^ Brooks CL, Chen J, Im W (2006). “Balancing solvation and intramolecular interactions: toward a consistent generalized born force field (CMAP opt. for GBSW)”. J Am Chem Soc 128 (11): 3728–3736. doi:10.1021/ja057216r. PMC 2596729. PMID 16536547 .
- ^ MacKerell, Jr. AD, Banavali N, Foloppe N (2001). “Development and current status of the CHARMM force field for nucleic acids”. Biopolymers 56 (4): 257–265. doi:10.1002/1097-0282(2000)56:4<257::AID-BIP10029>3.0.CO;2-W. PMID 11754339.
- ^ Vanommeslaeghe K, Hatcher E, Acharya C, Kundu S, Zhong S, Shim J, Darian E, Guvench O, Lopes P, Vorobyov I, Mackerell AD Jr (2009). “CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields”. J Comput Chem 31 (4): 671–90. doi:10.1002/jcc.21367. PMC 2888302. PMID 19575467 .
- ^ Patel S, Brooks CL 3rd (2004). “CHARMM fluctuating charge force field for proteins: I parameterization and application to bulk organic liquid simulations”. J Comput Chem 25 (1): 1–15. doi:10.1002/jcc.10355. PMID 14634989.
- ^ Patel S, Mackerell AD Jr, Brooks CL 3rd (2004). “CHARMM fluctuating charge force field for proteins: II protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model”. J Comput Chem 25 (12): 1504–1514. doi:10.1002/jcc.20077. PMID 15224394.
- ^ Lamoureux G, Roux B (2003). “Modeling induced polarization with classical Drude oscillators: Theory and molecular dynamics simulation algorithm”. J Chem Phys 119 (6): 3025–3039. Bibcode: 2003JChPh.119.3025L. doi:10.1063/1.1589749.
- ^ Lamoureux G, Harder E, Vorobyov IV, Roux B, MacKerell AD (2006). “A polarizable model of water for molecular dynamics simulations of biomolecules”. Chem Phys Lett 418: 245–249. Bibcode: 2006CPL...418..245L. doi:10.1016/j.cplett.2005.10.135.
- ^ “CHARMM Force Field Files”. MacKerell Lab. 2016年4月5日閲覧。
- ^ Karplus M (2006). “Spinach on the ceiling: a theoretical chemist's return to biology”. Annu Rev Biophys Biomol Struct 35 (1): 1–47. doi:10.1146/annurev.biophys.33.110502.133350. PMID 16689626.
- ^ “Versions”. CHARMM 2016. 2016年4月5日閲覧。