Revised descriptors for polycyclic aromatic and related hydrocarbons for the prediction of environmental properties using the solvation parameter model

Guidelines for descriptor assignments for the solvation parameter model by separation techniques

The solvation parameter model uses six compound descriptors to model equilibrium properties in biphasic systems formally defined as excess molar refraction, E, dipolarity/polarizability, S, overall hydrogen-bond acidity, A, overall hydrogen-bond basicity, B, McGowan's characteristic volume, V, and the gas-liquid partition constant on hexadecane at 25 °C, L. The V descriptor can be assigned from structure and the E descriptor for compounds liquid at 20 °C can be calculated from its refractive index and characteristic volume. The E descriptor for compounds solid at 20 °C and the S, A, B, and L descriptors are assigned from experimental properties traditionally obtained by chromatographic, liquid-liquid partition, and solubility measurements. Here I report an efficient experimental design using the Solver method for the accurate assignment of descriptors for neutral compounds that simultaneously minimizes laboratory resources. This multi-technique approach requires 3 retention factor measurements in a 60 °C temperature range per compound on four columns by gas chromatography, 3 retention factor measurements in a 30 % (v/v) acetonitrile composition range per compound on two columns by reversed-phase liquid chromatography, and eight partition constant measurements by liquid-liquid partition in totally organic and aqueous biphasic systems for a total of 26 experimental measurements. The accuracy of the descriptor assignments was validated by comparison with the values in the Wayne State University (WSU) descriptor database taken as the best estimate of the true descriptor values. The E, S, A, B and L descriptors were assigned simultaneously by the Solver method using the above approach without significant bias and with an average absolute deviation (AAD) of 0.054, 0.018, 0.015, 0.013, and 0.040, respectively, compared with the WSU database values, corresponding to a relative absolute average deviation in percent (RAAD) of 7.2, 1.9, 3.6, 5.1, and 0.84 %, respectively, for 32 varied compounds. This streamlined approach represents a significant improvement on earlier single-technique approaches used as the starting point for the development of the multi-technique approach. For compounds of variable hydrogen-bond basicity modifications to the multi-technique approach were implemented while maintaining the same number of experimental measurements. Acceptable descriptor assignments for B/B° were obtained for compounds liquid at 20 °C for which the E descriptor was available by calculation. For solid compounds at 20 °C the E and B/B° descriptors are restricted to qualitative application where approximate values may be acceptable.

Evaluation of the Goss-modified solvation parameter model for the characterization of biphasic systems and descriptor assignments

The solvation parameter model uses six descriptors identified as excess molar refraction, E, dipolarity/polarizability, S, overall hydrogen-bond acidity, A, overall hydrogen-bond basicity, B, McGowan's characteristic volume, V, and the gas-liquid partition constant on hexadecane at 25 °C, L to model the distribution of neutral compounds in biphasic systems. Abraham's version of this model uses all six descriptors with two separate linear free energy relationship models for the transfer of compounds from a gas phase to a condensed phase and between condensed phases. Goss proposed a modification to this model that uses a single calibration model regardless of the physical state for each phase and five of the descriptors employed in Abraham's model (E descriptor is eliminated). The capability of Abraham's model and the Goss-modified model to characterize the contribution of intermolecular interaction to retention for gas and reversed-phase liquid chromatographic systems and distribution in liquid-liquid partition systems is evaluated using the WSU compound descriptor database. These more accurate values for the Abraham descriptors have not been utilized previously for the evaluation of the Goss-modified model and should be more capable of discerning subtle differences in model performance. It is shown that model quality defined by statistical parameters favors Abraham's model over the Goss-modified model with differences in model quality greater for systems in which Abraham's model indicates a significant contribution from electron lone pair interactions and for systems in which one phase is a solvent containing perfluoroalkyl substituents. There is a small systematic difference for the terms describing the combined contributions of cavity formation and dispersion interactions and for interactions of a dipole-type. The contribution of hydrogen-bonding interactions is virtually identical for the two models. The model intercepts are generally different and potentially assigned to a larger contribution from lack-of-fit for the Goss-modified model. Although the Abraham model descriptors have been routinely employed for applications using the Goss-modified model the possibility that Goss-model specific descriptors should be employed was evaluated. Using the Solver method and Goss-model specific calibration models for chromatographic and liquid-liquid partition systems a new set of Goss-specific descriptors was calculated for 28 varied compounds. These descriptors show good statistical agreement with the Abraham descriptor values with an average deviation of 0.009, -0.003, -0.004, and -0.023, respectively, for the S, A, B, and L descriptors, corresponding to a relative absolute deviation in percent of 2.2 %, 3.9 %, 4.3 %, and 1.2 %, respectively.

Assessment of liquid-liquid partition for the assignment of descriptors for the solvation parameter model

