Modeling Aquatic Toxicity through Chromatographic Systems

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 descriptors for the principal flavor compounds of the cinnamons of commerce by gas chromatography and liquid-liquid partition

Descriptors for fourteen semivolatile organic compounds associated with the authenticity, botanical origin, and flavor potential of the cinnamons of commerce were determined using the Solver method and experimental retention factors determined by gas chromatography at several temperatures on a minimum of seven selectivity-selected, open-tubular columns and liquid-liquid partition constants in up to twenty totally organic biphasic systems. The six descriptors that encode the solvation properties of the compounds were used to predict water-gas, octanol-gas, and octanol-water partition constants commonly employed to assess environmental distribution properties. For octanol-water partition constants, log KOW, the predicted partition constants exhibited an average absolute deviation of 0.12 for log KOW experimental - log KOW predicted (n = 14). Soil-water, soil-air, urban aerosol-air, skin-water permeation, and non-specific toxicity to the fathead minnow were predicted for the same compounds to assess their potential environmental impact. The product terms of the solvation parameter model provide a useful insight into the contribution of individual intermolecular interactions to the distribution properties of the cinnamon compounds and their environmental impact.

Characterization of solute-solvent interactions in liquid chromatography systems: A fast method based on Abraham's linear solvation energy relationships

The Abraham's solvation parameter model, based on linear solvation energy relationships (LSER), allows the accurate characterization of the selectivity of chromatographic systems according to solute-solvent interactions (polarizability, dipolarity, hydrogen bonding, and cavity formation). However, this method, based on multilinear regression analysis, requires the measurement of the retention factors of a considerably high number of compounds, turning it into a time-consuming low throughput method. Simpler methods such as Tanaka's scheme are preferred. In the present work, the Abraham's model is revisited to develop a fast and reliable method, similar to the one proposed by Tanaka, for the characterization of columns employed in reversed-phase liquid chromatography and particularly in hydrophilic interaction liquid chromatography. For this purpose, pairs of compounds are carefully selected in order to have in common all molecular descriptors except for a specific one (for instance, similar molecular volume, dipolarity, polarizability, and hydrogen bonding basicity features, but different hydrogen bonding acidity). Thus, the selectivity factor of a single pair of test compounds can provide information regarding the extent of the dissimilar solute-solvent interactions and their influence on chromatographic retention. The proposed characterization method includes the determination of the column hold-up volume and Abraham's cavity term by means of the injection of four alkyl ketone homologues. Therefore, five chromatographic runs in a reversed-phase column (four pairs of test solutes and a mixture of four homologues) are enough to characterize the selectivity of a chromatographic system. Tanaka's method is also analyzed from the LSER point of view.

The effect of descriptor database selection on the physicochemical characterization and prediction of water-air, octanol-air and octanol-water partition constants using the solvation parameter model

