Prediction of cyclohexane-water distribution coefficients for the SAMPL5 data set using molecular dynamics simulations with the OPLS-AA force field

Abstract

International audienceAll-atom molecular dynamics (MD) simulations were used to predict water-cyclohexane distribution coefficients D cw of a range of small molecules as part of the SAMPL5 blind prediction challenge. Molecules were parameterized with the trans-ferable all-atom OPLS-AA force field, which required the derivation of new parameters for sulfamides and heterocycles and validation of cyclohexane parameters as a solvent. The distribution coefficient was calculated from the solvation free energies of the compound in water and cyclohexane. Absolute solvation free energies were computed by an established protocol using windowed alchemical free energy perturbation with thermodynamic integration. This protocol resulted in an overall root mean square error (RMSE) in log D cw of almost 4 log units and an overall signed error of −3 compared to experimental data. There was no substantial overall difference in accuracy between simulating in NV T and NPT ensembles. The signed error suggests a systematic error but the experimental D cw data on their own are insufficient to Manuscript Click here to download Manuscript sampl5-manuscript.pdf Click here to view linked References 2 Ian M. Kenney et al. uncover the source of this error. Preliminary work suggests that the major source of error lies in the hydration free energy calculations