QSAR analyses of the toxicity of aliphatic carboxylic acids and salts to Tetrahymena pyriformis

Recommendations for Improving Methods and Models for Aquatic Hazard Assessment of Ionizable Organic Chemicals

Ionizable organic chemicals (IOCs) such as organic acids and bases are an important substance class requiring aquatic hazard evaluation. Although the aquatic toxicity of IOCs is highly dependent on the water pH, many toxicity studies in the literature cannot be interpreted because pH was not reported or not kept constant during the experiment, calling for an adaptation and improvement of testing guidelines. The modulating influence of pH on toxicity is mainly caused by pH-dependent uptake and bioaccumulation of IOCs, which can be described by ion-trapping and toxicokinetic models. The internal effect concentrations of IOCs were found to be independent of the external pH because of organisms' and cells' ability to maintain a stable internal pH milieu. If the external pH is close to the internal pH, existing quantitative structure-activity relationships (QSARs) for neutral organics can be adapted by substituting the octanol-water partition coefficient by the ionization-corrected liposome-water distribution ratio as the hydrophobicity descriptor, demonstrated by modification of the target lipid model. Charged, zwitterionic and neutral species of an IOC can all contribute to observed toxicity, either through concentration-additive mixture effects or by interaction of different species, as is the case for uncoupling of mitochondrial respiration. For specifically acting IOCs, we recommend a 2-step screening procedure with ion-trapping/QSAR models used to predict the baseline toxicity, followed by adjustment using the toxic ratio derived from in vitro systems. Receptor- or plasma-binding models also show promise for elucidating IOC toxicity. The present review is intended to help demystify the ecotoxicity of IOCs and provide recommendations for their hazard and risk assessment. Environ Toxicol Chem 2020;39:269-286. © 2019 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals, Inc. on behalf of SETAC.

Comparative Quantitative Structure–Activity–Activity Relationships for Toxicity to Tetrahymena pyriformis and Pimephales promelas

An approach for predicting acute aquatic toxicity, in the form of a quantitative structure–activity–activity relationship (QSAAR), is described. This study assessed relative toxic effects to a fish, Pimephales promelas, and a ciliate, Tetrahymena pyriformis, and attempted to form relationships between them. A good agreement between toxic potencies ( R2 = 0.754) was found for a chemically diverse dataset of 364 compounds, when using toxicity to the ciliate as a surrogate to that for fish. This relationship was extended by adding three theoretical structural descriptors of the molecules. The inclusion of these descriptors improved the relationship further ( R2 = 0.824). The structural features that were found to improve the extrapolation between the toxicity to the two different species were related to the electron distribution of the carbon skeleton of the toxicant, its hydrogen-bonding ability, and its relative nitrogen content. Such a QSAAR approach provides a potential tool for predicting the toxicities of chemicals for environmental risk assessment and thus for reducing animal tests.

Impact of temperature, pH, and salinity changes on the physico-chemical properties of model naphthenic acids

The effects of temperature, pH, and salinity change on naphthenic acids (NAs) present in oil-sands process wastewater were modeled for 55 representative NAs. COSMO-RS was used to estimate octanol-water (KOW) and octanol-air (KOA) partition ratios and Henry's law constants (H). Validation with experimental carboxylic acid data yielded log KOW and log H RMS errors of 0.45 and 0.55 respectively. Calculations of log KOW, (or log D, for pH-dependence), log KOA and log H (or log HD, for pH-dependence) were made for model NAs between -20 °C and 40 °C, pH between 0 and 14, and salinity between 0 and 3 g NaCl L(-1). Temperature increase by 60 °C resulted in 3-5 log unit increase in H and a similar magnitude decrease in KOA. pH increase above the NA pKa resulted in a dramatic decrease in both log D and log HD. Salinity increase over the 0-3 g NaCl L(-1) range resulted in a 0.3 log unit increase on average for KOW and H values. Log KOW values of the sodium salt and anion of the conjugate base were also estimated to examine their potential for contribution to the overall partitioning of NAs. Sodium salts and anions of naphthenic acids are predicted to have on average 4 log units and 6 log units lower log KOW values, respectively, with respect to the corresponding neutral NA. Partitioning properties are profoundly influenced by the by the relative prevailing pH and the substance's pKa at the relevant temperature.

Toxicity and composition profiles of solid phase extracts of oil sands process-affected water

