Force Field Benchmark of Amino Acids. 2. Partition Coefficients Between Water and Organic Solvents
Haiyang Zhang,Yang Jiang,Ziheng Cui,Chunhua Yin
DOI: https://doi.org/10.1021/acs.jcim.8b00493
IF: 6.162
2018-01-01
Journal of Chemical Information and Modeling
Abstract:The partitioning of amino acids between water and apolar environments is of vital importance in protein function and drug delivery. Here we present an extensive benchmark for octanol/water (log P-oct), chloroform/water (log P-clf), and cyclohexane/water (log P-chr) partition coefficients of neutral amino acid side chain analogues (SCAs) with Amber families of ff99SB-ILDN, ff03, ff14SB, fb15, and ff15ipq, CHARMM 27, GROMOS 53A6, and OPLS-AA/L force fields. A root-mean-square error (RMSE) of 0.4-1.3 log units from experiment is observed for the tested FFs, of which Amber ff94 lineages of ff99SB-ILDN, ff14SB, and fb15 perform best with an RMSE and mean signed error (MSE) of about 0.5 and 0.2 log units, respectively, a performance comparable with quantum mechanical SMD calculations. This finding retains the possibility of modeling proteins in varied environments with one set of classical molecular mechanical force fields. All the FFs tend to overestimate log P, except for GROMOS 53A6 underestimating log P-clf and log P-chx These discrepancies are mainly due to the larger overestimated solvation free energies in water (Delta G(wat)) relative to that in organic solvents (Delta G(oct), Delta Gd-clf and Delta G(chx)); for GROMOS 53A6, it is due to the underestimated Delta G(wat) and Delta G(oct). The latest water models of "FB" and "OPC" families paired with the recent Amber fb15 do not show an obvious improvement for Delta G(wat) and log P calculations. The van der Waals interaction between amino acids and cyclohexane is found to be too strong (overestimated) systematically. Scaling protein water interactions lead to more favorable Delta G(wat), thereby lowering log P and resulting in a better performance for Amber ff03ws, while such scaling seems a bit too much for Amber ff99SBws. This, along with our previous work (Zhang et al. J. Chem. Inf. Model. 2018, 58, 1037-1052), may aid in the development and systematic improvements of classical force fields to model proteins in aqueous and nonaqueous phases accurately.