Systematic search for surrogate chromatographic models of biopartitioning processes

A systematic approach for identifing surrogate chromatographic models for biopartitioning processes is described. The method is based on a comparison of the system constant ratios of the solvation parameter model for biopartitioning processes and a database of system constant ratios for reversed-phase liquid chromatographic and micellar electrokinetic chromatographic systems compiled from literature sources. An acceptance filter of < or = 0.2 is applied for each difference in system constant ratio for the compared systems to provide a reasonable probability of success without outputting too many systems with limited predictive properties. Surrogate chromatographic models identified for the non-specific toxicity of neutral organic compounds to the fathead minnow and the soil-water distribution constant are tested by construction of a correlation model for the characteristic property of the biological process and the chromatographic retention factors for a structurally varied group of compounds. Although these models are not the best that could be obtained based on ranking of the differences in system constant ratios the predictive ability of the correlation models is suitable for typical applications and similar to the accepted uncertainty in the measurements of the biological property. Retention factors on the immobilized artificial membrane column (IAM PC DD 2) with 10% (v/v) methanol-water as mobile phase are able to estimate non-specific toxicity to the fathead minnow with a standard error (SE) of 0.22 log units and coefficient of determination (r2) of 0.97 for 31 compounds. Retention factors on a Bakerbond DIOL column with 20% (v/v) acetonitrile-water as mobile phase are able to estimate the soil-water distribution constant with an SE of 0.38 log units and r2 = 0.88 for 59 compounds. Other potential surrogate chromatographic models are identified for non-specific toxicity to the guppy, tadpole, Vibrio fischeri and Terahymena pyriformis as well as the plant cuticle matrix-water distribution constant. On the other hand reversed-phase chromatographic systems seem poorly suited for estimating intestinal absorption and the blood-brain distribution constant.     

Retention characteristics of porous graphitic carbon in subcritical fluid chromatography with carbon dioxide-methanol mobile phases

Numerous relationships usually used in high-performance liquid chromatography (HPLC) for describing the retention on porous graphitic carbon (PGC) have been applied in subcritical fluid chromatography, with CO2-methanol mobile phases. As reported in HPLC, octanol-water partition coefficient failed to fit the retention, whereas satisfactory results were obtained with the sum of partial negative charges. A better fit was reached by using the solvation parameter model, allowing a better understanding of the interactions developed between the solute, the stationary and the mobile phases. Results show that the dominant contribution to retention was given by the polarizability (E) and the volume (V), while the hydrogen-bond basicity (B) was not selected in the retention model, whatever the methanol content. The increase in methanol percentage favours the retention decrease, mainly through the volume for hydrophobic compounds, and through the hydrogen-bond acidity for polar compounds.     

Separation characteristics of phenyl-containing stationary phases for gas chromatography based on silarylene-siloxane copolymer chemistries

The solvation parameter model is used to characterize the retention properties of five open-tubular column stationary phases (ZB-5 ms, DB-5 ms, DB-XLB, DB-17 ms, and DB-35 ms) based on silarylene-siloxane copolymer chemistries at five equally spaced temperatures over the range 60-140 degrees C. System constant differences and regression models for varied compounds are used to establish the selectivity equivalence of the silarylene-siloxane copolymer stationary phases and to compare their separation characteristics with poly(dimethyldiphenylsiloxane) stationary phases containing a nominally similar concentration of phenyl groups. These studies demonstrate that ZB-5 ms and DB-5 ms are selectivity equivalent. DB-XLB is significantly more dipolar and polarizable than DB-5 ms. In general terms, the silarylenesiloxane copolymer stationary phases are slightly less cohesive and more dipolar and polarizable with similar hydrogen-bond basicity to the poly(dimethyldiphenylsiloxane) stationary phases they were designed to replace. None of the silarylenesiloxane copolymer or poly(dimethyldiphenylsiloxane) stationary phases are hydrogen-bond acidic. Selectivity differences between the two types of stationary phase are temperature dependent and tend to be smaller at higher temperatures within the temperature range studied. Consequently, selectivity differences cannot be globalized without reference to the temperature for the comparison.     

