COSMO-RS blind prediction of distribution coefficients and aqueous pKa values from the SAMPL8 challenge.

The SAMPL8 blind prediction challenge, which addresses the acid/base dissociation constants (pKa) and the distribution coefficients (logD), was addressed by the Conductor like Screening Model for Realistic Solvation (COSMO-RS). Using the COSMOtherm implementation of COSMO-RS together with a rigorous conformational sampling, yielded logD predictions with a root mean square deviation (RMSD) of 1.36 log units over all 11 compounds and seven bi-phasic systems of the data set, which was the most accurate of all contest submissions (logD).For the SAMPL8 pKa competition, participants were asked to report the standard state free energies of all microstates, which were then used to calculate the macroscopic pKa. We have used COSMO-RS based linear free energy fit models to calculate the requested energies. The assignment of the calculated and experimental pKa values was made on the basis of the popular transitions, i.e. the transition hat was predicted by the majority of the submissions. With this assignment and a model that covers both, pKa and base pKa, we achieved an RMSD of 3.44 log units (18 pKa values of 14 molecules), which is the second place of the six ranked submissions. By changing to an assignment that is based on the experimental transition curves, the RMSD reduces to 1.65. In addition to the ranked contribution, we submitted two more data sets, one for the standard pKa model and one or the standard base pKa model of COSMOtherm. Using the experiment based assignment with the predictions of the two sets we received a RMSD of 1.42 log units (25 pKa values of 20 molecules). The deviation mainly arises from a single outlier compound, the omission of which leads to an RMSD of 0.89 log units.

Quantifying Uncertainties in Solvation Procedures for Modeling Aqueous Phase Reaction Mechanisms.

Computational quantum chemistry provides fundamental chemical and physical insights into solvated reaction mechanisms across many areas of chemistry, especially in homogeneous and heterogeneous renewable energy catalysis. Such reactions may depend on explicit interactions with ions and solvent molecules that are nontrivial to characterize. Rigorously modeling explicit solvent effects with molecular dynamics usually brings steep computational costs while the performance of continuum solvent models such as polarizable continuum model (PCM), charge-asymmetric nonlocally determined local-electric (CANDLE), conductor-like screening model for real solvents (COSMO-RS), and effective screening medium method with the reference interaction site model (ESM-RISM) are less well understood for reaction mechanisms. Here, we revisit a fundamental aqueous hydride transfer reaction-carbon dioxide (CO2) reduction by sodium borohydride (NaBH4)-as a test case to evaluate how different solvent models perform in aqueous phase charge migrations that would be relevant to renewable energy catalysis mechanisms. For this system, quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations almost exactly reproduced energy profiles from QM simulations, and the Na+ counterion in the QM/MM simulations plays an insignificant role over ensemble averaged trajectories that describe the reaction pathway. However, solvent models used on static calculations gave much more variability in data depending on whether the system was modeled using explicit solvent shells and/or the counterion. We pinpoint this variability due to unphysical descriptions of charge-separated states in the gas phase (i.e., self-interaction errors), and we show that using more accurate hybrid functionals and/or explicit solvent shells lessens these errors. This work closes with recommended procedures for treating solvation in future computational efforts in studying renewable energy catalysis mechanisms.

Prediction of cyclohexane-water distribution coefficients with COSMO-RS on the SAMPL5 data set.

The Conductor-Like-Screening-Model for Real Solvents (COSMO-RS) method has been used for the blind prediction of cyclohexane-water distribution coefficients logD within the SAMPL challenge. The partition coefficient logP of the neutral species was calculated first and then corrected for dissociation or protonation, as appropriate for acidic or basic solutes, to obtain the cyclohexane-water logD. Using the latest version of the COSMOtherm implementation, this approach in combination with a rigorous conformational sampling yielded a predictive accuracy of 2.11 log units (RMSD) for the 53 compounds of the blind prediction dataset. By that it was the most accurate of all contest submissions and it also achieved the best rank order. The RMSD mainly arises from a group of outliers in the negative logD range, which at least partly may arise from dimerization or other experimental problems coming up for very polar molecules in very non-polar solvents.

