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

A Review of Mechanistic Models for Predicting Adverse Effects in Sediment Toxicity Testing

Since recognizing the importance of bioavailability for understanding the toxicity of chemicals in sediments, mechanistic modeling has advanced over the last 40 years by building better tools for estimating exposure and making predictions of probable adverse effects. Our review provides an up-to-date survey of the status of mechanistic modeling in contaminated sediment toxicity assessments. Relative to exposure, advances have been most substantial for non-ionic organic contaminants (NOCs) and divalent cationic metals, with several equilibrium partitioning-based (Eq-P) models having been developed. This has included the use of Abraham equations to estimate partition coefficients for environmental media. As a result of the complexity of their partitioning behavior, progress has been less substantial for ionic/polar organic contaminants. When the EqP-based estimates of exposure and bioavailability are combined with water-only effects measurements, predictions of sediment toxicity can be successfully made for NOCs and selected metals. Both species sensitivity distributions and toxicokinetic and toxicodynamic models are increasingly being applied to better predict contaminated sediment toxicity. Furthermore, for some classes of contaminants, such as polycyclic aromatic hydrocarbons, adverse effects can be modeled as mixtures, making the models useful in real-world applications, where contaminants seldomly occur individually. Despite the impressive advances in the development and application of mechanistic models to predict sediment toxicity, several critical research needs remain to be addressed. These needs and others represent the next frontier in the continuing development and application of mechanistic models for informing environmental scientists, managers, and decisions makers of the risks associated with contaminated sediments. Environ Toxicol Chem 2024;43:1778-1794. © 2023 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

In Silico Acute Aquatic Hazard Assessment and Prioritization Using a Grouped Target Site Model: A Case Study of Organic Substances Reported in Permian Basin Hydraulic Fracturing Operations

Hydraulic fracturing (HF) is commonly used to enhance onshore recovery of oil and gas during production. This process involves the use of a variety of chemicals to support the physical extraction of oil and gas, maintain appropriate conditions downhole (e.g., redox conditions, pH), and limit microbial growth. The diversity of chemicals used in HF presents a significant challenge for risk assessment. The objective of the present study is to establish a transparent, reproducible procedure for estimating 5th percentile acute aquatic hazard concentrations (e.g., acute hazard concentration 5th percentiles [HC5s]) for these substances and validating against existing toxicity data. A simplified, grouped target site model (gTSM) was developed using a database (n = 1696) of diverse compounds with known mode of action (MoA) information. Statistical significance testing was employed to reduce model complexity by combining 11 discrete MoAs into three general hazard groups. The new model was trained and validated using an 80:20 allocation of the experimental database. The gTSM predicts toxicity using a combination of target site water partition coefficients and hazard group-based critical target site concentrations. Model performance was comparable to the original TSM using 40% fewer parameters. Model predictions were judged to be sufficiently reliable and the gTSM was further used to prioritize a subset of reported Permian Basin HF substances for risk evaluation. The gTSM was applied to predict hazard groups, species acute toxicity, and acute HC5s for 186 organic compounds (neutral and ionic). Toxicity predictions and acute HC5 estimates were validated against measured acute toxicity data compiled for HF substances. This case study supports the gTSM as an efficient, cost-effective computational tool for rapid aquatic hazard assessment of diverse organic chemicals. Environ Toxicol Chem 2024;00:1-12. © 2024 ExxonMobil Petroleum and Chemical BV. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

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.

Predicting Solute Descriptors for Organic Chemicals by a Deep Neural Network (DNN) Using Basic Chemical Structures and a Surrogate Metric

Solute descriptors have been widely used to model chemical transfer processes through poly-parameter linear free energy relationships (pp-LFERs); however, there are still substantial difficulties in obtaining these descriptors accurately and quickly for new organic chemicals. In this research, models (PaDEL-DNN) that require only SMILES of chemicals were built to satisfactorily estimate pp-LFER descriptors using deep neural networks (DNN) and the PaDEL chemical representation. The PaDEL-DNN-estimated pp-LFER descriptors demonstrated good performance in modeling storage-lipid/water partitioning coefficient (log K_{storage-lipid/water}), bioconcentration factor (BCF), aqueous solubility (ESOL), and hydration free energy (freesolve). Then, assuming that the accuracy in the estimated values of widely available properties, e.g., logP (octanol-water partition coefficient), can calibrate estimates for less available but related properties, we proposed logP as a surrogate metric for evaluating the overall accuracy of the estimated pp-LFER descriptors. When using the pp-LFER descriptors to model log K_{storage-lipid/water}, BCF, ESOL, and freesolve, we achieved around 0.1 log unit lower errors for chemicals whose estimated pp-LFER descriptors were deemed "accurate" by the surrogate metric. The interpretation of the PaDEL-DNN models revealed that, for a given test chemical, having several (around 5) "similar" chemicals in the training data set was crucial for accurate estimation while the remaining less similar training chemicals provided reasonable baseline estimates. Lastly, pp-LFER descriptors for over 2800 persistent, bioaccumulative, and toxic chemicals were reasonably estimated by combining PaDEL-DNN with the surrogate metric. Overall, the PaDEL-DNN/surrogate metric and newly estimated descriptors will greatly benefit chemical transfer modeling.

Exploring QSPR models for predicting PUF-air partition coefficients of organic compounds with linear and nonlinear approaches

Partition coefficients are important parameters for measuring the concentration of chemicals by passive sampling devices. Considering the wide application of the polyurethane foam (PUF) in passive air sampling, an attempt for developing several quantitative structure-property relationship (QSPR) models was made in this work, to predict PUF-air partition coefficients (K_PUF-air) using linear (multiple linear regression, MLR) and non-linear (artificial neural network, ANN and support vector machine, SVM) methods by machine learning. All of the developed models were performed on a dataset of 170 compounds comprising 9 distinct classes. A series of statistical parameters and validation results showed that models had good prediction ability, robustness and goodness-of-fit. Furthermore, the underlying mechanisms of molecular descriptors emphasized that ionization potential, molecular bond, hydrophilicity, size of molecule and valence electron number had dominating influence on the adsorption process of chemicals. Overall, the obtained models were all established on the extensive applicability domains, and thus can be used as effective tools to predict the K_PUF-air of new organic compounds or those have not been synthesized yet which, in turn, could help researchers better understand the mechanistic basis of adsorption behavior of PUF.

Partition of Neutral Molecules and Ions from Water to o-Nitrophenyl Octyl Ether and of Neutral Molecules from the Gas Phase to o-Nitrophenyl Octyl Ether

We have set out an equation for partition of 87 neutral molecules from water to o-nitrophenyl octyl ether, NPOE, an equation for partition of the 87 neutral molecules and 21 ionic species from water to NPOE, and an equation for partition of 87 neutral molecules from the gas phase to NPOE. Comparison with equations for partition into other solvents shows that, as regards partition of neutral (nonelectrolyte) compounds, NPOE would be a good model for 1,2-dichloroethane and for nitrobenzene. In terms of partition of ions and ionic species, NPOE is quite similar to 1,2-dichloroethane and not far away from other aprotic solvents such as nitrobenzene.

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.