用途
為有機小分子添加原子電荷,並檢查是否存在Amber GAFF力場所沒有的拓撲參數。本工具提供兩種計算電荷的方法——RESP擬合靜電勢和AM1-BCC。
預備知識
原子電荷
原子電荷是用點電荷形式描述化學體系電荷分布的方法,常見電荷模型有:Mulliken、MK、MMFF94、AM1-BCC、Gasteiger、QEq;等等。涉及分子間相互作用的分子模擬計算(比如:分子對接、分子動力學模擬)都要求賦予原子電荷。採用合適的電荷模型對於計算的準確性是極為重要的。關於各類原子電荷計算方法的對比,可閱讀盧天的《原子電荷計算方法的對比》[1]。
擬合靜電勢電荷
在各種原子電荷計算方法中,擬合靜電勢電荷是極為重要的一類方法。其中,最為著名的要數Merz-Kollman (MK)、CHELPG和RESP。前兩種方法非常適合於模擬剛性分子,而後者則是最適合用於柔性小分子的分子模擬。
RESPRESP,全稱Restrained ElectroStatic Potential,是由Kollman等人[2]在1993年提出的計算原子電荷的方法。在Amber分子動力學模擬中,RESP是有機小分子計算電荷的標配方法。它需要配合使用量化計算軟體(如Gaussian、ORCA)獲得靜電勢,然後進行限制性擬合計算得到。
AM1-BCC與之比較接近的電荷是2000-2002年Bayly等人[3-4]提出的AM1-BCC電荷,在Amber中自帶計算程序。它首先在AM1水平下優化結構和計算靜電勢擬合得到電荷,然後進行鍵電荷校正(BCC)。AM1-BCC的計算開銷比RESP要小,準確性不比RESP更好,但適合於處理大批量的分子,或者沒有Gaussian等商業軟體版權的情況下使用。(不過,可以使用學術免費的ORCA程序,請看【量化計算(ORCA 4.2)】教程)。
淨電荷
分子結構中不能抵消的正電荷或負電荷稱為淨電荷(net charge),即各原子正負電荷的代數和。通常情況下,淨電荷總是整數。中性分子的淨電荷為0,帶電分子的淨電荷即為其含有的價電荷,如:乙酸淨電荷為0,去質子化(脫去羧基氫)後為-1,十八烷基胺的淨電荷為0,質子化後為1。
入口
平臺地址:https://cloud.yinfotek.com/
功能入口:左側菜單欄【計算方案】->【小工具】->【分子動力學】->【添加原子電荷】
步驟
1. 上傳分子結構文件
分子結構需為良好的3D結構,無氫原子缺失。可使用【準備化合物】小工具進行處理。
2. 設置淨電荷
該項默認為空,即由程序自動計算淨電荷。但對於特殊結構,可能會計算出錯。必要時,請用戶根據自己的化學知識判斷填寫。
3. 選擇電荷計算方法
RESP該方法需要用戶上傳包含靜電勢的量化計算文件(可使用【量化計算】方案進行計算),並指定該文件的計算程序。
AM1-BCC該方法無需提供額外數據文件,提交後將調用AmberTools 20的sqm半經驗程序即時計算。對於較大的柔性分子,計算時間較長,可能因超出限制(10 min)而失敗。同時,AM1-BCC的精度不如RESP電荷。因此,若非特殊情況,強烈建議使用RESP。
4. 點擊【生成】,等待返回結果
本工具除了計算原子電荷,還檢查該分子是否存在Amber GAFF力場所沒有的參數。該檢查主要用於使用amber力場的分子動力學模擬,若進行其他不需要該文件的分子模擬,請忽略之。
GAFF力場已經發展得相當好,幾乎覆蓋了大多數的有機小分子結構,不太容易出現參數缺失的情況。假若真的出現,請下載ligand.frcmod文件自行填補參數,在提交分子動力學計算時上傳該文件。
5. 查看電荷分布
從分子視圖上可查看生成的原子電荷,尤其注意極性原子及其附近的電荷分布是否合理。
6. 下載分子結構文件。
參考文獻
[1] 盧天, 陳飛武. 原子電荷計算方法的對比[J]. 物理化學學報, 2012, 28(1):1-18. doi: 10.3866/PKU.WHXB2012281
[2] Christopher I. Bayly, Piotr Cieplak, Wendy Cornell, and Peter A. Kollman. A well-behaved electrostatic potential based method using charge restraints for deriving atomic charges: the RESP model. The Journal of Physical Chemistry 1993, 97, 40, 10269-10280. doi: 10.1021/j100142a004
[3] Jakalian, A., Bush, B.L., Jack, D.B. and Bayly, C.I. (2000), Fast, efficient generation of highquality atomic charges. AM1BCC model: I. Method. J. Comput. Chem., 21: 132-146. https://doi.org/10.1002/(SICI)1096-987X(20000130)21:2%3C132::AID-JCC5%3E3.0.CO;2-P
[4] Jakalian, A., Jack, D.B. and Bayly, C.I. (2002), Fast, efficient generation of highquality atomic charges. AM1BCC model: II. Parameterization and validation. J. Comput. Chem., 23: 1623-1641. https://doi.org/10.1002/jcc.1012