Solvent-water Partition Research Articles

COSMO-RS solute partition ratios for solvent mixtures of unknown composition: Henry's law constants as descriptors for mixture sigma profiles

Henry's law constants (H) for selected probe molecules have been used as descriptors to estimate the COSMO-RS sigma profiles of solvents and solvent mixtures. Henry's law constants were calculated with COSMOtherm for small sets of probe molecules in 155 organic solvents (training set), and these constants subsequently used as descriptors to model the solvent sigma profiles with 61 multiple linear regression (MLR) equations. Subsequent input into COSMOtherm of weighted basis molecule solvent mixtures whose sigma profiles closely matched those modelled for the training set solvents allowed estimation of air-solvent and water-solvent partition ratios for solutes in solvents and solvent mixtures without input of the solvent or solvent mixture identity. The best performing model had 16 descriptors and gave both a training and test set average root-mean square error (RMSE) of 0.008 and an average relative square error (RSE) of 0.07. Partition ratios (K) were then generated for a test set of 251 additional organic solute molecules in solvent/water media where solvents were test set compounds and H constants for the same probe molecules were used as descriptors. The best performing sigma profile model yielded log K RMSE values ranging from 0.17 to 0.92. Finally, this approach was applied to several mixtures ranging from simple binary mixtures to two mixtures considered to be of unknown or variable composition, complex reaction productions or biological materials (UVCBs), namely gasoline and an essential oil mixture. Mixture/water partition ratios were estimated for 251 solutes giving log K RMSE values ranging from 0.24 to 0.88.

Modeling the Partitioning of Anionic Carboxylic and Perfluoroalkyl Carboxylic and Sulfonic Acids to Octanol and Membrane Lipid.

Perfluoroalkyl carboxylic and sulfonic acids (PFCAs and PFSAs, respectively) have low acid dissociation constant values and are, therefore, deprotonated under most experimental and environmental conditions. Hence, the anionic species dominate their partitioning between water and organic phases, including octanol and phospholipid bilayers which are often used as model systems for environmental and biological matrices. However, data for solvent-water (SW) and membrane-water partition coefficients of the anion species are only available for a few per- and polyfluoroalkyl substances (PFAS). In the present study, an equation is derived using a Born-Haber cycle that relates the partition coefficients of the anions to those of the corresponding neutral species. It is shown via a thermodynamic analysis that for carboxylic acids (CAs), PFCAs, and PFSAs, the log of the solvent-water partition coefficient of the anion, log KSW (A- ), is linearly related to the log of the solvent-water partition coefficient of the neutral acid, log KSW (HA), with a unity slope and a solvent-dependent but solute-independent intercept within a PFAS (or CA) family. This finding provides a method for estimating the partition coefficients of PFCAs and PFSAs anions using the partition coefficients of the neutral species, which can be reliably predicted using quantum chemical methods. In addition, we have found that the neutral octanol-water partition coefficient, log KOW , is linearly correlated to the neutral membrane-water partition coefficient, log KMW ; therefore, log KOW , being a much easier property to estimate and/or measure, can be used to predict the neutral log KMW . Application of this approach to KOW and KMW for PFCAs and PFSAs demonstrates the utility of this methodology for evaluating reported experimental data and extending anion property data for chain lengths that are unavailable. Environ Toxicol Chem 2023;00:1-12. © 2023 SETAC.

Abraham Model Descriptors for Melatonin; Prediction of Solution, Biological and Thermodynamic Properties

Literature solubilities have been used to obtain properties or descriptors of melatonin. These indicate the chemical nature of melatonin: it is dipolar and has moderate hydrogen bond acidity and hydrogen bond basicity. The descriptors can be combined with equations that we have previously constructed to estimate water–solvent partition coefficients and solubilities in a huge number of organic solvents. In the same way, a range of biological properties can be estimated. These include blood–tissue partitions, water–skin partition and permeability through skin.

Determination of Water–Solvent Partition Coefficients for Fluorescein: Evaluation of Descriptors for the Lactone Form and Prediction of Properties