After fractionation using sequential solid phase extraction, the presence of toxic components in oil sands process-affected water (OSPW) was detected by the Microtox® acute toxicity assay using effect-directed analysis. The composition of each fraction was determined by high-resolution electrospray ionization-Orbitrap mass spectrometry. Partial least-squares discriminant analysis (PLS-DA) was used to determine which chemical constituents in all seven fractions co-varied most strongly with toxicity. Although O2 compounds with double bond equivalence (DBE) between 3 and 9 positively correlated with toxicity, C15-C18 O2-NAs with DBE=4 (tricyclic structure), as well as C14-C17 O2-NAs with DBE=3 (bicyclic structure), were found to be most likely associated with OSPW toxicity, consistent with published toxicity studies of surrogate NAs. O4, many O3 (i.e. possibly hydroxylated O2 c-NAs) and a few O2 compounds were found to negatively correlate with toxicity. The results demonstrate the utility of the fractionation and the PLS-DA approach for evaluating composition-response relationships in a complex mixture and also contribute to a better understanding of the toxic compounds in OSPW. These findings will help to focus study on the most environmentally significant components in OSPW.

Modelling the behaviour of switchable-hydrophilicity solvents

Mathematical modelling of the chemical and phase equilibria of switchable-hydrophilicity solvents helps explain and predict their properties and performance.

The effect of oil sands process-affected water and naphthenic acids on the germination and development of Arabidopsis

Oil sands mining in the Athabasca region of northern Alberta results in the production of large volumes of oil sands process-affected water (OSPW). We have evaluated the effects of OSPW, the acid extractable organic (AEO) fraction of OSPW, and individual naphthenic acids (NAs) on the germination and development of the model plant, Arabidopsis thaliana (Arabidopsis). The surrogate NAs that were selected for this study were petroleum NAs that have been used in previous toxicology studies and may not represent OSPW NAs. A tricyclic diamondoid NA that was recently identified as a component of OSPW served as a model NA in this study. Germination of Arabidopsis seeds was not inhibited when grown on medium containing up to 75% OSPW or by 50mgL(-1) AEO. However, simultaneous exposure to three simple, single-ringed surrogate NAs or a double-ringed surrogate NA had an inhibitory effect on germination at a concentration of 10mgL(-1), whereas inhibition of germination by the diamondoid model NA was observed only at 50mgL(-1). Seedling root growth was impaired by treatment with low concentrations of OSPW, and exposure to higher concentrations of OSPW resulted in increased growth inhibition of roots and primary leaves, and caused bleaching of cotyledons. Treatment with single- or double-ringed surrogate NAs at 10mgL(-1) severely impaired seedling growth. AEO or diamondoid NA treatment was less toxic, but resulted in severely impaired growth at 50mgL(-1). At low NA concentrations there was occasionally a stimulatory effect on root and shoot growth, possibly owing to the broad structural similarity of some NAs to known plant growth regulators such as auxins. This report provides a foundation for future studies aimed at using Arabidopsis as a biosensor for toxicity and to identify genes with possible roles in NA phytoremediation.

From data point timelines to a well curated data set, data mining of experimental data and chemical structure data from scientific articles, problems and possible solutions

The scientific literature is important source of experimental and chemical structure data. Very often this data has been harvested into smaller or bigger data collections leaving the data quality and curation issues on shoulders of users. The current research presents a systematic and reproducible workflow for collecting series of data points from scientific literature and assembling a database that is suitable for the purposes of high quality modelling and decision support. The quality assurance aspect of the workflow is concerned with the curation of both chemical structures and associated toxicity values at (1) single data point level and (2) collection of data points level. The assembly of a database employs a novel "timeline" approach. The workflow is implemented as a software solution and its applicability is demonstrated on the example of the Tetrahymena pyriformis acute aquatic toxicity endpoint. A literature collection of 86 primary publications for T. pyriformis was found to contain 2,072 chemical compounds and 2,498 unique toxicity values, which divide into 2,440 numerical and 58 textual values. Every chemical compound was assigned to a preferred toxicity value. Examples for most common chemical and toxicological data curation scenarios are discussed.

Degradation and aquatic toxicity of naphthenic acids in oil sands process-affected waters using simulated wetlands