Interpretation of experimental hydrogen-bond enthalpies and entropies from COSMO polarisation charge densities.

In this work, experimental hydrogen-bond (HB) enthalpies measured in previous works for a wide range of acceptor molecules in dilute mixtures of 4-fluorophenol in non-polar solvents are quantified from COSMO polarisation charge densities σ of HB acceptors (HBA). As well as previously demonstrated for quantum chemically calculated HB enthalpies, a good correlation of the experimental data with the polarisation charge densities is observed, covering an extended range of HBA (O, N, S, π systems and halogens) ranging from very weak to strong hydrogen bonds. Furthermore, for the first time, a quantitative analysis of experimental HB entropies is performed for such a chemical diversity of HBA. A good quantification of these entropies is achieved using the polarisation charge density σ as a descriptor in combination with the logarithm of a directional partition function Ω(HB). This partition function covers the directional and multiplicity entropy of HBA and is based on the σ-proportional HB enthalpy expression taken from COSMO-RS. As a result, the experimental HB enthalpies and free energies of the ~300 HB complexes are quantified with an accuracy of ~2 kJ mol(-1) based on COSMO polarisation charge densities.

Polarization charge densities provide a predictive quantification of hydrogen bond energies.

A systematic density functional theory based study of hydrogen bond energies of 2465 single hydrogen bonds has been performed. In order to be closer to liquid phase conditions, different from the usual reference state of individual donor and acceptor molecules in vacuum, the reference state of donors and acceptors embedded in a perfect conductor as simulated by the COSMO solvation model has been used for the calculation of the hydrogen bond energies. The relationship between vacuum and conductor reference hydrogen bond energies is investigated and interpreted in the light of different physical contributions, such as electrostatic energy and dispersion. A very good correlation of the DFT/COSMO hydrogen bond energies with conductor polarization charge densities of separated donor and acceptor atoms was found. This provides a method to predict hydrogen bond strength in solution with a root mean square error of 0.36 kcal mol(-1) relative to the quantum chemical dimer calculations. The observed correlation is broadly applicable and allows for a predictive quantification of hydrogen bonding, which can be of great value in many areas of computational, medicinal and physical chemistry.

Prediction of acidity in acetonitrile solution with COSMO-RS.

The COSMO-RS method, a combination of the quantum chemical dielectric continuum solvation model COSMO with a statistical thermodynamics treatment for realistic solvation simulations, has been used for the prediction of pK(a) values in acetonitrile. For a variety of 93 organic acids, the directly calculated values of the free energies of dissociation in acetonitrile showed a very good correlation with the pK(a) values (r(2) = 0.97) in acetonitrile, corresponding to a standard deviation of 1.38 pK(a) units. Thus, we have a prediction method for acetonitrile pK(a) with the intercept and the slope as the only adjusted parameters. Furthermore, the pK(a) values of CH acids yielding large anions with delocalized charge can be predicted with a rmse of 1.12 pK(a) units using the theoretical values of slope and intercept resulting in truly ab initio pK(a) prediction. In contrast to our previous findings on aqueous acidity predictions the slope of the experimental pK(a) versus theoretical DeltaG(diss) was found to match the theoretical value 1/RT ln (10) very well. The predictivity of the presented method is general and is not restricted to certain compound classes. However, a systematic correction of -7.5 kcal mol(-1) is required for compounds that do not allow electron-delocalization in the dissociated anion. The prediction model was tested on a diverse test set of 129 complex multifunctional compounds from various sources, reaching a root mean square deviation of 2.10 pK(a) units.

Prediction of the free energy of hydration of a challenging set of pesticide-like compounds.