Determination of Water–Solvent Partition Coefficients for Fluorescein: Evaluation of Descriptors for the Lactone Form and Prediction of Properties

Descriptors for adamantane and some of its derivatives

Literature data on solubilities of adamantane in organic solvents have been used to obtain properties, or descriptors, of adamantane. There is much less data on substituted adamantanes but we have been able to obtain descriptors for some 40 substituted adamantanes. These descriptors can then be used to estimate a wide range of physicochemical, environmental and other properties of the adamantanes. For the first time, the water-solvent partition coefficient and the gas-solvent partition coefficient into a large range of solvents, can be estimated, the latter being equivalent to Henry's Law constants. A variety of other important properties can also be estimated. These include vapor pressures, enthalpies of vaporization and sublimation, partitions from air and from blood into biological tissues, and skin permeability from water. The descriptors themselves are not exceptional. Adamantane itself has a rather low dipolarity, zero hydrogen bond acidity and a very low hydrogen bond basicity, in common with other multicyclic aliphatic compounds. These lead to adamantane being a very hydrophobic compound, as is evident from our estimated water-octanol partition coefficient.

The assessment of intramolecular hydrogen bonding in ortho-substituted anilines by an NMR method

We describe the Δlog P method for the assessment of intramolecular hydrogen bonds (IMHBs), and show that it is not a very general method of distinguishing between molecules in which there is an IMHB and molecules in which there is no IMHB. The ‘double’ Δlog P method of Shalaeva et al. is a much more reliable method for the assessment of IMHB but requires the synthesis of a model compound and the determination of no less than four water-solvent partition coefficients. In addition, it is difficult to apply to compounds that contain more than one hydrogen bond acidic group capable of IMHB. We then describe our NMR method of assessing IMHB, based on 1H NMR chemical shifts in solvents DMSO and CDCl3. We have determined 1H NMR chemical shifts for a number of ortho-substituted anilines and show that the only compound we have studied that forms an IMHB is methyl 2-methylaminobenzoate though there is no IMHB present in methyl 2-aminobenzoate. This apparently anomalous result is supported by both MM and ab initio calculations.The NMR method is much simpler and less time consuming than other methods for the assessment of IMHB. It provides a quantitative assessment of IMHB and can be applied to molecules with more than one hydrogen bond acidic group.

Determination of the hydrogen-bond acidity and basicity for un-dissociated hydrazoic acid, isocyanic acid and isothiocyanic acid

Literature values of water-solvent partition coefficients for un-dissociated hydrazoic acid, isocyanic acid and isothiocyanic acid have been used to obtain properties or ‘descriptors’ of the three acids. These descriptors can be used in linear free energy relationships that we have already obtained to predict further values of partition coefficients for a very large number of solvents. The descriptors include the hydrogen-bond acidity, A, and hydrogen-bond basicity, B, of the acids. We have also used literature values of equilibrium constants for 1:1 complexation of the acids with a number of bases to obtain values of our 1:1 hydrogen-bond acidity, αH2. There is very good agreement between the two sets of hydrogen-bond acidity. Isothiocyanic acid is a strong hydrogen-bond acid about the same as chloroacetic acid, isocyanic acid is about as strong as acetic acid and hydrogen azide is a moderate-to-weak hydrogen bond acid. The hydrogen-bond basicity of the three acids is very small.

A new method for the determination of Henry's law constants (air-water-partition coefficients)

We set out a method for the determination of Henry's Law constants or air-water-partition coefficients at 298.15 K that does not require a knowledge of either the vapor pressure or the aqueous solubility of a compound, but is derived from experimental solubilities in non-aqueous solvents, or from experimental water-solvent partition coefficients. The method is applicable to a wide range of compound, from volatile organic compounds to very hydrophobic compounds with low vapor pressures and low aqueous solubilities. Using equations that we have developed previously, our obtained air-water partition coefficients at 298.15 K can be corrected to any temperature between 273 K and 573 K.

In Vitro Mechanistic Study of the Distribution of Lascufloxacin into Epithelial Lining Fluid.