Oil sands process-affected waters (OSPWs) produced during the extraction of bitumen at the Athabasca Oil Sands (AOS) located in northeastern Alberta, Canada, are toxic to many aquatic organisms. Much of this toxicity is related to a group of dissolved organic acids known as naphthenic acids (NAs). Naphthenic acids are a natural component of bitumen and are released into process water during the separation of bitumen from the oil sand ore by a caustic hot water extraction process. Using laboratory microcosms as an analogue of a proposed constructed wetland reclamation strategy for OSPW, we evaluated the effectiveness of these microcosms in degrading NAs and reducing the aquatic toxicity of OSPW over a 52-week test period. Experimental manipulations included two sources of OSPW (one from Syncrude Canada Ltd. and one from Suncor Energy Inc.), two different hydraulic retention times (HRTs; 40 and 400 d), and increased nutrient availability (added nitrate and phosphate). Microcosms with a longer HRT (for both OSPWs) showed higher reductions in total NAs concentrations (64-74% NAs reduction, p<0.05) over the test period, while nutrient enrichment appeared to have little effect. A 96 h static acute rainbow trout (Oncorhynchus mykiss) bioassay showed that the initial acute toxicity of Syncrude OSPW (LC50=67% v/v) was reduced (LC50>100% v/v) independent of HRT. However, EC20s from separate Microtox® bioassays were relatively unchanged when comparing the input and microcosm waters at both HRTs over the 52-week study period (p>0.05), indicating that some sub-lethal toxicity persisted under these experimental conditions. The present study demonstrated that given sufficiently long HRTs, simulated wetland microcosms containing OSPW significantly reduced total NAs concentrations and acute toxicity, but left behind a persistent component of the NAs mixture that appeared to be associated with residual chronic toxicity.

Predicted toxicity of naphthenic acids present in oil sands process-affected waters to a range of environmental and human endpoints

Naphthenic acids (NAs) are considered to be a major toxic component of oil sands process-affected waters (OSPW) and are also widely used for industrial processes. The effects of previously identified NAs (54 in total), together with six alkylphenols, were modelled for a range of environmental and human toxicity related endpoints using ADMET predictor™ software. In addition to the models, experimental CALUX® assays were performed on seven tricyclic diamondoid acids. Most of the NAs modelled were predicted to have lethal median concentrations (LC(50)) >100 μM for the three aquatic species modelled. Polycyclic acids containing a single aromatic ring were predicted to be the most toxic to fathead minnows with LC(50)s typically ca 1 μM. Some of these compounds were also predicted to be the most carcinogenic (based on rat and mouse models), possess human estrogenic and androgenic activity and potentially disrupt reproductive processes. Some aliphatic pentacyclic acids also were predicted to exhibit androgenic activity and, uniquely amongst the compounds tested, act as substrates for the cytochrome P450 enzyme CYP3A4. Consistent with the models' predictions for the tricyclic acids, no estrogenic or androgenic activity was detected by ER/AR CALUX®. Further experimental validation of the predictions should now be performed for the compounds highlighted by the models (e.g. priority should perhaps be focused on the polycyclic monoaromatic acids and the aliphatic pentacyclic acids). If shown to be accurate, these compounds can then be targeted for toxicity reduction remediation efforts.

Naphthenic acids degradation and toxicity mitigation in tailings wastewater systems and aquatic environments: a review

Naphthenic acids, NAs (classical formula C(n)H(2n+z)O(2), where n is the carbon numbers, z represents zero or negative even integers), found in oil sands process waters (OSPWs), are toxic to aquatic environments depending upon several factors such as pH, salinity, molecular size and chemical structure of NAs. Among various available methods, biodegradation seems to be generally the most cost-effective method for decreasing concentrations of NAs (n ≤ 21) and reducing their associated toxicity in OSPW, however the mechanism by which the biodegradation of NAs occurs are poorly understood. Ozonation is superior over biodegradation in decreasing higher molecular weight alkyl branched NAs (preferentially, n ≥ 22, -6 ≥ z ≥ -12) as well as enabling accelerated biodegradation and reducing toxicity. Photolysis (UV at 254 nm) is effective in cleaving higher molecular weight NAs into smaller fragments that will be easier for microorganisms to degrade, whereas photocatalysis can metabolize selective NAs (0 ≥ z ≥ -6) efficiently and minimize their associated toxicity. Phytoremediation is applicable for metabolizing specific NAs (O(2), O(3), O(4), and O(5) species) and minimizing their associated toxicities. Petroleum coke (PC) adsorption is effective in reducing the more structurally complex NAs (preferentially 12 ≥ n ≥ 18 and z = -10, -12) and their toxicity in OSPWs, depending upon the PC content, pH and temperature. Several factors have influence on the degradation of NAs in OSPWs and aquatic environments, which include molecular mass and chemical structure of NAs, sediment structure, temperature, pH, dissolved oxygen, nutrients, and bacteria types.

Repeatability analysis of the Tetrahymena pyriformis population growth impairment assay

