Combining MOSCED with molecular simulation free energy calculations or electronic structure calculations to develop an efficient tool for solvent formulation and selection.

Solubility parameter based methods have long been a valuable tool for solvent formulation and selection. Of these methods, the MOdified Separation of Cohesive Energy Density (MOSCED) has recently been shown to correlate well the equilibrium solubility of multifunctional non-electrolyte solids. However, before it can be applied to a novel solute, a limited amount of reference solubility data is required to regress the necessary MOSCED parameters. Here we demonstrate for the solutes methylparaben, ethylparaben, propylparaben, butylparaben, lidocaine and ephedrine how conventional molecular simulation free energy calculations or electronic structure calculations in a continuum solvent, here the SMD or SM8 solvation model, can instead be used to generate the necessary reference data, resulting in a predictive flavor of MOSCED. Adopting the melting point temperature and enthalpy of fusion of these compounds from experiment, we are able to predict equilibrium solubilities. We find the method is able to well correlate the (mole fraction) equilibrium solubility in non-aqueous solvents over four orders of magnitude with good quantitative agreement.

Predicting cyclohexane/water distribution coefficients for the SAMPL5 challenge using MOSCED and the SMD solvation model.

We present blind predictions using the solubility parameter based method MOSCED submitted for the SAMPL5 challenge on calculating cyclohexane/water distribution coefficients at 298 K. Reference data to parameterize MOSCED was generated with knowledge only of chemical structure by performing solvation free energy calculations using electronic structure calculations in the SMD continuum solvent. To maintain simplicity and use only a single method, we approximate the distribution coefficient with the partition coefficient of the neutral species. Over the final SAMPL5 set of 53 compounds, we achieved an average unsigned error of [Formula: see text] log units (ranking 15 out of 62 entries), the correlation coefficient (R) was [Formula: see text] (ranking 35), and [Formula: see text] of the predictions had the correct sign (ranking 30). While used here to predict cyclohexane/water distribution coefficients at 298 K, MOSCED is broadly applicable, allowing one to predict temperature dependent infinite dilution activity coefficients in any solvent for which parameters exist, and provides a means by which an excess Gibbs free energy model may be parameterized to predict composition dependent phase-equilibrium.