In a blind validation test the COSMO-RS method, a combination of the quantum chemical dielectric continuum solvation model COSMO with a statistical thermodynamics treatment for more realistic solvation (RS) simulations, has been used for the direct prediction of transfer free energies of 55 demanding pesticide-like compounds. Comparison with experimental data yields an rms deviation of approximately 2 kcal/mol, which is in the order of the estimated inaccuracy of the experimental data. A detailed comparison reveals experimental and calculation pitfalls on conformational flexible, multifunctional, polar compounds.

Application of the conductor-like screening model for real solvents for prediction of the aqueous solubility of chlorobenzenes depending on temperature and salinity.

The conductor-like screening model for real solvents (COSMO-RS) is applied to the prediction of the aqueous solubility of chlorobenzenes (CBZs) in a liquid-liquid and liquid-solid equilibrium. The solubilities of CBZs at temperatures ranging from 5 to 60 degrees C are reported, and the enthalpies of solution are derived from van't Hoff plots. The salting effect of 12 chlorobenzenes is determined for sodium, potassium, and calcium solutions at concentrations ranging from 0.01 to 3.5 mol/L. A new experimental aqueous solubility value is given for the isotopic labeled 1,2,4-trichlorobenzene determined at 25 degrees C by the generator column procedure. The COSMO-RS method presented in this study can be useful to estimate the magnitude of the salt effect and temperature influence on the behavior of any chlorobezene under environmental conditions.

Prediction of soil sorption coefficients with a conductor-like screening model for real solvents.

Using a general theory for partition coefficients based on a quantum chemically derived conductor-like screening model for real solvents sigma-moment descriptors, the logarithmic soil sorption coefficients log K(oc) of a database of 440 compounds has been successfully correlated, achieving a standard deviation (root-means-squared [RMS]) of 0.62 log-units on the training set and a predictive RMS of 0.72 log-units on a more demanding test set. The quality of this generally applicable predictive approach is almost the same as that of a regression of log K(oc) with experimental log K(ow) values, which are the best correlations currently available. The error of this new predictive method is only approximately 43% of the error of a recently published model using a different quantum chemically based approach.

COSMO-RS: an alternative to simulation for calculating thermodynamic properties of liquid mixtures.

The conductor-like screening model for realistic solvation (COSMO-RS) method has been established as a novel way to predict thermophysical data for liquid systems and has become a frequently used alternative to force field-based molecular simulation methods on one side and group contribution methods on the other. Through its unique combination of a quantum chemical treatment of solutes and solvents with an efficient statistical thermodynamics procedure for the molecular surface interactions, it enables the efficient calculation of many properties that other methods can barely predict. This review presents a short delineation of the theory, the application potential and limitations of COSMO-RS, and its most important application areas.

Probing carboxylate Gibbs transfer energies via liquid|liquid transfer at triple phase boundary electrodes: ion-transfer voltammetry versus COSMO-RS predictions.

Understanding liquid|liquid ion transfer processes is important in particular for naturally occurring species such as carboxylates. In this study electrochemically driven mono-, di-, and tri-carboxylate anion transfer at the 4-(3-phenylpropyl)pyridine|aqueous electrolyte interface is investigated experimentally for a triple phase boundary system at graphite electrodes. The tetraphenylporphyrinato-Mn(III/II) redox system (Mn(III/II)TPP) dissolved in the water-immiscible organic phase (4-(3-phenylpropyl)pyridine) is employed for the quantitative study of the structure-Gibbs transfer energy correlation and the effects of the solution pH on the carboxylate transfer process. For di- and tri-carboxylates the partially protonated anions are always transferred preferentially even at a pH higher than the corresponding pK(a). COSMO-RS computer simulations are shown to provide a quantitative rationalisation as well as a powerful tool for predicting Gibbs free energy of transfer data for more complex functionalised carboxylate anions. It is shown that the presence of water in the organic phase has a major effect on the calculated Gibbs free energies.

Accurate prediction of basicity in aqueous solution with COSMO-RS.