Assessments necessary to ensure the safety of both humans and the environment are challenged by the sheer number of chemicals in use today. Chemical legislation, such as REACH, aims to use alternative methods to reduce the reliance on in vivo animal testing. Consequently, databases such as the TETRATOX database, containing data from the Tetrahymena pyriformis population growth impairment assay, have been used extensively to develop computational models which aid in priority setting and initial hazard assessments. To use any toxicological data, an assessment of quality is required. One important aspect of quality is the repeatability of the assay. This study considered TETRATOX assay data for 85 structurally and mechanistically diverse compounds. The repeatability of replicate determinations was assessed and factors relating to repeatability are discussed. Despite the majority of compounds demonstrating excellent repeatability, it was found that the mechanism of action is likely to be a modulating factor, with compounds acting via electrophilic mechanisms being more likely to exhibit reduced repeatability than those acting via narcotic mechanisms. It is evident from this study that the TETRATOX assay is a robust and highly repeatable assay, suitable for use in toxicological modelling studies and priority setting.

Biological Read-Across: Mechanistically-Based Species–Species and Endpoint–Endpoint Extrapolations

This chapter describes the development and use of relationships between toxicity data, for the same chemicals, derived from different species. These relationships provide simple models to predict toxicity of, potentially, a higher species from that of a lower species. Approaches to the formation of these models are described, notably inter-species relationships, quantitative structure-activity-activity relationships and prediction models. It is noted that the best extrapolations are for closely related species i.e. within taxa. In addition, forming groups or categories of compounds according to common mechanisms of toxic action improves the correlation for extrapolations from lower to higher species. A freely available software package, Web-ICE, is introduced as a suitable tool to apply these methods.

Use of a (quantitative) structure-activity relationship [(Q)SAR] model to predict the toxicity of naphthenic acids

Naphthenic acids (NA) are a complex mixture of carboxylic acids that are natural constituents of oil sand found in north-eastern Alberta, Canada. NA are released and concentrated in the alkaline water used in the extraction of bitumen from oil sand sediment. NA have been identified as the principal toxic components of oil sands process-affected water (OSPW), and microbial degradation of lower molecular weight (MW) NA decreases the toxicity of NA mixtures in OSPW. Analysis by proton nuclear magnetic resonance spectroscopy indicated that larger, more cyclic NA contain greater carboxylic acid content, thereby decreasing their hydrophobicity and acute toxicity in comparison to lower MW NA. The relationship between the acute toxicity of NA and hydrophobicity suggests that narcosis is the probable mode of acute toxic action. The applicability of a (quantitative) structure-activity relationship [(Q)SAR] model to accurately predict the toxicity of NA-like surrogates was investigated. The U.S. Environmental Protection Agency (EPA) ECOSAR model predicted the toxicity of NA-like surrogates with acceptable accuracy in comparison to observed toxicity values from Vibrio fischeri and Daphnia magna assays, indicating that the model has potential to serve as a prioritization tool for identifying NA structures likely to produce an increased toxicity. Investigating NA of equal MW, the ECOSAR model predicted increased toxic potency for NA containing fewer carbon rings. Furthermore, NA structures with a linear grouping of carbon rings had a greater predicted toxic potency than structures containing carbon rings in a clustered grouping.

Aqueous in situ derivatization of carboxylic acids by an ionic carbodiimide and 2,2,2-trifluoroethylamine for electron-capture detection

We report a technique for the rapid, room temperature derivatization of aqueous carboxylic acids to the corresponding 2,2,2-trifluoroethylamide derivative. 3-Ethyl-1-[3-(dimethylamino)propyl]carbodiimide hydrochloride (EDC) and 2,2,2-trifluoroethylamine hydrochloride (TFEA) were added to aqueous samples of several acids of interest in environmental analytical chemistry, including benzoic acid, ibuprofen, clofibric acid, monochloroacetic acid, dicholoroacetic acid, bromoacetic acid, monochlorophenoxy acetic acid, dichlorophenoxy acetic acid, trichlorophenoxy acetic acid, acetylsalicylic acid, and cholorobenzoic acid. Amidization was essentially complete within ten minutes, and subsequent liquid-liquid extraction of the amides with methyl tert-butyl ether (MTBE) demonstrated recoveries of over 85%. The starting materials, both quaternary ammonium salts, were not co-extracted with the derivative, yielding much cleaner samples than historically obtained from carbodiimide based techniques. The fluorinated amides produced had excellent chromatographic characteristics for gas chromatography and were easily detected by electron-capture detection (ECD) or electron impact mass spectrometry. This method is suggested as a sensitive alternative to more traditional acidification, extraction, and ex situ derivatization techniques.

