Assessing the toxicity of ionic liquids - Application of the critical membrane concentration approach

Assessing Aquatic Baseline Toxicity of Plastic-Associated Chemicals: Development and Validation of the Target Plastic Model

We developed a Target Plastic Model (TPM) to estimate the critical plastic burden of organic toxicants in five types of plastics, namely, polydimethylsiloxane (PDMS), polyoxymethylene (POM), polyacrylate (PA), low-density polyethylene (LDPE), and polyurethane ester (PU), following the Target Lipid Model (TLM) framework. By substituting the lipid-water partition coefficient in the TLM with plastic-water partition coefficients to create TPM, we demonstrated that the biomimetic nature of these plastic phases allows for the calculation of critical plastic burdens of toxicants, similar to the notion of critical lipid burdens in TLM. Following this approach, the critical plastic burdens of baseline (n = 115), less-inert (n = 73), and reactive (n = 75) toxicants ranged from 0.17 to 51.33, 0.04 to 26.62, and 1.00 × 10-6 to 6.78 × 10-4 mmol/kg of plastic, respectively. Our study showed that PDMS, PA, POM, PE, and PU are similar to biomembranes in mimicking the passive exchange of chemicals with the water phase. Using the TPM, median lethal concentration (LC50) values for fish exposed to baseline toxicants were predicted, and the results agreed with experimental values, with RMSE ranging from 0.311 to 0.538 log unit. Similarly, for the same data set of baseline toxicants, other widely used models, including the TLM (RMSE: 0.32-0.34), ECOSAR (RMSE: 0.35), and the Abraham Solvation Model (ASM; RMSE: 0.31), demonstrated comparable agreement between experimental and predicted values. For less inert chemicals, predictions were within a factor of 5 of experimental values. Comparatively, ASM and ECOSAR showed predictions within a factor of 2 and 3, respectively. The TLM based on phospholipid had predictions within a factor of 3 and octanol within a factor of 4, indicating that the TPM's performance for less inert chemicals is comparable to these established models. Unlike these methods, the TPM requires only the knowledge of plastic bound concentration for a given plastic phase to calculate baseline toxic units, bypassing the need for extensive LC50 and plastic-water partition coefficient data, which are often limited for emerging chemicals. Taken together, the TPM can provide valuable insights into the toxicities of chemicals associated with environmental plastic phases, assisting in selecting the best polymeric phase for passive sampling and designing better passive dosing techniques for toxicity experiments.

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.

Screening the baseline fish bioconcentration factor of various types of surfactants using phospholipid binding data

Fish bioconcentration factors (BCFs) are commonly used in chemical hazard and risk assessment. For neutral organic chemicals BCFs are positively correlated with the octanol-water partition ratio (K_OW), but K_OW is not a reliable parameter for surfactants. Membrane lipid-water distribution ratios (D_MLW) can be accurately measured for all kinds of surfactants, using phospholipid-based sorbents. This study first demonstrates that D_MLW values for ionic surfactants are more than 100 000 times higher than the partition ratio to fish-oil, representing neutral storage lipid. A non-ionic alcohol ethoxylate surfactant showed almost equal affinity for both lipid types. Accordingly, a baseline screening BCF value for surfactants (BCF_baseline) can be approximated for ionic surfactants by multiplying D_MLW by the phospholipid fraction in tissue, and for non-ionic surfactants by multiplying D_MLW by the total lipid fraction. We measured D_MLW values for surfactant structures, including linear and branched alkylbenzenesulfonates, an alkylsulfoacetate and an alkylethersulfate, bis(2-ethylhexyl)-surfactants (e.g., docusate), zwitterionic alkylbetaines and alkylamine-oxides, and a polyprotic diamine. Together with sixty previously published D_MLW values for surfactants, structure-activity relationships were derived to elucidate the influence of surfactant specific molecular features on D_MLW. For 23 surfactant types, we established the alkyl chain length at which BCF_baseline would exceed the EU REACH bioaccumulation (B) threshold of 2000 L kg^-1, and would therefore require higher tier assessments to further refine the BCF estimate. Finally, the derived BCF_baseline are compared with measured literature in vivo BCF data where available, suggesting that refinements, most notably reliable estimates of biotransformation rates, are needed for most surfactant types.

Biotransformation Potential of Cationic Surfactants in Fish Assessed with Rainbow Trout Liver S9 Fractions

