Phospholipophilicity of CxHyN(+) amines: chromatographic descriptors and molecular simulations for understanding partitioning into membranes

Research Priorities for the Environmental Risk Assessment of Per- and Polyfluorinated Substances

Per- and polyfluorinated substances (PFAS) are a group of thousands of ubiquitously applied persistent industrial chemicals. The field of PFAS environmental research is developing rapidly, but suffers from substantial biases toward specific compounds, environmental compartments, and organisms. The aim of our study was therefore to highlight current developments and to identify knowledge gaps and subsequent research needs that would contribute to a comprehensive environmental risk assessment for PFAS. To this end, we consulted the open literature and databases and found that knowledge of the environmental fate of PFAS is based on the analysis of <1% of the compounds categorized as PFAS. Moreover, soils and suspended particulate matter remain largely understudied. The bioavailability, bioaccumulation, and food web transfer studies of PFAS also focus on a very limited number of compounds and are biased toward aquatic biota, predominantly fish, and less frequently aquatic invertebrates and macrophytes. The available ecotoxicity data revealed that only a few PFAS have been well studied for their environmental hazards, and that PFAS ecotoxicity data are also strongly biased toward aquatic organisms. Ecotoxicity studies in the terrestrial environment are needed, as well as chronic, multigenerational, and community ecotoxicity research, in light of the persistency and bioaccumulation of PFAS. Finally, we identified an urgent need to unravel the relationships among sorption, bioaccumulation, and ecotoxicity on the one hand and molecular descriptors of PFAS chemical structures and physicochemical properties on the other, to allow predictions of exposure, bioaccumulation, and toxicity. Environ Toxicol Chem 2023;42:2302-2316. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

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.

Automated Coarse-Grained Mapping Algorithm for the Martini Force Field and Benchmarks for Membrane-Water Partitioning

With a view to high-throughput simulations, we present an automated system for mapping and parameterizing organic molecules for use with the coarse-grained Martini force field. The method scales to larger molecules and a broader chemical space than existing schemes. The core of the mapping process is a graph-based analysis of the molecule's bonding network, which has the advantages of being fast, general, and preserving symmetry. The parameterization process pays special attention to coarse-grained beads in aromatic rings. It also includes a method for building efficient and stable frameworks of constraints for molecules with structural rigidity. The performance of the method is tested on a diverse set of 87 neutral organic molecules and the ability of the resulting models to capture octanol-water and membrane-water partition coefficients. In the latter case, we introduce an adaptive method for extracting partition coefficients from free-energy profiles to take into account the interfacial region of the membrane. We also use the models to probe the response of membrane-water partitioning to the cholesterol content of the membrane.

Titratable Martini model for constant pH simulations

In this work, we deliver a proof of concept for a fast method that introduces pH effects into classical coarse-grained (CG) molecular dynamics simulations. Our approach is based upon the latest version of the popular Martini CG model to which explicit proton mimicking particles are added. We verify our approach against experimental data involving several different molecules and different environmental conditions. In particular, we compute titration curves, pH dependent free energies of transfer, and lipid bilayer membrane affinities as a function of pH. Using oleic acid as an example compound, we further illustrate that our method can be used to study passive translocation in lipid bilayers via protonation. Finally, our model reproduces qualitatively the expansion of the macromolecule dendrimer poly(propylene imine) as well as the associated pKa shift of its different generations. This example demonstrates that our model is able to pick up collective interactions between titratable sites in large molecules comprising many titratable functional groups.

Tissue Distribution of Several Series of Cationic Surfactants in Rainbow Trout (Oncorhynchus mykiss) Following Exposure via Water

Bioaccumulation assessment is important for cationic surfactants in light of their use in a wide variety of consumer products and industrial processes. Because they sorb strongly to natural surfaces and to cell membranes, their bioaccumulation behavior is expected to differ from other classes of chemicals. Divided over two mixtures, we exposed rainbow trout to water containing 10 alkyl amines and 2 quaternary alkylammonium surfactants for 7 days, analyzed different fish tissues for surfactant residues, and calculated the tissues' contribution to fish body burden. Mucus, skin, gills, liver, and muscle each contributed at least 10% of body burden for the majority of the test chemicals. This indicates that both sorption to external surfaces and systemic uptake contribute to bioaccumulation. In contrast to the analogue alkylamine bases, the permanently charged quaternary ammonium compounds accumulated mostly in the gills and were nearly absent in internal tissues, indicating that systemic uptake of the charged form of cationic surfactants is very slow. Muscle-blood distribution coefficients were close to 1 for all alkyl amines, whereas liver-blood distribution coefficients ranged from 13 to 90, suggesting that the dominant considerations for sorption in liver are different from those in blood and muscle. The significant fraction of body burden on external surfaces can have consequences for bioaccumulation assessment.

Influence of in Vitro Assay Setup on the Apparent Cytotoxic Potency of Benzalkonium Chlorides

