Subscriber access provided by AUT Library
Computational Chemistry
Calculating the Partition Coefficients of Organic Solvents in Octanol/Water and Octanol/Air Miroslava Nedyalkova, Sergio Madurga, Marek Tobiszewski, and Vasil Simeonov J. Chem. Inf. Model., Just Accepted Manuscript • DOI: 10.1021/acs.jcim.9b00212 • Publication Date (Web): 01 May 2019 Downloaded from http://pubs.acs.org on May 6, 2019
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1
Journal of Chemical Information and Modeling
Calculating the Partition Coefficients of Organic Solvents in Octanol/Water and Octanol/Air
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
Miroslava A. Nedyalkova*†, Sergio Madurga‡, Marek Tobiszewski §, Vasil Simeonov ⊥
†
Inorganic Chemistry Department, Faculty of Chemistry and Pharmacy, University of Sofia, Sofia 1164, Bulgaria ‡ Departament de Ciència de Materials i Química Física & Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08028 Barcelona, Catalonia, Spain § Department of Analytical Chemistry, Faculty of Chemistry, Gdańsk University of Technology (GUT), 80-233 Gdańsk, Poland ⊥ Analytical Chemistry Department, Faculty of Chemistry and Pharmacy, University of Sofia, Sofia 1164, Bulgaria ABSTRACT Partition coefficients define how a solute is distributed between two immiscible phases at
19
equilibrium. The experimental estimation of partition coefficients in a complex system can be an
20
expensive, difficult and time-consuming process. Here, a computational strategy is presented to
21
predict the distribution of a set of solutes in two relevant phase equilibria. The octanol/water and
22
octanol/air partition coefficients are predicted for a group of polar solvents using density functional
23
calculations in combination with a solvation model based on density (SMD) and are in excellent
24
agreement with experimental data. Thus, the use of quantum chemical calculations to predict
25
partition coefficients from free energies should be a valuable alternative for unknown solvents.
26
The obtained results indicate that the SMD continuum model in conjunction with any of the three
27
DFT functionals (B3LYP, M06-2X and M11) agrees with the observed experimental values. The
28
highest correlation to experimental data for the octanol/water partition coefficients was reached by
29
the M11 functional; for the octanol/air partition coefficient, the M06-2X functional yielded the
30
best performance. To the best of our knowledge, this is the first computational approach for the
ACS Paragon Plus Environment
1
Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 2 of 18
31
prediction of octanol/air partition coefficients by DFT calculations, which has remarkable
32
accuracy and precision.
33
INTRODUCTION
34
Physical properties of molecules can be used, for example, to make predictions about the
35
environmental fate of unknown solvents. The lack of physical property data may be resolved by
36
the development of different computational methodologies for the prediction of the appropriate
37
physicochemical properties. In particular, a variety of methods have been developed to predict
38
partition coefficients from a chemical structure.1-9 One of the most important prediction methods
39
is based on the quantitative structure-property relationship (QSPR).10-12 Various algorithms and
40
online platforms based on QSPR, such as AlogPs, ADMET predictor, and ACD/logD, have been
41
developed. The main approach of these methods is based on finding the appropriate set of
42
molecular descriptors that allow the precise reproduction of a given physical property using a large
43
database of available experimental data. Accurate QSPR models are obtained when this method is
44
applied to molecules that resemble the ones in the database used to build the model. Thus, the
45
weakness of QSPR models is related to the prediction of properties of molecules that vary slightly
46
from those used in the database.13-15
47
Alternatively, partition coefficients can be predicted by taking into account that this property is
48
related to the free energy difference of a solute in different solvents. In the work of Bannan et al.16,
49
the computational scheme that was used consisted of molecular dynamics simulations with explicit
50
solvent molecules to obtain transfer free energies between the solvents, which in turn were used
51
to calculate their partition coefficients. The partition coefficients between octanol/water and
52
cyclohexane/water were obtained using the generalized AMBER force field (GAFF) and the
53
dielectric corrected GAFF (GAFF-DC).
ACS Paragon Plus Environment
2
Page 3 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Journal of Chemical Information and Modeling
54
Jones et al. used ab initio calculations to predict the cyclohexane/water partition coefficient for a
55
set of 53 compounds. The free energy of transfer was calculated with several DFT functionals in
56
combination with the solvation model based on density implicit solvation model, obtaining a good
57
estimation.17 Rayne et al. computed air-water partition coefficients for a dataset of 86 large
58
polycyclic aromatic hydrocarbons and their unsaturated relatives by means of high-level theory
59
G4(MP2) gas and aqueous phase calculations with the SMD, integral equation formalism variant
60
(IEFPCM-UFF), and conductor-like polarizable continuum model (CPCM) solvation models.18
61
The results of the three solvation models for a range of neutral and ionic compounds show accurate
62
air/water partition coefficients (Kaw). Better accuracy is obtained when using higher levels of
63
theory.
64
In the work of Michalik19, the octanol/water partition coefficients for 27 alkane alcohols were
65
predicted by quantum chemical calculations with three solvation models. The results were in rather
66
good agreement with the corresponding experimental values. When comparing their results with
67
other implicit solvent models (IEF-PCM or C-PCM), the authors observed deviation of the
68
linearity. This mixed quantum (DFT) – QSPR analysis has been recently implemented successfully
69
in the prediction of pKa values for carboxylic acids.20
70
In the present work, the prediction of octanol/water and air/water partition coefficients of a set of
71
55 organic solvents was carried out by means of DFT calculations. Solvation free energies were
72
computed at the DFT levels of theory to estimate the partition coefficients. Good correlation
73
coefficients were obtained between the calculated and experimentally measured values.
74
ACS Paragon Plus Environment
3
Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 4 of 18
75
DATASET:
76
The dataset consisted of 55 molecules that were selected by a previous study.21 In that study, 150
77
solvents were clustered on the basis of physicochemical properties—melting and boiling points,
78
density, water solubility, vapor pressure, Henry’s law constant, logP, logKoa and surface tension.
79
Molecules included in the present study belong to polar or nonpolar groups. The experimental
80
values for the octanol/water and air/water partition coefficients are presented in Supporting
81
information - S1. The majority of these 55 molecules have been considered to be green solvents
82
in various literature reports.21 Some of the included molecules, such as ether-alcohols, are poorly
83
characterized in terms of their hazards or physicochemical properties.22
84 85
COMPUTATIONAL METHODS
86
All calculations presented in this work were performed using the Gaussian 09 quantum chemistry
87
package.23 Molecular structures were generated in the more extended conformation using
88
GaussView 5.0. 24
89
Geometries of all 55 molecules were optimized with three levels of theory using M06-2X 24, M1126
90
and B3LYP27, with the 6-311+G** basis set using the continuum solvation model based on density,
91
SMD28. Hessian analysis indicated no existence of imaginary frequencies, proving that all
92
optimized structures were true minima. SMD can be used as a universal solvation model because
93
it can be applied to any charged or uncharged solute in any type of solvent. The parameters required
94
for the solvent are the dielectric constant, refractive index, bulk surface tension, and acidity and
ACS Paragon Plus Environment
4
Page 5 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Journal of Chemical Information and Modeling
95
basicity parameters. This model divides the solvation free energy into two main contributions—
96
the bulk electrostatic contribution and the cavity dispersion contribution.
97
To calculate the octanol/water partition coefficient, the obtained SMD free energies in both
98
solvents at 298.15 K were used to calculate the standard free energy associated with the transfer
99
of the solute from the aqueous phase (w) to octanol (o):
100 101
Δ𝐺0𝑜 𝑤 = Δ𝐺0𝑜 ― Δ𝐺0𝑤 (1) The octanol/water partition coefficient was then calculated according to ΔG0o w
102
log𝑃𝑜 𝑤 = ― 2.303𝑅𝑇 (2)
103
To calculate the octanol/air partition coefficient, Po/a, the SMD solvation free energy in octanol
104
was obtained from the free energies of the molecule in gas phase and in octanol:
105 106
Δ𝐺0𝑠𝑜𝑙𝑣,𝑜 = Δ𝐺0𝑜 ― Δ𝐺0𝑔𝑎𝑠 (3) The octanol/air partition coefficient was then calculated according to 𝛥𝐺0𝑆𝑜𝑙𝑣,𝑜
107
log𝑃𝑜 𝑎 = ― 2.303𝑅𝑇 (4)
108
RESULTS
109
The structures of the fifty-five (55) solvents under study are shown in Figure 1. All the
110
correlation coefficients, slopes, and intercepts for all molecules are collected in Table 1 and
111
Table 2. The linear correlation between the experimental and calculated values for octanol/water
112
and octanol/air partition coefficients for the observed organic solvent dataset are presented
113
graphically in Figure 2 and Figure 3.
114 115
ACS Paragon Plus Environment
5
Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 6 of 18
116 ACS Paragon Plus Environment
6
Page 7 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Journal of Chemical Information and Modeling
117 ACS Paragon Plus Environment
7
Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 8 of 18
118 119
Figure 1. The chemical structures of solvents.
120 121 122 123
ACS Paragon Plus Environment
8
Page 9 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Journal of Chemical Information and Modeling
124 125
Notably, a high linear correlation was obtained with the functional M11 for the LogP values. The
126
statistical variety of error metrics were calculated to compare the calculated LogP and logKoa,
127
values (illustrated in Tables 1 and 2) to the experimentally determined measures using root-mean-
128
squared error (RMSE), Pearson’s correlation coefficient (R), mean absolute deviation (MAD),
129
mean square error (MSE) and mean absolute percentage error (MAPE), along with standard errors.
130
The assumed outliers observed in the first data assessment (Supporting Information – 2) were
131
eliminated in the next level of linear regression by applying the 4-sigma rule30 for the detection of
132
outliers. Extremely large and extremely small values for both coefficients were observed in our
133
data. The 4-sigma region represents 99.99% of the values for a normal distribution and 97% for
134
symmetric unimodal distributions. In our case, the outliers were adequately identified and
135
excluded from further statistical treatment of the dataset. New plots were established and are
136
presented below.
137
These metrics allow us to make meaningful comparisons between experimental and calculated
138
data with the SMD solvent model and the three functionals. As previously mentioned, the best
139
results were obtained with the M11 functional with a RMSE of 0.72 and a Pearson’s correlation
140
coefficient (R) of 0.99.
141
The obtained statistics show a good accuracy of the computed partition coefficients when applying
142
the M11 functional. In general, the sign of the calculated and experimental data were in agreement.
143
Positive values indicate a preference for the organic phase, and negative values indicate a
144
preference for the aqueous phase. For DFT-based SMD solvation models, this approach appears
145
to be appropriate for obtaining good correlations with the experimental measurements in the
146
calculation of logP. Presumably, the accuracy of the results is due to the appropriated
147
parametrization of the SMD solvation model to yield accurate solvation free energies.
ACS Paragon Plus Environment
9
Journal of Chemical Information and Modeling
148
Level of theory – logP DFT-M062X DFT-M11 DFTB3LYP
149 150 151 152 153 154 155 156
Slope
Intercept
R2
R
MAD
MSE
RMSE
MAPE
Standard Error
0.9498
0.020
0.9726
0.9879
0.29
0.41
0.48
0.56
0.65
0.9816
0.5398
0.9841
0.9920
0.75
0.96
0.72
0.57
0.72
1.022
0.8597
0.9609
0.9802
0.86
0.96
0.90
0.58
0.86
Table 1. Linear regression parameters obtained for the calculated logP with respect to the experimental data, and statistical error assessment of the linear regression terms based on applied computational models.
Level of theory – logKoa
Slope
Intercept
R2
R
MAD
MSE
RMSE
MAPE
DFT-M062X
1.0111
0.1806
0.9780
0.9889
0.40
0.39
0.59
0.95
Standard Error
0.3 DFT-M11 DFTB3LYP
157 158 159 160 161 162
0.9816
0.5398
0.8891
0.9428
2.18
1.45
1.20
0.90
0.8410
0.9527
0.8618
0.9339
1.78
1.92
1.33
0.63
0.6 0.75
Table 2. Linear regression parameters obtained for the calculated logKoa with respect to the experimental data, and statistical error assessment of the linear regression terms based on applied computational models.
8 Calculated B3LYP
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 10 of 18
6 4 2 0
-2
0 -2
163 164 165
2
4
6
8
Experimentl logP
ACS Paragon Plus Environment
1 0
Page 11 of 18
6
M11
Calculated M11
8
4 2 0
-2
0
166 167
2
4
6
8
Experimental logP
-2 8
4
M11
Calculated M06-2X
6
2 0
-2
168 169 170 171 172 173
0
2
4
6
8
Experimental logP
-2
Figure 2. Comparison of the calculated logP to the experimental values. The upper graph shows the B3LYP model results, the middle graph shows M11, and the bottom graph shows M06-2X.
8 Calculated B3LYP
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Journal of Chemical Information and Modeling
6 4 2 0
-2
174 175 176
0 -2
2
4
6
8
Experimental logKoa
ACS Paragon Plus Environment
1 1
Journal of Chemical Information and Modeling
8
4
M11
Calculated M11
6
2 0
-2
0
2
4
6
8
Experimental logKoa
-2
177 178 179 180
8 6 4 2
M11
Calculated M06-2X
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 12 of 18
0 -2
0 -2
181 182 183 184 185 186 187 188
2
4
6
8
Experimental logKoa
Figure 3. Comparison of the calculated logKoa to the experimental values. The upper graph shows the B3LYP model results, the middle graph shows M11, and bottom graph shows M06-2X.
The predicted values for logKoa allowed us to evaluate the solvation free energy obtained with
189
the SMD model. The various error statistical metrics were calculated and are depicted in Table 2.
190
The best agreement between the calculated and experimental partition coefficients was obtained
191
with the M06-2X functional. This indicates a good estimation of the solvation free energy using
192
the M06-2X functional. These findings provide useful information for understanding the
ACS Paragon Plus Environment
1 2
Page 13 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Journal of Chemical Information and Modeling
193
partitioning, and the proposed computational scheme can be applied to other unknown solvents.
194
For the other functionals, the obtained results also indicate that the approach is valid.
195
In all calculations, a single conformation (the more extended conformation) was used for each
196
molecule in each phase. Although this conformation was verified to be a minimum in the gas phase
197
and in both solvents, it is possible that multiple conformations in the gas and/or solvents should be
198
taken into account to obtain accurate descriptions.
199 200
CONCLUSIONS
201
The accuracy of predicted partition coefficients, octanol/water and octanol/air were estimated
202
within the framework of DFT, with three functionals represented by B3LYP, and two popular (but
203
different) meta-GGA Minnesota functionals, M06-2X and M11. The best correlation with the
204
experiment data was achieved with the M11 and M06-2X functionals. The quality of the regression
205
model was significantly improved after excluding points that were defined as outliers. Thus, this
206
computational protocol could be used as a tool for newly synthesized solvents where there is a lack
207
of data regarding the octanol/water and octanol/air partition coefficients. The methodology is
208
useful for predicting missing data for envirometrics data interpretation and modeling.
209 210
Highlights
211 212 213 214 215 216 217
-
Prediction of the octanol/water partition coefficient for 55 green solvents Prediction of the octanol/air partition coefficient for 55 green solvents Accurate application of DFT for logP and logKoa Cross-validation of standard error for estimated models
AUTHOR INFORMATION
ACS Paragon Plus Environment
1 3
Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
218
Corresponding Author
219
*E-mail:
[email protected]. Phone: 00359 885 76 33 98
220
Author Contributions
221
All authors contributed to the writing of this manuscript. All authors have approved the
222
final version of the manuscript.
223
Notes
224
The authors declare no competing financial interests.
Page 14 of 18
225 226
ACKNOWLEDGMENTS
227
The authors are thankful for the support through the project BG05M2OP001-1.001-0004-C01/28.02.2018 (2018-
228
2023) and for the computational resources provided by the project HPC-EUROPA3 (INFRAIA-2016-1-730897), with
229
the support of the EC Research Innovation Action under the H2020 Programme. The author Miroslava Nedyalkova
230
gratefully acknowledges the computer resources and technical support provided by the HPC center - Barcelona
231
Supercomputing Center-Centro Nacional de Supercomputación (BSC-CNS). The financial support from Generalitat
232
de Catalunya (Grants 2017SGR1033) and the Spanish Structures of Excellence María de Maeztu program through
233
grant MDM-2017–0767 is fully acknowledged.
234 235 236 237 238 239 240 241 242 243 244
SUPPORTING INFORMATION Supporting Information Available: Additional Table with the experimental values for the octanol/water and air/water partition coefficients are presented in Supporting information – S1 Additional information about the applied Evaluation Metrics: MAD, MSE, RMSE, MAPE in Supporting information – S2 Additional plots with the deviated from the linearity points in Supporting information – S2
REFERENCES
ACS Paragon Plus Environment
1 4
Page 15 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
245
Journal of Chemical Information and Modeling
1.
Kohler, M. G.; Grigoras, S.; Dunn III, W. J. The Relationship between Chemical
246
Structure and the Logarithm of the Partition Coefficient. Quant. Struct. -Act. Relat. 1988,
247
150-159.
248
2.
Kasai, K.; Umeyama, H.; Tomonaga, A. The Study of Partition Coefficients. The
249
Prediction of Log P Value Based on Molecular Structure. Bull. Chem. Soc. Jpn. 1988, 61,
250
2701-2706.
251
3.
252 253
Bodor, N.; Gabayani, Z.; Wong, C. A New Method for the Estimation of Partition Coefficient. J. Am. Chem. Soc. 1989, 111, 3783-3786.
4.
Chou, J. T.; Jurs, P. C. Computer-Assisted Computation of Partition Coefficients from
254
Molecular Structure Using Fragment Constants. J. Chem. Inf. Comput. Sci. 1979, 19,
255
172-178.
256
5.
Sasaki, Y.; Kubodera, H.; Matuszaki, T.; Umeyama, H. Prediction of Octanol/Water
257
Partition Coefficients Using Parameters Derived from Molecular Structures. J.
258
Pharmacobio-Dyn. 1991, 14, 207-214.
259
6.
Klopman, G.; Li, J.-Y.; Wang, S.; Dimayuga, M. Computer Automated Log P
260
Calculation Based on an Extended Group Contribution Approach. J. Chem. Inf. Comput.
261
Sci. 1994, 34, 752-781.
262
7.
263 264
Maylan, W. M.; Howard, P. H. Atom/Fragment Contribution Method for Estimating Octanol-Water Partition Coefficients. J. Pharm. Sci. 1995, 84, 83-92.
8.
Ghose, A.; Crippen, G. M. Atomic Physicochemical Parameters for Three-Dimensional
265
Structure-Directed Quantitative Structure-Activity Relationships I. Partition Coefficients
266
as a Measure of Hydrophobicity. J. Comput. Chem. 1986, 4, 565-577.
267
9.
268 269
Biology; Wiley Interscience: New York, 1979; p 19 10.
270 271
Beck, B.; Breind,l.; Clark, T. QM/NN QSPR models with error estimation: vapor pressure and logP. Chem. Inf. Comput. Sci., 2000, 40, A. 1046–1051
11.
272 273
Hansch, C.; Leo, A. Substituent Constants for Correlation Analysis in Chemistry and
Šoškić, M.; Plavšić. D. Modeling the Octanol−Water Partition Coefficients by an Optimized Molecular Connectivity Index. J. Chem. Inf. Model. 2005, 4, 930-938
12.
Duchowicz, P.; Castro, E.; Toropov A.; Nesterov I.; and O. M. Nabiev. QSPR modeling
274
the aqueous solubility of alcohols by optimization of correlation weights of local graph
275
invariants. Mol. Diversity, 2004, 8, 325–330.
ACS Paragon Plus Environment
1 5
Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
276
13.
Page 16 of 18
Dyekjaer, J.; Rasmussen K.; Jonsdottir, S. QSPR models based on molecular mechanics
277
and quantum chemical calculations. 1. Construction of Boltzmann-averaged descriptors for
278
alkanes, alcohols, diols, ethers and cyclic compounds. J. Mol. Model.,2002, 8, 277–289.
279
14.
Doweyko, A. 3D-QSAR illusions. J. Comput.-Aided Mol. Des., 2004, 18, 587–596.
280
15.
L. D. Hughes, D. S. Palmer, F. Nigsch and J. B. Mitchell, J. Chem. Inf. Model., 2008, 48,
281 282
220–232. 16.
Bannan, C.; Calabró,G.; Kyu,D and Mobley, D Calculating Partition Coefficients of
283
Small Molecules in Octanol/Water and Cyclohexane/Water J. Chem. Theory Comput.
284
2016, 12, 4015−4024
285
17.
Jones, M.R.; Brooks, B.R.; Wilson, A.K., Partition Coefficients for the SAMPL5
286
Challenge Using Transfer Free Energies. J. Comput.-Aided Mol. Des. 2016, 30, 1079–
287
1086
288
18.
Rayne, S.; Forest, K. Air-water partition coefficients for a suite of polycyclic aromatic
289
and other C10 through C20 unsaturated hydrocarbons, J. Environ. Sci and health, A,
290
2016, 51, 938-953
291
19.
292 293
Michalík,M.; Lukeš,V.; The validation of quantum chemical lipophilicity prediction of alcohols Acta Chimica Slovaca, 2016, 2, 89-94
20.
Guillermo Caballero-García, Gustavo Mondragón-Solórzano, Raúl Torres-
294
Cadena, Marco Díaz-García, Jacinto Sandoval-Lira, Joaquín Barroso-Flores
295
“Calculation of VS,max and Its Use as a Descriptor for the Theoretical
296
Calculation of pKa Values for Carboxylic Acids” Molecules 2019, 24, 79-93
297
21.
Tobiszewski, M., Tsakovski, S., Simeonov, V., Namieśnik, J., Pena-Pereira, F A solvent
298
selection guide based on chemometrics and multicriteria decision analysis Green
299
Chemistry, 2015, 17, 4773-4785
300
22.
Pena-Pereira, F., Kloskowski, A., Namieśnik, J Perspectives on the replacement of
301
harmful organic solvents in analytical methodologies: a framework toward the
302
implementation of a generation of eco-friendly alternatives Green Chemistry, 2015, 17,
303
3687-3705
ACS Paragon Plus Environment
1 6
Page 17 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
304
Journal of Chemical Information and Modeling
23.
Gaussian 09, Revision A.02, Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G.
305
E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Petersson, G. A.;
306
Nakatsuji, H.; Li, X.; Caricato, M.; Marenich, A.; Bloino, J.; Janesko, B. G.; Gomperts,
307
R. B.; Mennucci, H.; Hratchian, P.; Ortiz, J. V.; Izmaylov, A. F.; Sonnenberg, J. L.;
308
Williams-Young, D.; Ding, F.; Lipparini, F.; Egidi, F.; Goings, J.;Peng, B.; Petrone, A.;
309
Henderson, T.; Ranasinghe, D.; Zakrzewski, V. G.; Gao, J.; Rega, N.; Zheng, G.; Liang,
310
W.; Hada, M.; Ehara, M.;Toyota, K.; Fukuda, R. Hasegawa, J.; Ishida, M.; Nakajima, T.;
311
Honda, Y.;Kitao, O.;Nakai, H.;Vreven, T.; Throssell, K.; Montgomery, J. A.; Peralta, Jr.,
312
J. E.; Ogliaro. F. ; Bearpark. M. ; Heyd, J.; Brothers, E.; Kudin. K. N.; Staroverov; Keith,
313
V. N.; Kobayashi,T.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar,
314
J. C.; Tomasi, J.; Cossi, M.; M. J.; Millam, Klene, M.; Adamo, C.; Cammi; Ochterski, J.;
315
Martin, R. L.; Morokuma, K.; Farkas, O.; Foresman, J. B.; Fox, D. J. Gaussian, Inc.,
316
Wallingford CT, 2016
317
24.
318 319
Dennington, R.; Keith, T.; Millam, J. Gaussview, Version 5; Semichem, Inc.: Shawnee Mission, KS, 2009
25.
Zhao, Y.;Truhlar, D. The M06 suite of density functionals for main group
320
thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and
321
transition elements: two new functionals and systematic testing of four M06-class
322
functionals and 12 other functionals Theor Chem Account 2008, 120, 215-241
323
26.
Peverati, R.; Truhlar, DG Performance of the M11 and M11-L density functionals for
324
calculations of electronic excitation energies by adiabatic time-dependent density
325
functional theory Phys Chem Chem Phys. 2012, 14, 1363-70
326
27.
327 328
Becke, A. D.; Density‐functional thermochemistry. III. The role of exact exchange J. Chem. Phys. 1993, 98, 5648-565
28.
Marenich, A. V., Cramer, C. J., Truhlar D. G Universal solvation model based on solute
329
electron density and a continuum model of the solvent defined by the bulk dielectric
330
constant and atomic surface tensions J. Phys. Chem. B. 2009, 113, 6378-96
331
29.
Hodges, G., Eadsfort,C., Bossuyt, B., Bouvy, A., Enrici, M-H., Geurts, M., Kotthoff, M.,
332
Michie, E., Miller, D., Müller,J., Oetter, G., JRoberts,J., Schowanek,D., Joachim
333
Venzmer A comparison of log Kow (n-octanol–water partition coefficient) values for
ACS Paragon Plus Environment
1 7
Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 18 of 18
334
non-ionic, anionic, cationic and amphoteric surfactants determined using predictions and
335
experimental methods Environ Sci Eur. 2019, 31, 1-18
336
30.
Sachs, L. (1984). Applied Statistics. New York, NY: Springer New York.
337 338
Graphical Representation:
339
340
ACS Paragon Plus Environment
1 8