Calculating the Partition Coefficients of Organic Solvents in Octanol

Publication Date (Web): May 1, 2019. Copyright © 2019 American Chemical Society. *E-mail: [email protected]. Phone: 00359 885 76 33 98. This article ...
1 downloads 0 Views 923KB Size
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