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

Predicting the Accumulation of Ionizable Pharmaceuticals and Personal Care Products in Aquatic and Terrestrial Organisms

The extent to which chemicals bioaccumulate in aquatic and terrestrial organisms represents a fundamental consideration for chemicals management efforts intended to protect public health and the environment from pollution and waste. Many chemicals, including most pharmaceuticals and personal care products (PPCPs), are ionizable across environmentally relevant pH gradients, which can affect their fate in aquatic and terrestrial systems. Existing mathematical models describe the accumulation of neutral organic chemicals and weak acids and bases in both fish and plants. Further model development is hampered, however, by a lack of mechanistic insights for PPCPs that are predominantly or permanently ionized. Targeted experiments across environmentally realistic conditions are needed to address the following questions: (1) What are the partitioning and sorption behaviors of strongly ionizing chemicals among species? (2) How does membrane permeability of ions influence bioaccumulation of PPCPs? (3) To what extent are salts and associated complexes with PPCPs influencing bioaccumulation? (4) How do biotransformation and other elimination processes vary within and among species? (5) Are bioaccumulation modeling efforts currently focused on chemicals and species with key data gaps and risk profiles? Answering these questions promises to address key sources of uncertainty for bioaccumulation modeling of ionizable PPCPs and related contaminants. Environ Toxicol Chem 2024;43:502-512. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

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.

Cocktails of pesticide residues in Prochilodus lineatus fish of the Salado River (South America): First record of high concentrations of polar herbicides

Muscle and viscera (gills-liver) of the fish Prochilodus lineatus were obtained from four sites of lower course of Salado river and one site at Santa Fe river near to its confluence with Salado river from Santa Fe (Argentina) between December 2021 and February 2022. Sediment samples were also obtained from the same sites. All samples were analyzed for pesticide residues following the QuEChERS method to quantify 136 compounds by UHPLC-ESI-MS/MS and GC-EI-MS/MS. Overall, muscle fish tissue showed very high concentrations (maximum concentrations detected) of the insecticide cypermethrin (204 μg/kg), polar herbicides (glyphosate; 187 μg/kg and its degradation product (aminomethylphosphonic acid) AMPA; 3116 μg/kg, and glufosinate-ammonium; 677 μg/kg), and the fungicide pyraclostrobin (50 μg/kg). In viscera samples, high values of cypermethrin (506 μg/kg), chlorpyrifos (78 μg/kg), and lambdacyhalothrin (73 μg/kg) were the main pesticides found. Mean residues concentrations detected among sites were not significantly different neither in muscle nor viscera of P. lineatus in most of the cases. Exceptionally, the southernmost studied site of the Lower Salado river showed significant differences in concentration of residues found in muscle, due to high concentrations of glyphosate and glufosinate-amonium (KW = 11.879 and KW = 13.013, respectively, P < 0.05). Other norther Lower Salado river site showed significant higher AMPA concentration in fish viscera than in the rest of the studied sites (KW = 12.86 P < 0.05). Some sediment samples showed low levels of herbicides such as glyphosate (24 μg/kg) and fungicides. However, the world highest levels of polar herbicides were recorded in fish muscle. The results of this study highlight the need for periodic monitoring due to the high concentration of pesticides and its potential risk in a very important commercial freshwater fish from Argentina, which is consumed locally and exported to other countries for human consumption.

Charged Small Molecule Binding to Membranes in MD Simulations Evaluated against NMR Experiments

Interactions of charged molecules with biomembranes regulate many of their biological activities, but their binding affinities to lipid bilayers are difficult to measure experimentally and model theoretically. Classical molecular dynamics (MD) simulations have the potential to capture the complex interactions determining how charged biomolecules interact with membranes, but systematic overbinding of sodium and calcium cations in standard MD simulations raises the question of how accurately force fields capture the interactions between lipid membranes and charged biomolecules. Here, we evaluate the binding of positively charged small molecules, etidocaine, and tetraphenylphosphonium to a phosphatidylcholine (POPC) lipid bilayer using the changes in lipid head-group order parameters. We observed that these molecules behave oppositely to calcium and sodium ions when binding to membranes: (i) their binding affinities are not overestimated by standard force field parameters, (ii) implicit inclusion of electronic polarizability increases their binding affinity, and (iii) they penetrate into the hydrophobic membrane core. Our results can be explained by distinct binding mechanisms of charged small molecules with hydrophobic moieties and monoatomic ions. The binding of the former is driven by hydrophobic effects, while the latter has direct electrostatic interactions with lipids. In addition to elucidating how different kinds of charged biomolecules bind to membranes, we deliver tools for further development of MD simulation parameters and methodology.

Human Apparent Volume of Distribution Predicts Bioaccumulation of Ionizable Organic Chemicals in Zebrafish Embryos

