Research article | Open | Published:
Predicting Abraham model solvent coefficients
Chemistry Central Journalvolume 9, Article number: 12 (2015)
The Abraham general solvation model can be used in a broad set of scenarios involving partitioning and solubility, yet is limited to a set of solvents with measured Abraham coefficients. Here we extend the range of applicability of Abraham’s model by creating open models that can be used to predict the solvent coefficients for all organic solvents.
We created open random forest models for the solvent coefficients e, s, a, b, and v that had out-of-bag R2 values of 0.31, 0.77, 0.92, 0.47, and 0.63 respectively. The models were used to suggest sustainable solvent replacements for commonly used solvents. For example, our models predict that propylene glycol may be used as a general sustainable solvent replacement for methanol.
The solvent coefficient models extend the range of applicability of the Abraham general solvation equations to all organic solvents. The models were developed under Open Notebook Science conditions which makes them open, reproducible, and as useful as possible.
The Abraham model was developed and is widely used to predict partition coefficients for both conventional organic solvents [1-11] and ionic liquid solvents [12,13], for the partitioning of drug molecules between blood and select body organs [14-18], and for partitioning into micelles  and for prediction of enthalpies of solvation in organic solvents  and ionic organic liquids . The Abraham model is based on the linear free energy relationship (LFER)
where logP is the solvent/water partition coefficient. Under reasonable conditions, this model can also be used to predict the solubility of organic compounds in organic solvents  as follows
where S s is the molar concentration of the solute in the organic solvent, S w is the molar concentration of the solute in water, (c, e, s, a, b) are the solvent coefficients, and (E, S, A, B, V) are the solute descriptors: E is the solute excess molar refractivity in units of (cm^3/mol)/10, S is the solute dipolarity/polarizability, A and B are the overall or summation hydrogen bond acidity and basicity, and V is the McGowan characteristic volume in units of (cm^3/mol)/100.
The solvent coefficients are obtained by linear regression using experimentally determined partitions and solubilities of solutes with known Abraham descriptors. Traditionally, the intercept c is allowed to float and is assumed to encode information not characterized by the other solvent-solute interaction terms. However, for some partitioning systems the value of c can vary greatly depending upon the training-set used . This makes it difficult to directly compare different solvents by examining their solvent coefficients. Van Noort has even suggested that the c-coefficient be derived directly from structure before the other coefficients are determined . A problem with this suggestion is that the c-coefficient depends on the standard state. Partition coefficients can be expressed in concentration units of molarity and mole fractions, and the numerical value of the c-coefficient will be different for each concentration unit. Abraham model correlations considered in this study have partition coefficients expressed in concentration units of molarity.
To date, solvent coefficients have been determined for over 90 commonly used solvents (Additional file 1), and group contribution methods have been developed to approximate all coefficients for certain classes of solvents that do not have published solvent coefficients [25,26]. The solvent coefficients in the supporting material pertain to dry solvents, or solvents that take up very little water (hexane, toluene, etc.). This study expands the applicability of the Abraham model by developing open models, using open descriptors from the Chemistry Development Kit (CDK)  that can be used to predict the Abraham solvent coefficients of any organic solvent directly from structure.
In order to directly compare various solvents, it is advantageous to first recalculate the solvent coefficients with the c-coefficient equal zero. This was accomplished by using equation (1) to calculate the log P values for 2144 compounds from our Open Data database of compounds with known Abraham descriptors  and then by regressing the results against the following equation
where the subscript-zero indicates that c = 0 has been used in the regression . As an informational note one could have set the c-coefficient of a given solvent equal to a calculated average value determined from numerical c-coefficients of solvents similar to the solvent under consideration. For example, the c-coefficient of all alkane solvents could be set equal to c = 0.225, which is the average value for the c-coefficients of the 13 alkane and cycloalkane solvents for which log P correlations have been determined. While average values could be used for several solvents, there is the problem of what value to use in the case of solvents for which a similar solvent log P solvent is not available. Abraham model correlations are available for two dialkyl ethers (e.g., diethyl ether and dibutyl ether) and for several alcohols, but not for alkoxyalcohols (e.g., 2-ethoxyethanol, 2-propoxyethanol, 2-butyoxyethanol) which contain both an ether and hydroxyl alcohol group. Our intended solvent set in the present communication includes the alternative “green” solvents, and there a number of solvents in this group that contain multi-functional groups. For several of the solvents on the list of alternative “green” solvents, such as 1,3-dioxan-5-ol, 1,3-dioxolane-4-methanol, 3-hydroxypropionic acid, 5-(hydroxymethyl)furfural, ethyl lactate, furfuryl alcohol, and other solvents, there are no similar solvents having a Abraham model log P correlation. To treat all solvents equally we have elected to set c = 0 in this study.
Table 1 lists the original solvent coefficients together with the c = 0 adjusted coefficients. Comparing the coefficients, we see, not surprisingly, the largest changes in coefficient values occur for solvents with c-values furthest away from zero (Additional file 1). What is intriguing is that all the coefficients move consistently the same way. That is, solvents with negative c-values all saw an increase in e and b (and a decrease in s, a, and v) when recalculated, whereas solvents with positive c-values all saw an increase in s, a, and v (and decrease in e and b).
One way to measure the effect of making c = 0 is to evaluate how the values of each solute-solvent term change as measured against the average solute descriptors (Eave = 0.884, Save = 1.002, Aave = 0.173, Bave = 0.486, Vave = 1.308). By multiplying the average absolute deviation of the solvent coefficients and the mean solute descriptor value, e.g. AAE(v) * Mean(Vave), the coefficients shifted from greatest to least in the following order v (0.124), s (0.043), e (0.013), b (0.011), a (0.010).
Results and discussion
We calculated CDK descriptors for each solvent using the cdkdescui  and then created five random forest models for e0, s0, a0, b0, and v0 using R. The resulting models had out of bag (OOB) R2 values ranging between the barely significant 0.31 for e0 to the very seignificant 0.92 for a0, see the Open Notebook page for more details . It is important to note that due to the limited number of data points, we decided not to split the data into training and test sets and instead use the OOB values which are automatically generated with random forest models as our means of validation. A summary of the modeling results can be found in Table 2.
Quite why some endpoints are more difficult to model than others is not known. Comparing the OOB R2 values with the standard deviation of the endpoints (e0: 0.31, s0: 0.77, a0: 0.92, b0:0.47, and v0: 0.63) we see no negative correlation between the range of a given endpoint and the actual prediction performances of the associated models as one would possibly suspect. It is our conjecture that as more measured values become available that refined models will have better performance. For now, these models should be used only as an initial starting point for exploring the wider solvent chemical space.
Errors in the predications of the coefficients for new solvents are not equivalent because when used to predict partition coefficients they are scaled by their corresponding Abraham descriptors, see equation 3. Thus, on average, when predicting solvent coefficients for new solvents, the errors in predicting v and s are more significant that errors in predicting a and b due to the difference in the sizes of average values for the solute descriptors. Multiplying the OOB-RMSE for each coefficient by the corresponding average descriptors value we see the following scaled RMSE values for e0, s0, a0, b0, and v0 of 0.16, 0.33, 0.08, 0.23, and 0.30 respectively. Thus the poor OOB R2 values for e0 (0.31) and b0 (0.47) seem not to be as detrimental to the applicability of the model as suggested by a first glance.
To analyze the modeling results further and to investigate model outliers we calculated an adjusted error D, the distance between the observed values and the predicted values scaled by the average descriptor values, for each solvent using the following equation:
where the superscript p indicates the predicted value. These distances were then plotted as colors on a graph with the x and y axes corresponding to the first two principal components of the measured values for e0, s0, a0, b0, and v0, see Figure 1. Those solvents colored red have higher calculated distances between their measured and predicted values [Figure 1].
As we can see from the figure, model outliers include: formamide, trifluoroethanol, carbon disulfide, and DMSO. These solvents are on the outskirts of the chemical space. In fact, we can clearly see that the model makes far better predictions for solvents towards the center of the chemical space with particular success in predicting the coefficients for series such as alkanes and alcohols. These observations should give us caution when using the models to predict the solvent coefficients for novel solvents, especially when they do not lie within the chemical space established by solvents with known coefficients.
These Open Models (CC0) can be downloaded from the Open Notebook pages [29,31] and can be used to predict the solvent coefficients for any organic solvent; either with the view of predicting partition coefficients or other partitioning processes including solubilities via equation (1); or with the view of finding replacement and novel solvents for current syntheses, recrystallization procedures, and other solvent dependent processes . As an informational note we remind readers that solute solubility and partitioning are only two of the considerations in finding an appropriate replacement solvent. Other considerations include the toxicity and the purchase price of the solvent, disposal costs of the solvent, physical properties of the solvent, and whether or not the solvent undergoes any undesired chemical reactions with other chemical compounds that might be present in the solution. For example, some chemical reactions take place at elevated temperatures and here one would want to use a solvent having a sufficiently high boiling point temperature that it would not vaporize under the experimental conditions.
As an example of the application of our models, we used our models to calculate the solvent descriptors for a list of sustainable solvents from a paper by Moity et. al. . The resulting coefficients for 119 select novel sustainable solvents are presented in Table 3. A complete set of coefficients for all 293 solvents (sustainable, classic, and measured) can be found in Additional file 2. These values should be used in light of the limitation of the model as described above, as possible starting places for further investigation, and not as gospel.
By comparing the predicted solvent coefficients to that of solvents with measured coefficients, we can make solvent replacement suggestions both in general and in particular. In general, the distance between solvents can be measured as the difference in predicted solubilities for the average compound.
Using this method we found several possible replacements. For example, 1,2-propylene glycol (e0 = 0.387, s0 = −0.447, a0 = 0.259, b0 = −3.447, v0 = 3.586) and methanol (e0 = 0.312, s0 = −0.649, a0 = 0.330, b0 = −3.355, v0 = 3.691) have a d-value of 0.07. This suggests that 1,2-propylene glycol may be a general sustainable solvent replacement for methanol. To confirm our model’s suggestion, we compared the solubilities of compounds from the Open Notebook Science Challenge solubility database  that had solubility values for both 1,2-propylene glycol and methanol, see [Figure 2].
Examining Figure 2, we see that solubility values are of the same order in most cases. The biggest discrepancy being for dimethyl fumerate. The measured solubility values are reported to be 0.182 M and 0.005 M for methanol and propylene glycol respectively , whereas the predicted solubilities are 0.174 M for methanol and 0.232 M for propylene glycol based upon the Abraham descriptors: E = 0.292, S = 1.511, A = 0.000, B = 0.456, V = 1.060 . This suggests that the reported value for the solubility of dimethyl fumerate in ethylene glycol may be incorrect and that, in general, 1,2-propylene glycol is a sustainable solvent replacement for methanol.
Other strongly suggested general replacements include: dimethyl adipate for hexane, ethanol/water(50:50)vol for o-dichlorobenzene, and alpha-pinene for 1,1,1-trichloroethane. Many more replacement suggestions can be generated by this technique.
In a similar manner to the above procedure for general solvent replacement for all possible solutes, one can easily compare partition and solvation properties across all solvents for a specific solute (or set of solutes) with known or predicted Abraham descriptors (E, S, A, B, V). For example, using descriptors E = 0.730, S = 0.90, A = 0.59, B = 0.40, V = 0.9317 for benzoic acid (and using d = 0.001), we can make several benzoic acid-specific solvent replacement recommendations, see Table 4. These replacement suggestions do not seem unreasonable chemically and several examples can be explicitly verified by comparing actual measured solubility values . Such a procedure can easily be done for other specific compounds with known or predicted Abraham descriptors to find alternative green solvents in varying specific circumstances (solubility, partition, etc.).
In addition to sustainable solvents, we also considered the list of commonly used solvents in the pharmaceutical industry . Of all the solvents listed, the only one not covered previously by this work (Additional file 2) was 4-methylpent-3-en-2-one which has SMILES: O = C(\\C = C(/C)C)C and predicted solvent coefficients: e0 = 0.269, s0 = −0.362, a0 = −0.610, b0 = −4.830, v0 = 4.240.
We have provided a set of Open Models that can be used to predict the Abraham coefficients for any organic solvent. These coefficients can then in turn be used to predict various partition processes and solubilities of compounds with known or predicted Abraham descriptors. We illustrated the usefulness of the models by demonstrating how one can compare solvent coefficients both in general and in particular for specific solutes or sets of solutes to find solvent replacement leads.
Linear free energy relationship
Chemistry development kit
Average absolute error
Out of bag
Simplified molecular-input line-entry system
Open Notebook Science
Sprunger LM, Achi SS, Pointer R, Acree Jr WE, Abraham MH. Development of Abraham model correlations for solvation characteristics of secondary and branched alcohols. Fluid Phase Equilibr. 2010;288:121–7.
Sprunger LM, Achi SS, Pointer R, Blake-Taylor BH, Acree Jr WE, Abraham MH. Development of Abraham model correlations for solvation characteristics of linear alcohols. Fluid Phase Equilibr. 2009;286:170–4.
Sprunger LM, Proctor A, Acree Jr WE, Abraham MH, Benjelloun-Dakhama N. Correlation and prediction of partition coefficient between the gas phase and water, and the solvents dry methyl acetate, dry and wet ethyl acetate, and dry and wet butyl acetate. Fluid Phase Equilibr. 2008;270:30–44.
Abraham MH, Acree Jr WE, Leo AJ, Hoekman D. The partition of compounds from water and from air into wet and dry ketones. New J Chem. 2009;33:568–73.
Abraham MH, Acree Jr WE, Cometto-Muniz JE. Partition of compounds from water and from air into amides. New J Chem. 2009;33:2034–43.
Abraham MH, Acree Jr WE, Leo AJ, Hoekman D. Partition of compounds from water and from air into the wet and dry monohalobenzenes. New J Chem. 2009;33:1685–92.
Abraham MH, Acree Jr WE. Comparison of solubility of gases and vapours in wet and dry alcohols, especially octan-1-ol. J Phys Org Chem. 2008;21:823–32.
Sprunger LM, Achi SS, Acree Jr WE, Abraham MH, Leo AJ, Hoekman D. Correlation and prediction of solute transfer to chloroalkanes from both water and the gas phase. Fluid Phase Equilibr. 2009;281:144–62.
Sprunger LM, Gibbs J, Acree Jr WE, Abraham MH. Correlation and prediction of partition coefficients for solute transfer to 1,2-dichloroethane from both water and from the gas phase. Fluid Phase Equilibr. 2008;273:78–86.
Abraham MH, Zissimos AM, Acree Jr WE. Partition of solutes into wet and dry ethers; an LFER analysis. New J Chem. 2003;27:1041–4.
Abraham MH, Zissimos AM, Acree Jr WE. Partition of solutes from the gas phase and from water to wet and dry di-n-butyl ether: a linear free energy relationship analysis. Phys Chem Chem Phys. 2001;3:3732–6.
Acree Jr WE, Abraham MH. The analysis of solvation in ionic liquids and organic solvents using the Abraham linear free energy relationship. J Chem Technol Biotechnol. 2006;81:1441–6.
Abraham MH, Acree Jr WE. Comparative analysis of solvation and selectivity in room temperature ionic liquids using the Abraham linear free energy relationship. Green Chem. 2006;8:906–15.
Abraham MH, Ibrahim A, Acree Jr WE. Air to liver partition coefficients for volatile organic compounds and blood to liver partition coefficients for volatile organic compounds and drugs. Eur J Med Chem. 2007;42:743–51.
Abraham MH, Ibrahim A, Acree Jr WE. Air to brain, blood to brain and plasma to brain distribution of volatile organic compounds: linear free energy analyses. Eur J Med Chem. 2008;41:494–502.
Abraham MH, Ibrahim A, Zhao Y, Acree Jr WE. A database for partition of volatile organic compounds and drugs from blood/plasma/serum to brain, and an LFER analysis of the data. J Pharm Sci. 2006;95:2091–100.
Abraham MH, Ibrahim A. Air to fat and blood to fat distribution of volatile organic compounds and drugs: Linear free energy analyses. Eur J Med Chem. 2006;41:1430–8.
Abraham MH, Ibrahim A, Acree Jr WE. Air to muscle and blood/plasma to muscle distribution of volatile organic compounds and drugs: Linear free energy analyses. Chem Res Toxicol. 2006;19:801–8.
Sprunger LM, Gibbs J, Acree Jr WE, Abraham MH. Linear free energy relationship correlation of the distribution of solutes between water and cetyltrimethylammonium bromide (CTAB) micelles. QSAR Comb Sci. 2009;28:72–88.
Wilson A, Tian A, Dabadge N, Acree Jr WE, Varfolomeev MA, Rakipov IT, et al. Enthalpy of solvation correlations for organic solutes and gases in dichloromethane and 1,4-dioxane. Struct Chem. 2013;24:1841–53.
Stephens TW, Chou V, Quay AN, Shen C, Dabadge N, Tian A, et al. Thermochemical investigations of solute transfer into ionic liquid solvents: Updated Abraham model equation coefficients for solute activity coefficient and partition coefficient predictions. Phys Chem Liquids. 2014;52:488–518.
Abraham MH, Smith RE, Luchtefeld R, Boorem AJ, Luo R, Acree Jr WE. Prediction of solubility of drugs and other compounds in organic solvents. J Pharm Sci. 2010;99:1500–15.
Bradley J-C, Acree Jr. WE, Lang ASID. An open logP model based upon Abraham Descriptors. ONS Challenge. Open Notebook. [http://onschallenge.wikispaces.com/ADlogP001]
Van Noort PCM. Solvation thermodynamics and the physical–chemical meaning of the constant in Abraham solvation equations. Chemosphere. 2012;87:125–31.
Grubbs LM, Saifullaha M, De La Rosa NE, Yea S, Achia SS, Acree Jr WE, et al. Mathematical correlations for describing solute transfer into functionalized alkane solvents containing hydroxyl, ether, ester or ketone solvents. Fluid Phase Equilib. 2010;298(1:15):48–53.
Sprunger LM, Achi SS, Acree Jr WE, Abraham MH. Reply to comments of Endo and Goss concerning “development of correlations for describing solute transfer into acyclic alcohol solvents based on the Abraham model and fragment-specific equation coefficients. Fluid Phase Equilibr. 2010;295:148–50.
The Chemistry Development Kit [http://sourceforge.net/projects/cdk]
Bradley J-C, Acree Jr. WE, Lang ASID. Compounds with known Abraham descriptors. figshare 2014. [http://dx.doi.org/10.6084/m9.figshare.1176994]
Bradley J-C, Lang ASID. Solvent coefficients for alternative ‘safe’ solvents. ONS Challenge. Open Notebook. [http://onschallenge.wikispaces.com/AbrahamSolventsModel003]
Guha R. CDK Descriptor UI. [https://github.com/rajarshi/cdkdescui]
Lang ASID. ONS Models in R. ONS Challenge. Open Notebook. [http://onschallenge.wikispaces.com/ONSModels-R]
Bradley J-C, Acree Jr. WE, Lang ASID. Solvent coefficients for sustainable solvents. ONS Challenge. Open Notebook. [http://onschallenge.wikispaces.com/AbrahamSolventsModel004]
Moity L, Durand M, Benazzouz A, Pierlot C, Molinier V, Aubry J-M. Panorama of sustainable solvents using the COSMO-RS approach. Green Chem. 2012;14:4.
Bradley J-C, Guha R, Hooker B, Koch SJ, Lang ASID, Neylon C, et al.: Open Notebook Science Solubility Challenge. Open Notebook [http://onschallenge.wikispaces.com/]
Bradley J-C, Lang ASID. Real-time prediction of Abraham descriptors by Chemspider ID for measured solubilities in the Open Notebook Science Challenge solubility database. Webservice. [http://showme.physics.drexel.edu/onsc/models/solventselector.php?csids=553171&limreact=0&limprod=0&bp=0&washes=0]
Grodowska K, Parczewski A. Organic solvents in the pharmaceutical industry. Acta Pol Pharm. 2010;67:3.
The authors declare that they have no competing interests.
J-CB conceived of the study and led the initiative to collect and model solubility values using ONS, MHA provided the physical-chemical interpretation of the models referenced in this paper and aided in manuscript preparation, WEA provided the Open Data list of compounds with known Abraham descriptors and aided in manuscript preparation, ASIDL performed all the modeling and aided in manuscript preparation. All authors read and approved the final manuscript.
Current list of dry solvents with known Abraham descriptors together with their c = 0 predicted values, SMILES, melting points, boiling points, and ChemSpider ID (CSID).
Predicted solvent coefficients for all 293 solvents considered in this study: sustainable, classic, and measured.