The COSMO-RS method, a combination of the quantum chemical dielectric continuum solvation model COSMO with a statistical thermodynamics treatment for realistic solvation simulations, has been used for the prediction of base pK(a) constants. For a variety of 43 organic bases the directly calculated values of the free energies of dissociation in water showed a very good correlation with experimental base pK(a) values (r2 = 0.98), corresponding to a standard deviation of 0.56 pK(a) units. Thus, we have an a priori prediction method for base pK(a) with the regression constant and the slope as only adjusted parameters. In accord with recent findings for pK(a) acidity predictions, the slope of pK(a) vs. DeltaG(diss) was significantly smaller than the theoretically expected value of 1/RTln(10). The predictivity of the presented method is general and not restricted to certain compound classes, but systematic corrections of 1 and 2 pKa units for secondary and tertiary aliphatic amines are required, respectively. The pK(a) prediction method was validated on a set of 58 complex multifunctional drug-like compounds, yielding an RMS accuracy of 0.66 pK(a) units.

First Principles Calculations of Aqueous pKa Values for Organic and Inorganic Acids Using COSMO-RS Reveal an Inconsistency in the Slope of the pKa Scale.

The COSMO-RS method, a combination of the quantum chemical dielectric continuum solvation model COSMO with a statistical thermodynamics treatment for more realistic solvation (RS) simulations, has been used for the direct prediction of pKa constants of a large variety of 64 organic and inorganic acids. A highly significant correlation of r(2) = 0.984 with a standard deviation of only 0.49 between the calculated values of the free energies of dissociation and the experimental pKa values was found, without any special adjustment of the method. Thus, we have a theoretical a priori prediction method for pKa, which has the regression constant and the slope as only adjusted parameters. Such a method can be of great value in many areas of physical chemistry, especially in pharmaceutical and agrochemical industry. To our surprise, the slope of pKa vs ΔGdiss is only 58% of the theoretically expected value of 1/RTln(10). A careful analysis with respect to different contributions as well as a comparison with the work of other authors excludes the possibility that the discrepancy is due to weaknesses of the calculation method. Hence, we must conclude that the experimental pKa scale depends differently on the free energy of dissociation than generally assumed.

A comparison between the two general sets of linear free energy descriptors of Abraham and Klamt.

Two sets of molecular descriptors, the five experimental Abraham, and the five COSMOments of Klamt's COSMO-RS, have been compared for a data set of 470 compounds. Both sets are considered as almost complete sets of LFER. The two sets of descriptors are shown to exhibit a large overlap as far as their chemical content. The chemical information however is distributed differently in each set with the Abraham set incorporating extra information in the excess molar refraction descriptor E. Regression equations have been constructed to predict the experimental Abraham descriptors from theoretically calculated COSMOments. The chemical interpretation of these equations is however difficult because of the lack of clustering which characterizes the distribution of chemical information through the two sets of descriptors. The predictability of the regression equations is tested successfully using a reasonably large set of data, and the method is compared to recent attempts to calculate the Abraham descriptors from various theoretical bases.

Prediction of aqueous solubility of drugs and pesticides with COSMO-RS.

The COSMO-RS method, originally developed for the prediction of liquid-liquid and liquid-vapor equilibrium constants based on quantum chemical calculations, has been extended to solid compounds by addition of a heuristic expression for the Gibbs free energy of fusion. By this addition, COSMO-RS is now capable of a priori prediction of aqueous solubilities of a wide range of typical neutral drug and pesticide compounds. Only three parameters in the heuristic expression have been fitted on a data set of 150 drug-like compounds. On these data an rms deviation of 0.66 log-units was achieved. Later, the model was tested on a set of 107 pesticides, which have been critically selected based on two experimental data sources and by a crosscheck with an independent HQSAR model. On this data set an rms of 0.61 log-units was achieved, without any adjustments to the structurally extremely diverse pesticides. This result verifies the ability of this extended COSMO-RS to predict aqueous solubilities of drugs and pesticides of almost arbitrary structural classes. The new method is COSMO-RSol.