The present study aimed to clarify the mechanism underlying the high distribution of lascufloxacin in epithelial lining fluid (ELF). Involvement of transporters was examined by transcellular transport across Calu-3 and transporter-overexpressing cells; the binding of lascufloxacin to ELF components was examined by an organic solvent-water partitioning system that employed pulmonary surfactant and phospholipids. Transcellular transport across the transporter-overexpressing cells indicated lascufloxacin to be a substrate of both P-glycoprotein (P-gp) and breast cancer resistance protein (BCRP); therefore, its transport across Calu-3 cells was inhibited by P-gp and BCRP inhibitors. However, permeability and efflux ratios of lascufloxacin were similar to those of the other quinolones with relatively low ELF distribution, indicating the existence of another mechanism for lascufloxacin distribution in ELF. Amongst pulmonary surfactants, which are a primary component of ELF, lascufloxacin preferentially bound to phosphatidylserine (PhS) from several phospholipids, and the binding was significantly greater than that for other quinolones. This binding was saturable with two apparent classes of binding sites and inhibited by some weakly basic drugs, indicating the presence of an ionic bond. In conclusion, the results of this study suggest that the binding of lascufloxacin to PhS in the pulmonary surfactant is the major mechanism of the high distribution of lascufloxacin in the ELF.

Solvation Descriptors for Zwitterionic α-Aminoacids; Estimation of Water-Solvent Partition Coefficients, Solubilities, and Hydrogen-Bond Acidity and Hydrogen-Bond Basicity.

The literature data on solubilities and water–solvent partition coefficients have been used to obtain properties or “Absolv descriptors” for zwitterionic α-aminoacids: glycine, α-alanine (α-aminopropanoic acid), α-aminobutanoic acid, norvaline (α-aminopentanoic acid), norleucine (α-aminohexanoic acid), valine (α-amino-3-methylbutanoic acid), leucine (α-amino-4-methylpentanoic acid), and α-phenylalanine. Together with equations that we have previously constructed, these descriptors can be used to estimate further solubilities and partition coefficients in a variety of organic solvents and in water–methanol and water–ethanol mixtures. It is shown that equations for neutral solutes are inadequate for the description of solubilities and partition coefficients for these α-aminoacids, and our equations developed for use with both neutral and ionic solutes must be used. The Absolv descriptors include those for hydrogen-bond acidity, A, and hydrogen-bond basicity, B. We find that both of these descriptors are far smaller in value than those for compounds that contain the corresponding ionic groups. Thus, A for α-alanine is 0.28, but A for the ethylammonium cation is 1.31; B for α-alanine is 0.83, and yet B for the acetate anion is no less than 2.93. The additional descriptors that we developed for equations that involve ions, J+ and J–, are very significant for the α-aminoacids, although numerically smaller than for ionic species such as EtNH3+ and CH3CO2–.

Descriptors for the α,ω-dicarboxylic acids from oxalic acid to sebacic acid

We use literature data on solubilities in water and solvents, water-solvent partition coefficients and vapor pressures at 298 K to obtain Abraham (Absolv) descriptors for the neutral, unionized, forms of oxalic acid, malonic acid, glutaric acid, succinic acid, pimelic acid, suberic acid, azaleic acid and sebacic acid. With a knowledge of equations for partition coefficients for numerous water-wet solvents, water to dry solvents and water to ionic liquid systems, we can then predict partition coefficients of the various acids in a huge number of systems. We set out the solubility of the acids in water and these, together with estimated partition coefficients yield solubilities of the acids into the corresponding solvents. The descriptors themselves show that the acids are dipolar, are strong hydrogen bond acids and reasonably strong hydrogen bond bases, although their hydrogen bond acidity is less than expected for compounds containing two (independent) carboxylic acid groups.

Estimating system parameters for solvent-water and plant cuticle-water using quantum chemically estimated Abraham solute parameters.