The nominal concentration is generally used to express concentration–effect relationships in in vitro toxicity assays. However, the nominal concentration does not necessarily represent the exposure concentration responsible for the observed effect. Surfactants accumulate at interphases and likely sorb to in vitro system components such as serum protein and well plate plastic. The extent of sorption and the consequences of this sorption on in vitro readouts is largely unknown for these chemicals. The aim of this study was to demonstrate the effect of sorption to in vitro components on the observed cytotoxic potency of benzalkonium chlorides (BAC) varying in alkyl chain length (6–18 carbon atoms, C_6–18) in a basal cytotoxicity assay with the rainbow trout gill cell line (RTgill-W1). Cells were exposed for 48 h in 96-well plates to increasing concentration of BACs in exposure medium containing 0, 60 μM bovine serum albumin (BSA) or 10% fetal bovine serum (FBS). Before and after exposure, BAC concentrations in exposure medium were analytically determined. Based on freely dissolved concentrations at the end of the exposure, median effect concentrations (EC₅₀) decreased with increasing alkyl chain length up to 14 carbons. For BAC with alkyl chains of 12 or more carbons, EC₅₀’s based on measured concentrations after exposure in supplement-free medium were up to 25-times lower than EC₅₀’s calculated using nominal concentrations. When BSA or FBS was added to the medium, a decrease in cytotoxic potency of up to 22 times was observed for BAC with alkyl chains of eight or more carbons. The results of this study emphasize the importance of expressing the in vitro readouts as a function of a dose metric that is least influenced by assay setup to compare assay sensitivities and chemical potencies.

Membrane-Water Partition Coefficients to Aid Risk Assessment of Perfluoroalkyl Anions and Alkyl Sulfates

This study determined the sorption affinity to artificial phospholipid membranes ( K_MW) for series of perfluorinated carboxylates (PFCAs), perfluorinated sulfonates (PFSAs), alkyl sulfates (C _xSO₄), and 1-alkanesulfonates (C _xSO₃). A sorbent dilution assay with solid supported lipid membranes (SSLM) showed consistent CF₂ unit increments of 0.59, and CH₂ unit increments of 0.53, for the log K_MW of perfluorinated and hydrogenated anions, respectively. PFSAs sorbed 0.90 log units stronger than analogue PFCAs; C _xSO₄ sorbed 0.75 log units stronger than analogue C _xSO₃ anions. The log K_MW values for the octyl analogues increase in the order H(CH₂)₈SO₃^- (1.74) < H(CH₂)₈SO₄^- (2.58) < F(CF₂)₈CO₂^- (PFNA, 4.04) < F(CF₂)₈SO₃^- (PFOS, 4.88). Intrinsic partition ratios determined on a phospholipid coated HPLC column (IAM-HPLC) closely aligned with SSLM K_MW values. COSMO-RS based molecular calculations of K_MW aligned with SSLM K_MW values for hydrogenated anions with C₈-C₁₄ alkyl chains but strongly underestimated CF₂ and CH₂ unit increments for C₄-C₈ based anions. Dividing the critical narcotic membrane burden of 100 mmol/kg by the experimental K_MW predicts lethal baseline toxicity concentrations (LC_50,narc). The LC_50,narc coincides with the lowest reported acute LC₅₀ values for several anionic surfactants but were on average about an order of magnitude lower.

Sorption of Cationic Surfactants to Artificial Cell Membranes: Comparing Phospholipid Bilayers with Monolayer Coatings and Molecular Simulations

This study reports the distribution coefficient between phospholipid bilayer membranes and phosphate buffered saline (PBS) medium (D_MW,PBS) for 19 cationic surfactants. The method used a sorbent dilution series with solid supported lipid membranes (SSLMs). The existing SSLM protocol, applying a 96 well plate setup, was adapted to use 1.5 mL glass autosampler vials instead, which facilitated sampling and circumvented several confounding loss processes for some of the cationic surfactants. About 1% of the phospholipids were found to be detached from the SSLM beads, resulting in nonlinear sorption isotherms for compounds with log D_MW values above 4. Renewal of the medium resulted in linear sorption isotherms. D_MW values determined at pH 5.4 demonstrated that cationic surfactant species account for the observed D_MW,PBS. Log D_MW,PBS values above 5.5 are only experimentally feasible with lower LC-MS/MS detection limits and/or concentrated extracts of the aqueous samples. Based on the number of carbon atoms, dialkylamines showed a considerably lower sorption affinity than linear alkylamine analogues. These SSLM results closely overlapped with measurements on a chromatographic tool based on immobilized artificial membranes (IAM-HPLC) and with quantum-chemistry based calculations with COSMOmic. The SSLM data suggest that IAM-HPLC underestimates the D_MW of ionized primary and secondary alkylamines by 0.8 and 0.5 log units, respectively.

Which Molecular Features Affect the Intrinsic Hepatic Clearance Rate of Ionizable Organic Chemicals in Fish?

Greater knowledge of biotransformation rates for ionizable organic compounds (IOCs) in fish is required to properly assess the bioaccumulation potential of many environmentally relevant contaminants. In this study, we measured in vitro hepatic clearance rates for 50 IOCs using a pooled batch of liver S9 fractions isolated from rainbow trout (Oncorhynchus mykiss). The IOCs included four types of strongly ionized acids (carboxylates, phenolates, sulfonates, and sulfates), three types of strongly ionized bases (primary, secondary, tertiary amines), and a pair of quaternary ammonium compounds (QACs). Included in this test set were several surfactants and a series of beta-blockers. For linear alkyl chain IOC analogues, biotransformation enzymes appeared to act directly on the charged terminal group, with the highest clearance rates for tertiary amines and sulfates and no clearance of QACs. Clearance rates for C₁₂-IOCs were higher than those for C₈-IOC analogues. Several analogue series with multiple alkyl chains, branched alkyl chains, aromatic rings, and nonaromatic rings were evaluated. The likelihood of multiple reaction pathways made it difficult to relate all differences in clearance to specific molecular features the tested IOCs. Future analysis of primary metabolites in the S9 assay is recommended to further elucidate biotransformation pathways for IOCs in fish.