The distribution of neutral compounds in biphasic separation systems can be described by the solvation parameter model using six solute properties, or descriptors. These descriptors (McGowan's characteristic volume, excess molar refraction, dipolarity/polarizability, hydrogen-bond acidity and basicity, and the gas-liquid partition constant on n-hexadecane at 298.15 K) are curated in two publicly accessible databases for hundreds (WSU compound descriptor database) or thousands (Abraham compound descriptor database). These databases were developed independently using different approaches resulting in descriptor values that vary for many compounds. Previously, it was shown that the two descriptor databases are not interchangeable, and the WSU descriptor database consistently demonstrated improved model performance for chromatographic systems where the uncertainty in the dependent variable was minimized by suitable quality control and calibration procedures. In this report we wish to evaluate whether the same conclusions are true for models with a dependent variable containing significant measurement uncertainty. To evaluate this hypothesis, we assembled databases for water-air, octanol-air, and octanol-water partition constants reported by multiple laboratories using various measurement methods. It was found that database selection has little effect on model quality or model predictive capability but significantly affects the assignment of the contribution of individual intermolecular interactions to the dependent variable. The latter information is database specific, and a quantitative comparison of system constants should be restricted to models using the same compound descriptor database.

Evaluation of the Ability of PAMPA Membranes to Emulate Biological Processes through the Abraham Solvation Parameter Model

Two parallel artificial membrane permeability assay (PAMPA) systems intended for emulating skin permeability have been characterized through the solvation parameter model of Abraham using multilinear regression analysis. The coefficients of the obtained equations have been compared to the ones already established for other PAMPA membranes using statistical tools. The results indicate that both skin membranes are similar to each other in their physicochemical properties. However, they are different from other PAMPA membranes (e.g., intestinal absorption and blood-brain PAMPAs), mainly in terms of hydrophobicity and hydrogen bonding properties. Next, all PAMPA membranes have been compared to relevant biological processes also characterized through the solvation parameter model. The results highlight that skin-PAMPA membranes are a very good choice to emulate skin permeability.

Using membrane-water partition coefficients in a critical membrane burden approach to aid the identification of neutral and ionizable chemicals that induce acute toxicity below narcosis levels

The risk assessment of thousands of chemicals used in our society benefits from adequate grouping of chemicals based on the mode and mechanism of toxic action (MoA). We measure the phospholipid membrane-water distribution ratio (D_MLW) using a chromatographic assay (IAM-HPLC) for 121 neutral and ionized organic chemicals and screen other methods to derive D_MLW. We use IAM-HPLC based D_MLW as a chemical property to distinguish between baseline narcosis and specific MoA, for reported acute toxicity endpoints on two separate sets of chemicals. The first set comprised 94 chemicals of US EPA's acute fish toxicity database: 47 categorized as narcosis MoA, 27 with specific MoA, and 20 predominantly ionic chemicals with mostly unknown MoA. The narcosis MoA chemicals clustered around the median narcosis critical membrane burden (CMB_narc) of 140 mmol kg^-1 lipid, with a lower limit of 14 mmol kg^-1 lipid, including all chemicals labelled Narcosis_I and Narcosis_II. This maximum 'toxic ratio' (TR) between CMB_narc and the lower limit narcosis endpoint is thus 10. For 23/28 specific MoA chemicals a TR >10 was derived, indicative of a specific adverse effect pathway related to acute toxicity. For 10/12 cations categorized as "unsure amines", the TR <10 suggests that these affect fish via narcosis MoA. The second set comprised 29 herbicides, including 17 dissociated acids, and evaluated the TR for acute toxic effect concentrations to likely sensitive aquatic plant species (green algae and macrophytes Lemna and Myriophyllum), and non-target animal species (invertebrates and fish). For 21/29 herbicides, a TR >10 indicated a specific toxic mode of action other than narcosis for at least one of these aquatic primary producers. Fish and invertebrate TRs were mostly <10, particularly for neutral herbicides, but for acidic herbicides a TR >10 indicated specific adverse effects in non-target animals. The established critical membrane approach to derive the TR provides for useful contribution to the weight of evidence to bin a chemical as having a narcosis MoA or less likely to have acute toxicity caused by a more specific adverse effect pathway. After proper calibration, the chromatographic assay provides consistent and efficient experimental input for both neutral and ionizable chemicals to this approach.

Recent advances for estimating environmental properties for small molecules from chromatographic measurements and the solvation parameter model

The transfer of neutral compounds between immiscible phases in chromatographic or environmental systems can be described by six solute properties (solute descriptors) using the solvation parameter model. The solute descriptors are size (McGowan's characteristic volume), V, excess molar refraction, E, dipolarity/polarizability, S, hydrogen-bond acidity and basicity, A and B, and the gas-liquid partition constant on n-hexadecane at 298.15 K, L. V and E for liquids are accessible by calculation but the other descriptors and E for solids are determined experimentally by chromatographic, liquid-liquid partition, and solubility measurements. These solute descriptors are available for several thousand compounds in the Abraham solute descriptor databases and for several hundred compounds in the WSU experimental solute descriptor database. In the first part of this review, we highlight features important in defining each descriptor, their experimental determination, compare descriptor quality for the two organized descriptor databases, and methods for estimating Abraham solute descriptors. In the second part we focus on recent applications of the solvation parameter model to characterize environmental systems and its use for the identification of surrogate chromatographic models for estimating environmental properties.

A joint optimization QSAR model of fathead minnow acute toxicity based on a radial basis function neural network and its consensus modeling

Acute toxicity of the fathead minnow (Pimephales promelas) is an important indicator to evaluate the hazards and risks of compounds in aquatic environments. The aim of our study is to explore the predictive power of the quantitative structure-activity relationship (QSAR) model based on a radial basis function (RBF) neural network with the joint optimization method to study the acute toxicity mechanism, and to develop a potential acute toxicity prediction model, for fathead minnow. To ensure the symmetry and fairness of the data splitting and to generate multiple chemically diverse training and validation sets, we used a self-organizing mapping (SOM) neural network to split the modeling dataset (containing 955 compounds) characterized by PaDEL-descriptors. After preliminary selection of descriptors via the mean decrease impurity method, a hybrid quantum particle swarm optimization (HQPSO) algorithm was used to jointly optimize the parameters of RBF and select the key descriptors. We established 20 RBF-based QSAR models, and the statistical results showed that the 10-fold cross-validation results (R cv10 2) and the adjusted coefficients of determination (R adj 2) were all great than 0.7 and 0.8, respectively. The Q ext 2 of these models was between 0.6480 and 0.7317, and the R ext 2 was between 0.6563 and 0.7318. Combined with the frequency and importance of the descriptors used in RBF-based models, and the correlation between the descriptors and acute toxicity, we concluded that the water distribution coefficient, molar refractivity, and first ionization potential are important factors affecting the acute toxicity of fathead minnow. A consensus QSAR model with RBF-based models was established; this model showed good performance with R 2 = 0.9118, R cv10 2 = 0.7632, and Q ext 2 = 0.7430. A frequency weighted and distance (FWD)-based application domain (AD) definition method was proposed, and the outliers were analyzed carefully. Compared with previous studies the method proposed in this paper has obvious advantages and its robustness and external predictive power are also better than Xgboost-based model. It is an effective QSAR modeling method.

Estimation of the octanol-water distribution coefficient of basic compounds by a cationic microemulsion electrokinetic chromatography system

The octanol-water partition coefficient (P_o/w), or the octanol-water distribution coefficient (D_o/w) for ionized compounds, is a key parameter in the drug development process. In a previous work, this parameter was estimated through the retention factor measurements in a sodium dodecyl sulfate (SDS) - microemulsion electrokinetic chromatography (MEEKC) system for acidic compounds. Nonetheless, when ionized basic compounds were analyzed, undesirable ion pairs were formed with the anionic surfactant and avoided a good estimation of log D_o/w. For this reason, an alternative MEEKC system based on a cationic surfactant has been evaluated to estimate P_o/w or D_o/w of neutral compounds and ionized bases. To this end, it has been characterized through the solvation parameter model (SPM) and compared to the octanol-water partition system. Results pointed out that both systems show a similar partition behavior. Hence, the log P_o/w of a set of neutral compounds has been successfully correlated against the logarithm of the retention factor (log k) determined in this MEEKC system. Then, the log D_o/w of 6 model bases have been estimated at different pH values and they have been compared to data from the literature, determined by the reference shake-flask and potentiometric methods. Good agreement has been observed between the literature and the estimated values when the base is neutral or partially ionized (up to 99% of ionization).

Synthesis of 1-Substituted 1H-1,2,3,4-Tetrazoles Using Biosynthesized Ag/Sodium Borosilicate Nanocomposite

An expedient solvent-free methodology has been developed to produce 1-substituted 1H-1,2,3,4-tetrazoles using sodium borosilicate glass-supported silver nanoparticles (Ag NPs) as a novel heterogeneous catalyst. A cost-efficient, facile, and greener method was deployed for the creation of Ag/sodium borosilicate nanocomposite (ASBN) catalyst by using Aleurites moluccana leaf extract as a stabilizing and reducing agent. The ASBN catalyst was identified using the latest microscopic and spectroscopic techniques such as FT-IR, TEM, FESEM, XRD, EDS, and elemental mapping. The deployment of this new catalyst enables the preparation of assorted 1-substituted tetrazoles in good to high yields via an easy work-up procedure in a relatively short reaction time under environmentally friendly conditions without using harmful and toxic reducing agents. The ASBN catalyst can be recycled and reused multiple times without meaningful loss of activity. To extend the application of the ASBN, the performance of the quantitative structure-activity relationships model was investigated for protein binding and toxicity hazard considerations.