Chromatographic models for the sorption of neutral organic compounds by soil from water and air.

The solvation parameter model is used to construct models for the estimation of the soil-water and soil-air distribution constants and to characterize the contribution of fundamental intermolecular interactions to the underlying sorption processes. Wet soil is shown to be quite cohesive and polar but relatively non-selective for dipole-type, lone-pair electron and hydrogen-bond interactions. Using a comparison of system constant ratios chromatographic systems employing reversed-phase liquid chromatography on polar bonded phases are shown to provide suitable models for estimating soil-water distribution constants. No suitable gas chromatographic models were found for the soil-air distribution constant but the requirements for such a system are indicated. Models are also provided for adsorption at the air-water interface. Estimation methods based on either the solvation parameter model or chromatographic model reproduce experimental distribution constants for a wide variety of compounds with a similar error (0.2-0.3 log units) to that expected in the experimental data.     

有机化合物脂水分配系数和溶解度的计算方法

有机化合物的脂水分配系数和溶解度在药物化学以及环境化学的研究中是十分重要的物理化学性质,在许多涉及有机化合物在生物体内的吸收、转运、代谢以及在环境中迁移等过程的定量构效关系研究(QSAR)中发挥着不可替代的作用.目前在实践中应用较多的计算有机化合物脂水分配系数和溶解度的理论方法主要有片段加合法和基于描述符的方法.本文总结了这两大类方法的优缺点以及在该领域中未来可能的发展方向.     

Atomic physicochemical parameters for three dimensional structure directed quantitative structure-activity relationships III: Modeling hydrophobic interactions

In an earlier article⁸ the need was demonstrated for atomic physicochemical properties for three dimensional structure directed quantitative structure-activity relationships, and it was shown how atomic parameters can be developed for successfully evaluating the molecular octanol-water partition coefficient, which is a measure of hydrophobicity. In this work we report more refined atomic values of octanol-water partition coefficients derived from nearly twice the number of compounds. Carbon, hydrogen, oxygen, nitrogen, sulfur and halogens are divided into 110 atom types of which 94 atomic values are evaluated from 830 molecules by least squares. These values gave a standard deviation of 0.470 and a correlation coefficient of 0.931. These parameters predicted the octanol-water partition coefficient of 125 compounds with a standard deviation of 0.520 and a correlation coefficient of 0.870. There is only a correlation coefficient of 0.432 between the atomic octanol-water partition coefficients and the atomic contributions to molar refractivity over the 93 atom types used for both the properties. This suggests that both parameters can be used simultaneously to model intermolecular interactions. We evaluated the CNDO/2 gross atomic charge distribution over several molecules to check the validity of our classification. We found that the charge density on the heteroatoms in conjugated systems is strongly affected by the presence of similar atoms in the conjugation which suggests it should be incorporated as a separate parameter in evaluating the partition coefficient.     

Retention properties of acetone-water mobile phases on a biphenylsiloxane-bonded silica stationary phase in reversed-phase liquid chromatography

The solvation parameter model system constants and retention factors were used to interpret retention properties of 39 calibration compounds on a biphenylsiloxane-bonded stationary phase (Kinetex biphenyl) for acetone-water binary mobile phase systems containing 30–70% v/v. Variation in system constants, phase ratios, and retention factors of acetone-water binary mobile phases systems were compared with more commonly used acetonitrile and methanol mobile phase systems. Retention properties of acetone mobile phases on a Kinetex biphenyl column were more similar to that of acetonitrile than methanol mobile phases except with respect to selectivity equivalency. Importantly, selectivity differences arising between acetone and acetonitrile systems (the lower hydrogen-bond basicity of acetone-water mobile phases and differences in hydrogen-bond acidity, cavity formation and dispersion interactions) could be exploited in reversed-phase liquid chromatography method development on a Kinetex biphenyl stationary phase.     