Polyparameter Linear Free Energy Relationships (pp-LFERs) using Abraham system parameters have many useful applications. However, developing the Abraham system parameters depends on the availability and quality of the Abraham solute parameters. Using Quantum Chemically estimated Abraham solute Parameters (QCAP) is shown to produce pp-LFERs that have lower root mean square errors (RMSEs) of predictions for solvent-water partition coefficients than parameters that are estimated using other presently available methods. pp-LFERs system parameters are estimated for solvent-water, plant cuticle-water systems, and for novel compounds using QCAP solute parameters and experimental partition coefficients. Refitting the system parameter improves the calculation accuracy and eliminates the bias. Refitted models for solvent-water partition coefficients using QCAP solute parameters give better results (RMSE = 0.278 to 0.506 log units for 24 systems) than those based on ABSOLV (0.326 to 0.618) and QSPR (0.294 to 0.700) solute parameters. For munition constituents and munition-like compounds not included in the calibration of the refitted model, QCAP solute parameters produce pp-LFER models with much lower RMSEs for solvent-water partition coefficients (RMSE = 0.734 and 0.664 for original and refitted model, respectively) than ABSOLV (4.46 and 5.98) and QSPR (2.838 and 2.723). Refitting plant cuticle-water pp-LFER including munition constituents using QCAP solute parameters also results in lower RMSE (RMSE = 0.386) than that using ABSOLV (0.778) and QSPR (0.512) solute parameters. Therefore, for fitting a model in situations for which experimental data exist and system parameters can be re-estimated, or for which system parameters do not exist and need to be developed, QCAP is the quantum chemical method of choice.

Quantum Chemically Estimated Abraham Solute Parameters Using Multiple Solvent-Water Partition Coefficients and Molecular Polarizability.

Polyparameter Linear Free Energy Relationships (pp-LFERs), also called Linear Solvation Energy Relationships (LSERs), are used to predict many environmentally significant properties of chemicals. A method is presented for computing the necessary chemical parameters, the Abraham parameters (AP), used by many pp-LFERs. It employs quantum chemical calculations and uses only the chemical's molecular structure. The method computes the Abraham E parameter using density functional theory computed molecular polarizability and the Clausius-Mossotti equation relating the index refraction to the molecular polarizability, estimates the Abraham V as the COSMO calculated molecular volume, and computes the remaining AP S, A, and B jointly with a multiple linear regression using sixty-five solvent-water partition coefficients computed using the quantum mechanical COSMO-SAC solvation model. These solute parameters, referred to as Quantum Chemically estimated Abraham Parameters (QCAP), are further adjusted by fitting to experimentally based APs using QCAP parameters as the independent variables so that they are compatible with existing Abraham pp-LFERs. QCAP and adjusted QCAP for 1827 neutral chemicals are included. For 24 solvent-water systems including octanol-water, predicted log solvent-water partition coefficients using adjusted QCAP have the smallest root-mean-square errors (RMSEs, 0.314-0.602) compared to predictions made using APs estimated using the molecular fragment based method ABSOLV (0.45-0.716). For munition and munition-like compounds, adjusted QCAP has much lower RMSE (0.860) than does ABSOLV (4.45) which essentially fails for these compounds.

Gas-solvent and water-solvent partition of trans-stilbene at 298 K

We start with literature data on gas to solvent partition of trans-stilbene, as log K, for 52 solvents. The log K values can be converted into corresponding water to solvent partitions, as log P, through the stilbene gas to water partition coefficient, and lead, with four additional equations, to a set of 108 simultaneous equations that can be solved to yield descriptors for stilbene with a standard deviation of no more than 0.096 log units. This represents statistically the best analysis of stilbene solubility to date. The stilbene descriptors can be combined with equation coefficients to give predictions of stilbene gas to solvent partition and water to solvent partition for numerous dry solvents, wet (water saturated) solvents and ionic liquid solvents, as well as predictions of environmental and biological processes involving styrene.

Predicting solvent-water partitioning of charged organic species using quantum-chemically estimated Abraham pp-LFER solute parameters

