Predictions of the physicochemical properties of amino acid side chain analogs using molecular simulation

A candidate drug compound is released for clinical trails (in vivo activity) only if its physicochemical properties meet desirable bioavailability and partitioning criteria. Amino acid side chain analogs play vital role in the functionalities of protein and peptides and as such are important in drug discovery. We demonstrate here that the predictions of solvation free energies in water, in 1-octanol, and self-solvation free energies computed using force field-based expanded ensemble molecular dynamics simulation provide good accuracy compared to existing empirical and semi-empirical methods. These solvation free energies are then, as shown here, used for the prediction of a wide range of physicochemical properties important in the assessment of bioavailability and partitioning of compounds. In particular, we consider here the vapor pressure, the solubility in both water and 1-octanol, and the air-water, air-octanol, and octanol-water partition coefficients of amino acid side chain analogs computed from the solvation free energies. The calculated solvation free energies using different force fields are compared against each other and with available experimental data. The protocol here can also be used for a newly designed drug and other molecules where force field parameters and charges are obtained from density functional theory.