Biotransformation may substantially reduce the extent to which organic environmental contaminants accumulate in fish. Presently, however, relatively little is known regarding the biotransformation of ionized chemicals, including cationic surfactants, in aquatic organisms. To address this deficiency, a rainbow trout liver S9 substrate depletion assay (RT-S9) was used to measure in vitro intrinsic clearance rates (CL_int ; ml min^-1 g liver^-1 ) for 22 cationic surfactants that differ with respect to alkyl chain length and degree of methylation on the charged nitrogen atom. None of the quaternary N,N,N-trimethylalkylammonium compounds exhibited significant clearance. Rapid clearance was observed for N,N-dimethylalkylamines, and slower rates of clearance were measured for N-methylalkylamine analogs. Clearance rates for primary alkylamines were generally close to or below detectable levels. For the N-methylalkylamines and N,N-dimethylalkylamines, the highest CL_int values were measured for C₁₀ -C₁₂ homologs; substantially lower clearance rates were observed for homologs containing shorter or longer carbon chains. Based on its cofactor dependency, biotransformation of C₁₂ -N,N-dimethylamine appears to involve one or more cytochrome P450-dependent reaction pathways, and sulfonation. On a molar basis, N-demethylation metabolites accounted for up to 25% of the N,N-dimethylalkylamines removed during the 2-h assay, and up to 55% of the removed N-methylalkylamines. These N-demethylation products possess greater metabolic stability in the RT-S9 assay than the parent structures from which they derive and may contribute to the overall risk of ionizable alkylamines. The results of these studies provide a set of consistently determined CL_int values that may be extrapolated to whole trout to inform in silico bioaccumulation assessments. Environ Toxicol Chem 2021;40:3123-3136. © 2021 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Review of the toxic effects of ionic liquids

Interest in ionic liquids (ILs), called green or designer solvents, has been increasing because of their excellent properties such as thermal stability and low vapor pressure; thus, they can replace harmful organic chemicals and help several industrial fields e.g., energy-storage materials production and biomaterial pretreatment. However, the claim that ILs are green solvents should be carefully considered from an environmental perspective. ILs, given their minimal vapor pressure, may not directly cause atmospheric pollution. However, they have the potential to cause adverse effects if leaked into the environment, for instance if they are spilled due to human mistakes or technical errors. To estimate the risks of ILs, numerous ILs have had their toxicity assessed toward several micro- and macro-organisms over the past few decades. Since the toxic effects of ILs depend on the method of estimating toxicity, it is necessary to briefly summarize and comprehensively discuss the biological effects of ILs according to their structure and toxicity testing levels. This can help simplify our understanding of the toxicity of ILs. Therefore, in this review, we discuss the key findings of toxicological information of ILs, collect some toxicity data of ILs to different species, and explain the influence of IL structure on their toxic properties. In the discussion, we estimated two different sensitivity values of toxicity testing levels depending on the experiment condition, which are theoretical magnitudes of the inherent sensitivity of toxicity testing levels in various conditions and their changes in biological response according to the change in IL structure. Finally, some perspectives, future research directions, and limitations to toxicological research of ILs, presented so far, are discussed.

Environmental Sorption Behavior of Ionic and Ionizable Organic Chemicals

Traditionally our tools for environmental risk assessment of organic chemicals have been developed for neutral chemicals. However, many commercial chemicals are ionic or ionizable and require different tools and approaches for their assessment. In recent years this task starts to obtain increasing attention but our understanding for their environmental fate is still far behind that for neutral chemicals. This review first gives an overview on the principles that govern ionic partitioning in environmental systems which are more complex than the simple partition processes of neutral chemicals. Second, a summary of our current knowledge on various topics such as bioaccumulation, sorption in soils, and nonspecific-toxicity reveals that ionic species can actually be quite hydrophobic contrary to commonly held beliefs. Eventually, we discuss existing models for the quantitative prediction of organic ions' sorption in soils and biota. We have to assert that the available model tools are quite restricted in their application range compared to neutral chemicals which is due to the higher complexity of the various ionic sorption processes. In order to further advance our understanding more high-quality sorption data are needed with a focus on multivalent and zwitterionic ions in all partition systems as well as cations in biological matrices.

Benzalkonium ion sorption to peat and clays: Relative contributions of ion exchange and van der Waals interactions