Evaluation of the separation characteristics of application-specific (volatile organic compounds) open-tubular columns for gas chromatography

The solvation parameter model is used to characterize the separation characteristics of two application-specific open-tubular columns (Rtx-Volatiles and Rtx-VGC) and a general purpose column for the separation of volatile organic compounds (DB-WAXetr) at five equally spaced temperatures over the range 60-140 degrees C. System constant differences and retention factor correlation plots are then used to determine selectivity differences between the above columns and their closest neighbors in a large database of system constants and retention factors for forty-four open-tubular columns. The Rtx-Volatiles column is shown to have separation characteristics predicted for a poly(dimethyldiphenylsiloxane) stationary phase containing about 16% diphenylsiloxane monomer. The Rtx-VGC column has separation properties similar to the poly(cyanopropylphenyldimethylsiloxane) stationary phase containing 14% cyanopropylphenylsiloxane monomer DB-1701 for non-polar and dipolar/polarizable compounds but significantly different characteristics for the separation of hydrogen-bond acids. For all practical purposes the DB-WAXetr column is shown to be selectivity equivalent to poly(ethylene glycol) columns prepared using different chemistries for bonding and immobilizing the stationary phase. Principal component analysis and cluster analysis are then used to classify the system constants for the above columns and a sub-database of eleven open-tubular columns (DB-1, HP-5, DB-VRX, Rtx-20, DB-35, Rtx-50, Rtx-65, DB-1301, DB-1701, DB-200, and DB-624) commonly used for the separation of volatile organic compounds. A rationale basis for column selection based on differences in intermolecular interactions is presented as an aid to method development for the separation of volatile organic compounds.     

Structure-toxicity relationships for aliphatic compounds encompassing a variety of mechanisms of toxic action to Vibrio fischeri

QSARs based upon the logarithm of the octanol-water partition coefficient, log P, and energy of the lowest unoccupied molecular orbital, ELUMO were developed to model the toxicity of aliphatic compounds to the marine bacterium Vibrio fischeri. Statistically robust, hydrophobic-dependent QSARs were found for chloroalcohols and haloacetonitriles. Modelling of the toxicity of the haloesters and the diones required the use of terms to describe both hydrophobicity and electrophilicity. The differences in intercepts, slopes, and fit of these models suggest different electrophilic mechanisms occur between classes, as well as within the diones and haloesters. In order to model globally the toxicity of aliphatic compounds to V. fischeri, all the data determined in this study were combined with those determined previously for alkanones, alkanals, and alkenals. A highly predictive two-parameter QSAR [pT15 = 0.760(log P) - 0.625(ELUMO) - 0.466; n = 63, s = 0.462, r2 = 0.846, F = 171, Pr > F = 0.0001] was developed for the combined data that models across classes and is independent of mechanisms of action. The toxicity of these compounds to V. fischeri compares well to the toxicity (50% population growth inhibition) to the ciliate Tetrahymena pyriformis (r2 = 0.850).     

Selectivity equivalence of poly(ethylene glycol) stationary phases for gas chromatography

The solvation parameter model is used to study differences in selectivity for poly(ethylene glycol) stationary phases for packed column (Carbowax 20M) and fused-silica, open-tubular column (HP-20M, AT-Wax, HP-INNOWax and DB-FFAP) gas chromatography. All phases are dipolar, strongly hydrogen-bond basic with no hydrogen-bond acidity and of moderate cohesion. No two phases are exactly alike, however, and selectivity differences identified with cavity formation and dispersion interactions, n- and pi-electron pair interactions, dipole-type interactions and hydrogen-bond interactions are quantified by differences in the system constants at a fixed temperature where retention occurs solely by gas-liquid partitioning. The system constants vary linearly with temperature over the range 60-140 degrees C (except for n- and pi-electron pair interactions which are temperature invariant) facilitating a general comparison of the importance of temperature on selectivity differences for compared phases. From a mechanistic point of view it is demonstrated that selectivity differences can result from chemical differences between the poly(ethylene glycol) stationary phases and from differences in the relative contribution of interfacial adsorption to the retention mechanism. The latter depends on both system properties and solute characteristics.     