Methods for obtaining accurate predictions of solvent-water partitioning for neutral organic chemicals (e.g., Kow) are well-established. However, methods that provide comparable accuracy are not available for predicting the solvent-water partitioning of ionic species. Previous methods for addressing charge contributions to solvent-water partitioning rely on charged solute descriptors which are obtained from regressions to neutral species descriptors as well as charged descriptors which are specific to unique charge-functionalities and structural moieties. This paper presents a method for obtaining Abraham poly-parameter linear free energy relationship (pp-LFER) descriptors using quantum chemical calculations and molecular structure, only. The method utilizes a large number of solvent-water systems to overcome large errors in individual quantum chemical computations of ionic solvent-water partition coefficients. The result is a single set of quantum-chemically estimated Abraham solute parameters (QCAP) which are solvent-independent, and can be used to predict the solvent-water partitioning of ionic species.Predictions of solvent-water partition coefficients for ionic species using quantum-chemically estimated Abraham parameters (QCAPs) are shown to provide improved accuracy compared over both existing Absolv-estimated Abraham solute parameters (AAP) as well as direct a priori quantum chemical (QC) calculations for partitioning of anionic solutes in 4 organic solvent-water systems (RMS = 0.740, 2.48 and 0.426 for the Absolv, QC and QCAP methods, respectively). For quaternary amine cations in the octanol-water system the RMS errors of the solvent-water partition coefficients were larger and similar between the two Abraham models (RMSE = 0.997 and 1.16, for the AAP and QCAP methods, respectively). Both methods showed significant improvement over direct QC calculations (RMSE = 2.82).

Experimental determination of solvent-water partition coefficients and Abraham parameters for munition constituents

There is concern about the environmental fate and effects of munition constituents (MCs). Polyparameter linear free energy relationships (pp-LFERs) that employ Abraham solute parameters can aid in evaluating the risk of MCs to the environment. However, poor predictions using pp-LFERs and ABSOLV estimated Abraham solute parameters are found for some key physico-chemical properties. In this work, the Abraham solute parameters are determined using experimental partition coefficients in various solvent-water systems. The compounds investigated include hexahydro-1,3,5-trinitro-1,3,5-triazacyclohexane (RDX), octahydro-1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX), hexahydro-1-nitroso-3,5-dinitro-1,3,5-triazine (MNX), hexahydro-1,3,5-trinitroso-1,3,5-triazine (TNX), hexahydro-1,3-dinitroso-5- nitro-1,3,5-triazine (DNX), 2,4,6-trinitrotoluene (TNT), 1,3,5-trinitrobenzene (TNB), and 4-nitroanisole. The solvents in the solvent-water systems are hexane, dichloromethane, trichloromethane, octanol, and toluene. The only available reported solvent-water partition coefficients are for octanol-water for some of the investigated compounds and they are in good agreement with the experimental measurements from this study. Solvent-water partition coefficients fitted using experimentally derived solute parameters from this study have significantly smaller root mean square errors (RMSE = 0.38) than predictions using ABSOLV estimated solute parameters (RMSE = 3.56) for the investigated compounds. Additionally, the predictions for various physico-chemical properties using the experimentally derived solute parameters agree with available literature reported values with prediction errors within 0.79 log units except for water solubility of RDX and HMX with errors of 1.48 and 2.16 log units respectively. However, predictions using ABSOLV estimated solute parameters have larger prediction errors of up to 7.68 log units. This large discrepancy is probably due to the missing R2NNO2 and R2NNO2 functional groups in the ABSOLV fragment database.

Equations for water-triolein partition coefficients for neutral species; comparison with other water-solvent partitions, and environmental and toxicological processes

Linear free energy relationships, LFERs, have been constructed for water-triolein partition coefficients for neutral species. It is shown that separate equations are required for wet and dry triolein. From a comparison of the equation coefficients for water-wet triolein with those for 52 other water-solvent systems it is shown that there is little correspondence between triolein and any of the 52 other solvents - only the water-isopropyl myristate system is close to the water-wet triolen system. A comparison of equation coefficients for the water-wet triolein system with LFER coefficients of 16 environmentally important processes shows that wet triolein is not a suitable model for any of the processes, although a number of other water-solvent systems are possible models for some of the environmental processes. A comparison of LFER coefficients with those of 17 aqueous toxicological processes reveals that most of the water-solvent systems, including water-wet triolein, will be poor models for any of the toxicological systems, but the water-lower alcohol systems show promise as models for a number of the toxicological systems. Our method of comparison of coefficients for LFERs that have exactly the same independent variables can be extended to various other types of system.