Due to their use in various domestic and industrial formulations, benzalkonium compounds have been isolated in many environmental matrices. Sorption to soil components has been shown to play a key role in their environmental fate. Whereas sorption of benzalkonium compounds to soils is attributed to cation exchange and van der Waals forces, the relative contributions of these two mechanisms at environmental levels have not been clearly defined. In this study a previously reported algal toxicity assay-based method was employed to determine the distribution coefficients (K_d) of benzalkonium compounds between water and soil components, at environmental concentrations. Cation exchange capacity corrected K_d values for organic matter and clays were all within one order of magnitude. This implies that ion exchange is the dominant mechanism of sorption for benzalkonium compounds. When the sorption data were used to compute sorption energies for four homologues of benzalkonium ions, the magnitude of the free energy change of sorption increased with size of the molecule. The increase in sorption energy could be partly explained by increased energy of hydration with addition of methylene groups to the alkyl chain. A model that predicts sorption coefficients of benzalkonium compounds to soils using organic carbon content and cation exchange capacity was also defined. When tested using an artificial soil, the model estimates were all within one order of magnitude of the experimental values.

Ionic Liquid Enriches the Antibiotic Resistome, Especially Efflux Pump Genes, Before Significantly Affecting Microbial Community Structure

An expanding list of chemicals may permeabilize bacterial cells and facilitate horizontal gene transfer (HGT), which enhances propagation of antibiotic resistance genes (ARGs) in the environment. Previous studies showed that 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIm][PF6]), an ionic liquid, can facilitate HGT of some ARGs among bacteria. However, the dynamic response of a wider range of ARGs and associated mobile genetic elements (MGEs) in different environments is unknown. Here, we used metagenomic tools to study shifts of the resistome and microbiome in both sediments and freshwater microcosms exposed to [BMIm][PF6]. Exposure for 16 h to 0.1 or 1.0 g/L significantly enriched more than 207 ARG subtypes primarily encoding efflux pumps in freshwater microcosms as well as cultivable antibiotic-resistant bacteria. This resistome enrichment was attributed to HGT facilitated by MGEs (428 plasmids, 61 integron-integrase genes, and 45 gene cassettes were enriched) as well as to HGT-related functional genes. Interestingly, resistome enrichment occurred fast (within 16 h) after [BMIm][PF6] exposure, before any significant changes in bacterial community structure. Similar ARG enrichment occurred in sediment microcosms exposed to [BMIm][PF6] for 28 d, and this longer exposure affected the microbial community structure (e.g., Proteobacteria abundance increased significantly). Overall, this study suggests that [BMIm][PF6] releases could rapidly enrich the antibiotic resistome in receiving environments by increasing HGT and fortuitously selecting for efflux pump genes, thus contributing to ARG propagation.

A perspective on the role of fugacity and activity for evaluating the PBT properties of organic chemicals and providing a multi-media synoptic indicator of environmental contamination

The hazard and risk from organic chemicals present in the environment are routinely evaluated using P (persistence), B (bioaccumulation) and T (toxicity) criteria. We present a perspective on how models based on the equilibrium criteria of fugacity and chemical activity can contribute to all three evaluations, thus providing a consistent and seamless assessment process. Persistence and its closely related, but fundamentally different property, chemical residence time can be determined from degradation half-lives (typically obtained from monitoring data) and multi-media mass balance models describing chemical fate and transport in real or evaluative environments. Bioaccumulation is driven by equilibrium partitioning processes that can be estimated from fugacity models treating uptake by respiration and diet in single organisms and in food webs, most commonly for aquatic systems for which confirmatory monitoring data can be obtained. Biomagnification is readily evaluated both experimentally and as a fugacity (or activity) increase from prey to predator. Toxicokinetic evaluations of chemical uptake and disposition in a variety of organisms, including humans, can be obtained using fugacity- and physiologically-based pharmacokinetic models. Toxicodynamic processes of chemical interactions with organisms leading to adverse toxic outcomes are less obviously amenable to fugacity assessment with the notable exception of baseline toxicity (narcosis), the most common mode of action of commercial chemicals. It is shown that a range of fugacities can be defined and correlated for specific species thus enabling direct comparison of fugacities estimated or monitored in the environment with fugacities resulting in baseline toxicity. Finally, it is shown that fugacity and activity can serve as multi-media metrics of chemical contamination in ecosystems, thus enabling identification of specific regions and species at greatest risk in a pictorial format that is readily understandable to scientific and regulatory communities and to the general public.

Mode of Action Classifications in the EnviroTox Database: Development and Implementation of a Consensus MOA Classification