Selectivity equivalence of two poly(methylphenylsiloxane) open-tubular columns prepared with different deactivation techniques for gas chromatography

The solvation parameter model is used to characterize the retention properties of a poly(methylphenylsiloxane) column Rxi-50 over the temperature range 60-240 degrees C. The smooth variation of the system constants with temperature affords a general picture of how the relative importance of the different intermolecular interactions change with temperature. The system constants and retention factors for varied compounds are compared with those for Rtx-50 prepared with a similar stationary phase but using a different surface deactivation technique. The two columns are shown to be nearly selectivity equivalent. The Rtx-50 column is slightly more cohesive, dipolar/polarizable and hydrogen-bond basic than Rxi-50, while Rxi-50 is slightly more electron lone pair attractive and hydrogen-bond acidic. Only the difference in hydrogen-bond acidity can be identified with some certainty as related to the difference in deactivation processes. For compounds with a separation greater than 0.2 retention factor units on Rtx-50, it should be relatively straightforward to achieve an acceptable separation for the same compounds on Rxi-50.     

The application of the dielectric constant measurements to determination of the hydrogen-bond energy of hydrophobic organic compounds

Infrared absorption spectroscopy belongs to these analytical methods, which have wide application in estimating hydrogen-bond energies (1). In the present paper an effort was undertaken to apply the dielectric constant measurements to evaluate the mean hydrogen-bond energy levels of some selected homologous groups of chemical compounds. The chemical compounds examined in our experiment were higher fatty alcohols and acids, which showed hydrophobic properties.     

Retention characteristics of an immobilized artificial membrane column in reversed-phase liquid chromatography

Retention for a varied group of compounds on an immobilized artificial membrane column (IAM PC DD2) with a methanol-water mobile phase is shown to fit a second-order model for the retention factor (log k) as a function of the volume fraction of organic solvent. The numerical value of the intercept obtained by linear extrapolation to zero organic solvent (log k(w)) is shown to depend on the range of mobile phase composition used for the extrapolation. Each series of intercepts so obtained represents a different hypothetical distribution system as identified by the system constants of the solvation parameter model. Although a linear model is a poor fit for isocratic retention data, the linear solvent strength gradient model provides a reasonable estimate of isocratic retention factor values that are (slightly) larger than experimental values, but provide the same chemical information for the system. These preliminary results suggest that gradient elution may prove to be a rapid and useful method for creating system maps for column characterization and method development. In this work a system map is provided for methanol-water compositions from 0 to 60% (v/v) methanol and additional system constants for acetonitrile-water compositions containing 20 and 30% (v/v) acetonitrile. It is shown that the main factors contributing to retention on the IAM PC DD2 column are favorable cavity formation and dispersion interactions, electron lone pair interactions and the hydrogen-bond basicity of the sorbent. The latter feature more than any other distinguishes the IAM column from conventional chemically bonded phases. Interactions of a dipole-type (weakly) and inability to compete with the mobile phase as a hydrogen-bond acid reduce retention. A comparison of system constant ratios is used to demonstrate that the retention properties of the IAM column are not easily duplicated by conventional chemically bonded phases. The retention characteristics of the IAM column, however, are strongly correlated with the retention properties of pseudostationary phases used for micellar electrokinetic chromatography, which provide a suitable alternative to IAM columns for physical property estimations. By the same comparative method it is shown that retention on the IAM column possesses some similarity to biomembrane absorption processes, allowing suitable correlation models to be developed for the estimation of certain biopartitioning properties.