A linear free energy analysis of PAMPA models for biological systems

Equation coefficients in the Abraham linear free energy relationship have been obtained for a number of parallel artificial membrane permeability assays, PAMPA, and have been compared to equation coefficients for several biological processes and a large number of water-solvent partitions. From a principal component analysis of the equation coefficients we show that one particular PAMPA (Wohnsland and Faller, 2001) is a good model for human intestinal absorption, but we could find no PAMPA that was a good model for rat brain permeation, mouse brain permeation, blood–brain distribution, skin permeation or water–skin partition, although models that include PAMPA and other descriptors may be useful. The water–octanol system is a poor model for all the biological process we have studied, but water-nonpolar solvent partitions are good models for permeation through egg lecithin membranes.

Using water-solvent systems to estimate in vivo blood-tissue partition coefficients.

BackgroundBlood–tissue partition coefficients indicate how a chemical will distribute throughout the body and are an important part of any pharmacokinetic study. They can be used to assess potential toxicological effects from exposure to chemicals and the efficacy of potential novel drugs designed to target certain organs or the central nervous system. In vivo measurement of blood–tissue partition coefficients is often complicated, time-consuming, and relatively expensive, so developing in vitro systems that approximate in vivo ones is desirable. We have determined such systems for tissues such as brain, muscle, liver, lung, kidney, heart, skin, and fat.ResultsSeveral good (p < 0.05) blood–tissue partition coefficient models were developed using a single water–solvent system. These include blood–brain, blood–lung, blood–heart, blood–fat, blood–skin, water–skin, and skin permeation. Many of these partition coefficients have multiple water–solvent systems that can be used as models. Several solvents—methylcyclohexane, 1,9-decadiene, and 2,2,2-trifluoroethanol—were common to multiple models and thus a single measurement can be used to estimate multiple blood–tissue partition coefficients. A few blood–tissue systems require a combination of two water–solvent partition coefficient measurements to model well (p < 0.01), namely: blood–muscle: chloroform and dibutyl ether, blood–liver: N-methyl-2-piperidone and ethanol/water (60:40) volume, and blood–kidney: DMSO and ethanol/water (20:80) volume.ConclusionIn vivo blood–tissue partition coefficients can be easily estimated through water–solvent partition coefficient measurements.Graphical abstract:Predicted blood-brain barrier partition coefficients coloured by measured log BB valueElectronic supplementary materialThe online version of this article (doi:10.1186/s13065-015-0134-z) contains supplementary material, which is available to authorized users.

Determination of Abraham model solute descriptors for the monomeric and dimeric forms of trans-cinnamic acid using measured solubilities from the Open Notebook Science Challenge.

BackgroundCalculating Abraham descriptors from solubility values requires that the solute have the same form when dissolved in all solvents. However, carboxylic acids can form dimers when dissolved in non-polar solvents. For such compounds Abraham descriptors can be calculated for both the monomeric and dimeric forms by treating the polar and non-polar systems separately. We illustrate the method of how this can be done by calculating the Abraham descriptors for both the monomeric and dimeric forms of trans-cinnamic acid, the first time that descriptors for a carboxylic acid dimer have been obtained.ResultsAbraham descriptors were calculated for the monomeric form of trans-cinnamic acid using experimental solubility measurements in polar solvents from the Open Notebook Science Challenge together with a number of water-solvent partition coefficients from the literature. Similarly, experimental solubility measurements in non-polar solvents were used to determine Abraham descriptors for the trans-cinnamic acid dimer.ConclusionAbraham descriptors were calculated for both the monomeric and dimeric forms of trans-cinnamic acid. This allows for the prediction of further solubilities of trans-cinnamic acid in both polar and non-polar solvents with an error of about 0.10 log units.Graphical abstractMolar concentration of trans-cinnamic acid in various polar and non-polar solvents.