Multiple mode of action (MOA) frameworks have been developed in aquatic ecotoxicology, mainly based on fish toxicity. These frameworks provide information on a key determinant of chemical toxicity, but the MOA categories and level of specificity remain unique to each of the classification schemes. The present study aimed to develop a consensus MOA assignment within EnviroTox, a curated in vivo aquatic toxicity database, based on the following MOA classification schemes: Verhaar (modified) framework, Assessment Tool for Evaluating Risk, Toxicity Estimation Software Tool, and OASIS. The MOA classifications from each scheme were first collapsed into one of 3 categories: non-specifically acting (i.e., narcosis), specifically acting, or nonclassifiable. Consensus rules were developed based on the degree of concordance among the 4 individual MOA classifications to attribute a consensus MOA to each chemical. A confidence rank was also assigned to the consensus MOA classification based on the degree of consensus. Overall, 40% of the chemicals were classified as narcotics, 17% as specifically acting, and 43% as unclassified. Sixty percent of chemicals had a medium to high consensus MOA assignment. When compared to empirical acute toxicity data, the general trend of specifically acting chemicals being more toxic is clearly observed for both fish and invertebrates but not for algae. EnviroTox is the first approach to establishing a high-level consensus across 4 computationally and structurally distinct MOA classification schemes. This consensus MOA classification provides both a transparent understanding of the variation between MOA classification schemes and an added certainty of the MOA assignment. In terms of regulatory relevance, a reliable understanding of MOA can provide information that can be useful for the prioritization (ranking) and risk assessment of chemicals. Environ Toxicol Chem 2019;38:2294-2304. © 2019 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals, Inc. on behalf of SETAC.

Application of an algal growth inhibition assay to determine distribution coefficients of benzalkonium ions between kaolinite and water

Benzalkonium compounds are widely used and found in environmental samples. Due to their amphiphilic nature, it is important to know sorption coefficients to account their bioavailability. However, currently available models describing their partitioning were developed using low molecular weight homologues and it cannot be ascertained whether they are applicable to their higher molecular weight homologues. Reasons for the scarcity of data on highly sorptive compounds include the lack of reliable quantification techniques for analyzing these chemicals at environmentally relevant levels. This study, therefore, reports on an algal growth inhibition assay-based method for the determination of kaolinite/water distribution coefficients for benzalkonium compounds at their environmentally relevant concentration range. Sorption to clay was computed using the difference between median effective concentration determined in a culture with kaolinite and that derived from a culture grown in standard medium. A kinetic model was used to account for uptake into algal cells and to calculate free concentrations. Due to the sensitivity of the algal species, Pseudokirchneriella subcapitata, it was possible to determine distribution coefficients below micromole per liter concentrations. The computed distribution coefficients showed a linear increase with number of carbon atoms in the alkyl chain up to 14. The proposed bioassay-based method should be applicable to determine distribution coefficients for highly hydrophobic chemicals and ionic liquids at a concentration range lower than typical analytical limits.

Passive Permeability of Planar Lipid Bilayers to Organic Anions

The membrane permeability P of organic ions was reported to be governed by the structure of the permeating molecule. Thus far, it is unclear whether the ion structure alters membrane partition or translocation proper across the membrane. Here, we obtained P values for 24 anionic compounds (18 concrete values, 6 upper limits) measuring the current that they carry through folded planar lipid bilayers. The P values range over more than 10 log units. Our measured permeability values correlate well (r = 0.95; logRMSE 0.74) with the hexadecane/water partition coefficients of the respective chemicals predicted by the COSMO-RS theory. Other attempts to predict P from the partition coefficient of the neutral molecule and from the solvation energy (Born energy) that opposes transfer into the membrane once the molecule is charged were unsuccessful. The uncertainties in assigning an effective radius to nonspherical molecules were much too large. The observation underlines that the actual structure of the molecules needs to be considered to predict partition and thus P by the solubility-diffusion model.

Screening tools for the bioconcentration potential of monovalent organic ions in fish

Currently the bioaccumulation potential of organic chemicals is assessed in a first tier approach via their octanol-water partition coefficient. This approach has been developed for neutral chemicals and cannot work for ionizable and ionic chemicals because the latter have different sorption-mechanisms and -preferences. Thus, suitable screening tools for the bioconcentration potential of ionic and ionizable chemicals need to be developed because it cannot be expected that these chemicals are non-bioaccumulative per se. Here, we present such screening tools for monovalent ions and ionizable chemicals based on calibrated sorption models for membrane lipids, structural proteins and albumin. The molecular descriptors used for these models arise from quantum chemical calculations and are based on COSMO-RS theory. When we applied our screening tools to 1839 preselected chemicals from the REACH registration data base, we identified 187 chemicals as potentially bioconcentrating (still ignoring any kind of metabolism). Among these were carbon and sulphur based aromatic and aliphatic acids mostly with a rather high molecular surface area. We hope that this outcome will trigger further research on ion specific sorption mechanisms and lead to a re-evaluation of the bioconcentration potential of ionic chemicals.