Chemicals with elevated bioaccumulation profiles present potential hazards to public health and the environment. Ionizable organic compounds (IOCs) increasingly represent a large proportion of commercial chemicals; however, historical approaches for bioaccumulation determinations are mainly developed for neutral chemicals, which were not appropriate for IOCs. Herein, we employed the zebrafish embryo, a common vertebrate model in environmental and biomedical studies, to elucidate toxicokinetics and bioconcentration of eight IOCs with diverse physicochemical properties and pharmacokinetic parameters. At an environmentally relevant pH (7.5), most IOCs exhibited rapid uptake and depuration in zebrafish, suggesting the ionized forms of IOCs are readily bioavailable. Bioconcentration factors (BCF) of these IOCs ranged from 0.0530 to 250 L·kg^-1 wet weight. The human pharmacokinetic proportionality factor, apparent volume of distribution (V_D), better predicted the BCF of selected IOCs than more commonly used hydrophobicity-based parameters (e.g., pH-dependent octanol-water distribution ratio, D_ow). Predictive bioaccumulation models for IOCs were constructed and validated using V_D alone or with D_ow. Significant relationships between fish BCF and human V_D, which is readily available for pharmaceuticals, highlighted the utility of biologically based "read-across" approaches for predicting bioaccumulative potential of IOCs. Our novel findings thus provided an understanding of the partitioning behavior and improved predictive bioconcentration modeling for IOCs.

Evaluating phospholipid- and protein-water partitioning of two groups of chemicals of emerging concern: Diastereo- and enantioselectivity

The partitioning between phospholipids/proteins and water can be used to predict the bioaccumulation potential of chemicals with better accuracy compared with n-octanol-water partition coefficient. However, such partitioning is poorly understood for chiral chemicals, many of which exhibit differential bioaccumulation and toxicity potential between enantiomers. In this study, the enantiospecific liposome-water and bovine serum albumin (BSA)-water partition coefficients (Klip/w and KBSA/w, determined at 25 ℃ and 37 ℃, respectively) were measured by equilibrium dialysis for α-, β-, and γ-hexabromocyclododecane (HBCD) and three β-blockers (propranolol, metoprolol, and sotalol). Raman and fluorescence analyses and molecular docking were conducted to provide additional insights into the partitioning process. Results showed α- and β-HBCD displayed stronger enantioselective partitioning to liposomes with the (-)-form, while (-)-α-HBCD, R-(+)-propranolol, R-(+)-metoprolol, and E2-sotalol favored partitioning to BSA compared with their antipodes. Raman spectra revealed α- and γ-HBCD enhanced and reduced the organization of liposome acyl chains, respectively, and polar interactions enhanced the liposome partitioning of β-blockers. Fluorescence spectra indicated the changed tryptophan microenvironment might influence the BSA steric effect toward HBCD, and electrostatic interactions dominated the formation of BSA-β-blocker complexes. Molecular docking results supported the difference in the thermodynamic nature of interaction between the studied enantiomers and BSA.

Concentration of Antioxidant Compounds from Calendula officinalis through Sustainable Supercritical Technologies, and Computational Study of Their Permeability in Skin for Cosmetic Use

The growing interest in the cosmetic industry in using compounds of natural and sustainable origin that are safe for humans is encouraging the development of processes that can satisfy these needs. Chlorogenic acid (CHA), caffeic acid (CAF) and ferulic acid (FA) are three compounds widely used within the cosmetic industry due to their functionalities as antioxidants, collagen modifiers or even as radiation protectors. In this work, two advanced separation techniques with supercritical CO2 are used to obtain these three compounds from Calendula officinalis, and these are then evaluated using a computational skin permeability model. This model is encompassed by the COSMO-RS model, the calculations of which make it possible to study the behaviour of the compounds in the epidermis. The results show that both CAF and FA are retained in the stratum corneum, while CHA manages to penetrate to the stratum spinosum. These compounds were concentrated by antisolvent fractionation with super-critical CO2 using a Response Surface Methodology to study the effect of pressure and CO2 flow rate. CHA, CAF and FA were completely retained in the precipitation vessel, with concentrations between 40% and 70% greater than in the original extract. The conditions predicted that the optimal overall yield and enrichment achieved would be 153 bar and 42 g/min.

Screening of 6000 Compounds for Uncoupling Activity: A Comparison Between a Mechanistic Biophysical Model and the Structural Alert Profiler Mitotox

Protonophoric uncoupling of phosphorylation is an important factor when assessing chemicals for their toxicity, and has recently moved into focus in pharmaceutical research with respect to the treatment of diseases such as cancer, diabetes, or obesity. Reliably identifying uncoupling activity is thus a valuable goal. To that end, we screened more than 6000 anionic compounds for in vitro uncoupling activity, using a biophysical model based on ab initio COSMO-RS input parameters with the molecular structure as the only external input. We combined these results with a model for baseline toxicity (narcosis). Our model identified more than 1250 possible uncouplers in the screening dataset, and identified possible new uncoupler classes such as thiophosphoric acids. When tested against 423 known uncouplers and 612 known inactive compounds in the dataset, the model reached a sensitivity of 83% and a specificity of 96%. In a direct comparison, it showed a similar specificity than the structural alert profiler Mitotox (97%), but much higher sensitivity than Mitotox (47%). The biophysical model thus allows for a more accurate screening for uncoupling activity than existing structural alert profilers. We propose to use our model as a complementary tool to screen large datasets for protonophoric uncoupling activity in drug development and toxicity assessment.

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.

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.