Acute toxicity of fatty acids to the freshwater green alga Selenastrum capricornutum

The acute toxicity of fatty acids (C14 to C18) commonly found in wood was determined by the standard algal growth inhibition test using the freshwater green alga Selenastrum capricornutum. Toxicity, quantified as IC50 values, varied depending on the number of total carbons and double bonds. Of the tested acids, oleic (cis-9-octadecenoic) acid showed the highest toxicity (72-h IC50 = 0.47 mg/L) to the alga, and triolein, a triglyceride of oleic acid, showed no apparent toxicity. Further examination of a series of C18:1 acids with a double bond at the 6, 11, or 12 position revealed that both double-bond position and cis or trans configuration affected toxicity. The 72-h IC50 data for these fatty acids and related compounds seemed to correlate well with the melting point (mp), showing two separate linear relationships: at mp < 35 degrees C toxicity increased with increasing melting point, and at mp > 40 degrees C toxicity decreased with melting point.

Computer-modeling-based QSARs for analyzing experimental data on biotransformation and toxicity

Over the past decades the description of quantitative structure-activity relationships (QSARs) has been undertaken in order to find predictive models and/or mechanistic explanations for chemical as well as biological activities. This includes QSAR studies in toxicology. In an approach beyond the classical QSAR approaches, attempts have been made to define parameters for the QSAR studies on the basis of quantum mechanical computer calculations. The conversion of relatively small xenobiotics within the active sites of biotransformation enzymes can be expected to follow the general rules of chemistry. This makes the description of QSARs on the basis of only one parameter, chosen on the basis of insight in the mechanism, feasible. In contrast, toxicological endpoints can very often be the result of more than one physico-chemical interaction of the compound with the model system of interest. Therefore the description of quantitative structure-toxicity relationships often does not follow a one-descriptor mechanistic approach but starts from the other end, describing QSARs by multi-parameter approaches. The present paper focuses on the possibilities and restrictions of using computer-based QSAR modeling for analyzing experimental toxicological data, with emphasis on examples from the field of biotransformation and toxicity.

Reproducibility of toxicity across mode of toxic action in the Tetrahymena population growth impairment assay

Toxicity data collected in a laboratory setting are the primary source of potency information used for regulatory, modeling, or risk assessment purposes. However, the relative reproducibility of such toxicity data is rarely discussed. This study investigated the reproducibility of growth impairment data for the freshwater ciliate Tetrahymena pyriformis exposed to a structurally diverse group of chemicals of varying hydrophobicity within different modes of toxic action, either non-covalent narcosis or covalent electro(nucleo)philicity. The proportions of chemicals representing each mode of toxic action, or mechanism of action within each mode, were not chosen to emulate the occurrence of manufactured chemicals or chemicals within the TETRATOX database. Chemicals for which prior toxicity data existed were re-tested and reproducibility was evaluated. The toxic potency values of the selected chemicals were largely reproducible after re-testing of the toxic potency, as 98% of the chemicals had re-test toxicity values within one log unit of the original potency value. To further scrutinize the reproducibility of toxicity values, differences between values were investigated by mode of toxic action. A stringent criterion for reproducibility was enforced, which dictated that the re-tested toxicity value must be encompassed by the fiducial interval (FI) of the original toxicity value and vice versa for the chemical to be considered reproducible. Toxicity values of 28 of the 50 re-tested chemicals conformed to the criterion set for reproducible values. Of the nonreproducible chemicals, seven were narcotics: four nonpolar or neutral narcotics and three other narcotics (e.g. polar narcotics). However, four of these seven narcotics did have toxicity values encompassed by one FI, but not the other FI. The remaining chemicals that did not have reproducible potency measurements were electro(nucleo)philic in nature. Certain toxicophores were highly represented among these chemicals. These included quinone derivatives, electron releasing amino and hydroxyl moieties, and electron withdrawing nitro substituents, often in tandem with strong leaving groups (i.e. halogens), and unsaturated alcohols. Lack of reproducibility was common among the chemicals that elicited toxicity after either abiotic or biotic transformation. There was no clear trend between hydrophobicity and lack of reproducibility. While data are limited, these results suggest that toxic potency values of chemicals acting via the electro(nucleo)philic mode of toxic action could be more susceptible to non-reproducibility. Ramifications of such lack of reproducibility could manifest in predictive toxicology models and their use in regulatory and risk assessment endeavors.