The solvation parameter model uses five system independent descriptors to characterize compound properties defined as excess molar refraction, E, dipolarity/polarizability, S, hydrogen-bond acidity, A, hydrogen-bond basicity, B, and McGowan's characteristic volume, V, to model transfer properties between condensed phases. The V descriptor is assigned from structure. For compounds liquid at 20 °C the E descriptor can be assigned from the characteristic volume and its refractive index. The E descriptor for compounds solid at 20 °C and the S, A, and B descriptors are experimental properties traditionally assigned from chromatographic, liquid-liquid partition, and solubility measurements. In this report liquid-liquid partition constants in totally organic and aqueous biphasic systems are evaluated as a standalone technique for descriptor assignments. Using six totally organic biphasic systems the S, A, and B descriptors were assigned with an average absolute deviation (AAD) of about 0.04, 0.03, and 0.04, respectively, compared with the best estimate of the true descriptor values for 65 compounds. The E descriptor for compounds solid at 20 °C can only be estimated with an AAD of approximately 0.1. For six aqueous biphasic systems the B descriptor is assigned with a lower AAD of 0.028 and higher AAD of 0.08 and 0.05 for the S and A descriptors, respectively, than for the totally organic biphasic systems for compounds with a reliable value for the E descriptor. The preferred system for descriptor assignments utilizes both totally organic biphasic systems (heptane-1,1,1-trifluoroethanol, isopentyl ether-propylene carbonate, isopentyl ether-ethanolamine, heptane-ethylene glycol, heptane-formamide, and 1,2-dichloroethane-ethylene glycol) and aqueous biphasic systems (octanol-water, cyclohexane-water) with the possible substitution of some systems with alternative systems of similar selectivity. For 55 varied compounds this combination of eight organic and aqueous biphasic systems resulted in an AAD of approximately 0.03, 0.02, and 0.02 for the S, A, and B descriptors compared to the best estimate of the true descriptor value. For 30 compounds solid at 20 °C the AAD for the E descriptor of 0.11 is poorly assigned. The relative average absolute deviation in percent (RAAD) corresponds to 9.7 %, 3.1 %. 4.0 % and 8.3 % for E, S, A, and B, respectively, for the eight biphasic systems. Liquid-liquid partition is compared to reversed-phase liquid and gas chromatography as a standalone technique for descriptor assignments.

Determination of solvation parameter model compound descriptors by gas chromatography

The solvation parameter model uses five system independent descriptors to characterize compound properties defined as excess molar refraction, E, dipolarity/polarizability, S, hydrogen-bond acidity, A, hydrogen-bond basicity, B, and the gas-liquid partition constant at 25 °C on n-hexadecane, L, to model transfer properties in gas-condensed phase biphasic systems. The E descriptor for compounds liquid at 20 °C is available by calculation using a refractive index value while E for solid compounds at 20 °C and the S, A, B, and L descriptors are determined by experiment. As a single-technique approach, it is shown that with up to 20 retention factor measurements on four columns comprising a poly(siloxane) containing methyloctyl or dimethyldiphenylsiloxane monomers (SPB-Octyl or HP-5), a poly(siloxane) containing methyltrifluoropropylsiloxane monomers (Rtx-OPP or DB-210), a poly(siloxane) containing bis(cyanopropylsiloxane) monomers (HP-88 or SGE BPX-90), and a poly(ethylene glycol) stationary phase (DB-WAXetr or HP-INNOWAX) are suitable for assigning the S, A, and L descriptors. Using the descriptors in the updated WSU compound descriptor database as target values the average absolute error in the descriptor assignments for 52 varied compounds in the temperature range 60-140 °C was 0.072 for E, 0.016 for S, 0.008 for A, and 0.022 for L corresponding to 30 %, 3.5 %, and 0.6 % as a relative average absolute error for E, S, and L, respectively. For the higher temperature range of 160-240 °C and 34 varied compounds that are liquid at 20 °C the average absolute error for the S, A and L descriptors was 0.026, 0.020, and 0.031, respectively, with the largest relative average absolute error for S of 3.2 % (< 1 % for the L descriptor). For 35 varied compounds that are solid at 20 °C the relative absolute error for the E, S, A, and L descriptors in the higher temperature range was 0.068, 0.035, 0.020, and 0.020, respectively, with a relative average absolute error for E (6.5 %), S (3.5 %) and L (0.88 %). The S, A, and L descriptor can be accurately assigned on the four-column system over a wide temperature range. The E descriptor for solid compounds at 20 °C exhibits greater variability than desirable. The B descriptor cannot be assigned by the four-column system, which lack hydrogen-bond acid functional groups, and is only poorly assigned on the weak hydrogen-bond acid ionic liquid column SLB-IL100.

Determination of the hydrogen-bond basicity descriptor by reversed-phase liquid chromatography

Except for alkanes, most organic compounds are hydrogen-bond bases. The B° descriptor of the solvation parameter model provides a convenient measure of the effective (or summation) hydrogen-bond basicity of organic compounds. A fast and convenient method to assign the B° descriptor is required to support studies of hydrogen-bonding in separation systems. A two-column system with acetonitrile-water mobile phase compositions and the measurement of up to eleven isocratic retention factors is proposed for this purpose. Several reversed-phase column chemistries and mobile phases were evaluated with the two-column system consisting of a pentafluorophenylpropylsiloxane-bonded and octadecylsiloxane-bonded silica columns recommended for this purpose. To assess the accuracy of the method values for B° were taken from the Wayne State University (WSU) compound descriptor database, which were assigned using conventional multi-technique methods and large datasets. The two-column systems provided an unbiased assignment of B° with an average deviation of 0.008 and an average absolute deviation of 0.021 compared with the target value for 55 varied compounds. The two-column system is unsuitable for assigning the other descriptors used in the solvation parameter model and results in erroneous assignments of B° for nitrogen-containing compounds capable of electrostatic interactions on silica-based reversed-phase columns.