Development of liposome/water partition coefficients predictive models for neutral and ionogenic organic chemicals

Membrane/water partition coefficient (K_m/w) is a vital parameter used to characterize the membrane permeability of compounds. Considering the K_m/w value is difficult to observe experimentally for real biological membranes, liposome/water partition coefficient (K_lip/w) is employed to approximate K_m/w. Here, quantitative structure property relationship (QSPR) models for logK_lip/w of the neutral organic chemicals and the neutral form of ionogenic organic chemicals (IOCs) (logK_{lip/w-neutral}), ionic form of IOCs (logK_lip/w-ionic), the speciation-corrected liposome-water distribution ratios at a pH = 7.40 (logD_{lip/w-(pH=7.40)}) were developed. In the modeling, two modeling methods (multiple linear regressions (MLR) and k-nearest neighbor (kNN)) were used. The predictive variables employed here could be calculated from the molecular structure directly. For logK_{lip/w-neutral} and logD_{lip/w-(pH=7.40)}, the logK_OW and logD_OW-based, non-logK_OW and non-logD_OW-based kNN-QSPR and MLR-QSPR models were developed, respectively. The evaluation results implied that the predictive performance of kNN-QSPR models is better than that of MLR-QSPR models. For logK_lip/w-ionic, only one acceptable MLR-QSPR model was developed for cation and anion, respectively. The model quality of the derived models was evaluated following the OECD QSPR models validation guideline. The determination coefficient (R²), leave-one-out cross validation Q² (Q²_LOO) and bootstrapping coefficient (Q²_BOOT), the external validation coefficient (Q²_EXT) of all the models met the acceptable criteria (Q² > 0.600, R² > 0.700); while the root-mean-square error (RMSE) range from 0.351 to 0.857. All the results implied that the models had good goodness-of-fit, robustness and predictive ability. Therefore, the developed models could be used to fill the data gap for substances within the applicability domain on their missing logK_{lip/w-neutral}, logK_lip/w-ionic, logD_{lip/w-(pH=7.40)} values.

Combined Ion-Trapping and Mass Balance Models To Describe the pH-Dependent Uptake and Toxicity of Acidic and Basic Pharmaceuticals in Zebrafish Embryos ( Danio rerio)

The aim of the current study was to understand and develop models to predict the pH-dependent toxicity of ionizable pharmaceuticals in embryos of the zebrafish Danio rerio. We found a higher uptake and toxicity with increasing neutral fraction of acids (diclofenac, genistein, naproxen, torasemide, and warfarin) and bases (metoprolol and propranolol). Simple mass balance models accounting for the partitioning to lipids and proteins in the zebrafish embryo were found to be suitable to predict the bioconcentration after 96 h of exposure if pH values did not differ much from the internal pH of 7.55. For other pH values, a kinetic ion-trap model for the zebrafish embryo explained the pH dependence of biouptake and toxicity. The total internal lethal concentrations killing 50% of the zebrafish embryos (ILC₅₀) were calculated from the measured BCF and LC₅₀. The resulting ILC₅₀ were independent of external pH. Critical membrane concentrations were deduced by an internal mass balance model, and apart from diclofenac, whose specific toxicity in fish had already been established, all pharmaceuticals were confirmed to act as baseline toxicants in zebrafish.

Phospholipid Levels Predict the Tissue Distribution of Poly- and Perfluoroalkyl Substances in a Marine Mammal

Exposure to poly- and perfluoroalkyl substances (PFASs) has been linked to many negative health impacts in humans and wildlife. Unlike neutral hydrophobic organic pollutants, many PFASs are ionic and have been hypothesized to accumulate in both phospholipids and protein-rich tissues. Here we investigate the role of phospholipids for PFAS accumulation by analyzing associations among concurrent measurements of phospholipid, total protein, total lipid and 24 PFASs in the heart, muscle, brain, kidney, liver, blubber, placenta and spleen of North Atlantic pilot whales (Globicephala melas). The sum of 24 PFASs ( ∑ 24 PFAS ) was highest in the liver (median 260 ng g^-1; interquartile range (IQR) 216-295 ng g^-1) and brain (86.0; IQR 54.5-91.3 ng g^-1), while phospholipid levels were highest in brain. The relative abundance of PFASs in the brain greatly increases with carbon chain lengths of 10 or greater, suggesting shorter-chained compounds may cross the blood-brain barrier less efficiently. Phospholipids were significant predictors of the tissue distribution of the longest-chained PFASs: perfluorodecanesulfonate (PFDS), perfluorododecanoate (PFDoA), perfluorotridecanoate (PFTrA), and perfluorotetradecanoic acid (PFTA) (r_s = 0.5-0.6). In all tissues except the brain, each 1 mg g^-1 increase in phospholipids led to a 12%-25% increase in the concentration of each PFAS. We conclude that partitioning to phospholipids is an important mechanism of bioaccumulation for long-chained PFASs in marine mammals.