MS2Tox Machine Learning Tool for Predicting the Ecotoxicity of Unidentified Chemicals in Water by Nontarget LC-HRMS

Prioritization strategies for non-target screening in environmental samples by chromatography - High-resolution mass spectrometry: A tutorial

Non-target screening (NTS) using chromatography coupled to high-resolution mass spectrometry (HRMS), has become fundamental for detecting and prioritizing chemicals of emerging concern (CECs) in complex environmental matrices. The vast number of generated features (m/z, retention time, and intensity) necessitate effective prioritization strategies to identify environmentally and toxicologically relevant CECs. Since compound identification remains a major bottleneck in NTS, prioritization is critical to focus identification efforts where they matter most. This tutorial presents seven prioritization strategies: (1) Target and suspect screening for identifying known or suspected compounds using reference libraries. (2) Data quality filtering to apply quality control measures to reduce noise and the number of false positives. (3) Chemistry-driven prioritization using HRMS data properties to prioritize specific compound classes (e.g., halogenated substances, transformation products). (4) Process-driven - using spatial, temporal, or process-based comparisons (pre- and post-technical processes) to identify key features. (5) Effect-Directed Analysis (EDA) and Virtual Effect-Directed Analysis (vEDA) prioritization to link chemical features to biological effects. (6) Prediction-based prioritization such as quantitative structure-property relationships (QSPR) and machine learning to estimate risk or concentration levels, and (7) Pixel- or tile-based analysis where the chromatographic image (2D data) is used to pin-point regions of interest or for comparison of larger sample sets. By integrating these prioritization strategies, this tutorial provides a structured foundation to evaluate both identified and unidentified features, prioritize high-risk compounds, and advance environmental risk assessment and regulatory decision-making.

Xenometabolomics in Ecotoxicology: Concepts and Applications

Nontargeted high-resolution mass spectrometry (HRMS) allows for the characterization of a large fraction of the exposome, i.e., the entirety of chemicals an organism is exposed to, and helps detect important exogenous chemical compounds that could be key drivers of toxicological impact. Along with these chemical compounds occur endogenous metabolites that are essential for the health of the host organism. Chemical compounds derived from the biotransformation of xenobiotics present in the exposome are referred to as the xenometabolome, while endogenous metabolites derived from the host organism are referred to as the endometabolome. Recent advancements in HRMS technology allow for the detection of chemical features of biological and ecological importance in the context of chemical safety assessments with unprecedented sensitivity and resolution. In this perspective, we highlight the application of HRMS-based metabolomics of organisms in the context of ecotoxicology, the complexity of comprehensively characterizing the endometabolome, and distinguishing chemical compounds of the xenometabolome.

Prioritization of Unknown LC-HRMS Features Based on Predicted Toxicity Categories

Complex environmental samples contain a diverse array of known and unknown constituents. While liquid chromatography coupled with high-resolution mass spectrometry (LC-HRMS) nontargeted analysis (NTA) has emerged as an essential tool for the comprehensive study of such samples, the identification of individual constituents remains a significant challenge, primarily due to the vast number of detected features in each sample. To address this, prioritization strategies are frequently employed to narrow the focus to the most relevant features for further analysis. In this study, we developed a novel prioritization strategy that directly links fragmentation and chromatographic data to aquatic toxicity categories, bypassing the need for identification of individual compounds. Given that features are not always well-characterized through fragmentation, we created two models: (1) a Random Forest Classification (RFC) model, which classifies fish toxicity categories based on MS1, retention, and fragmentation data─expressed as cumulative neutral losses (CNLs)─when fragmentation information is available, and (2) a Kernel Density Estimation (KDE) model that relies solely on retention time and MS1 data when fragmentation is absent. Both models demonstrated accuracy comparable to that of structure-based prediction methods. We further tested the models on a pesticide mixture in a tea extract measured by LC-HRMS, where the CNL-based RFC model achieved 0.76 accuracy and the KDE model reached 0.61, showcasing their robust performance in real-world applications.

FT-GNN Tool for Bridging HRMS Features and Bioactivity: Uncovering Unidentified Estrogen Receptor Agonists in Sewage

Identifying primary estrogen receptor (ER) agonists in municipal sewage is essential for ensuring the health of aquatic environments. Given the complex and variable chemical composition of sewage, the predominant ER agonists remain unclear. High-resolution mass spectrometry (HRMS)-based models have been developed to predict compound bioactivity in complex matrices, but further optimization is needed to effectively bridge HRMS features with ER agonists. To address this challenge, an FT-GNN (fragmentation tree-based graph neural network) model was proposed. Given limited data and class imbalance, data augmentation was performed using model predictions within the applicability domain (AD) and oversampling technique (OTE). Model development results demonstrated that integrating the FT-GNN with data augmentation improved the balanced accuracy (bACC) value by 6%-31%. The developed model, with a high bACC to identify more true ER agonists, efficiently classified tens of thousands of unidentified HRMS features in sewage, reducing postprocessing workload in nontargeted screening. Analysis of ER agonist transformation during sewage treatment revealed the anaerobic stage as key to both their removal and formation. Estrogenic effect balance analysis suggests that α-E2 and 9,11-didehydroestriol may be two previously overlooked key ER agonists. Collectively, the development and application of the FT-GNN model are crucial advancements toward credible tracking and efficient control of estrogenic risks in water.

Advancing the Spatiotemporal Dimension of Wildlife-Pollution Interactions

Chemical pollution is one of the fastest-growing agents of global change. Numerous pollutants are known to disrupt animal behavior, alter ecological interactions, and shift evolutionary trajectories. Crucially, both chemical pollutants and individual organisms are nonrandomly distributed throughout the environment. Despite this fact, the current evidence for chemical-induced impacts on wildlife largely stems from tests that restrict organism movement and force homogeneous exposures. While such approaches have provided pivotal ecotoxicological insights, they overlook the dynamic spatiotemporal interactions that shape wildlife-pollution relationships in nature. Indeed, the seemingly simple notion that pollutants and animals move nonrandomly in the environment creates a complex of dynamic interactions, many of which have never been theoretically modeled or experimentally tested. Here, we conceptualize dynamic interactions between spatiotemporal variation in pollutants and organisms and highlight their ecological and evolutionary implications. We propose a three-pronged approach-integrating in silico modeling, laboratory experiments that allow movement, and field-based tracking of free-ranging animals-to bridge the gap between controlled ecotoxicological studies and real-world wildlife exposures. Advances in telemetry, remote sensing, and computational models provide the necessary tools to quantify these interactions, paving the way for a new era of ecotoxicology that accounts for spatiotemporal complexity.

Advancing non-target analysis of emerging environmental contaminants with machine learning: Current status and future implications

Emerging environmental contaminants (EECs) such as pharmaceuticals, pesticides, and industrial chemicals pose significant challenges for detection and identification due to their structural diversity and lack of analytical standards. Traditional targeted screening methods often fail to detect these compounds, making non-target analysis (NTA) using high-resolution mass spectrometry (HRMS) essential for identifying unknown or suspected contaminants. However, interpreting the vast datasets generated by HRMS is complex and requires advanced data processing techniques. Recent advancements in machine learning (ML) models offer great potential for enhancing NTA applications. As such, we reviewed key developments, including optimizing workflows using computational tools, improved chemical structure identification, advanced quantification methods, and enhanced toxicity prediction capabilities. It also discusses challenges and future perspectives in the field, such as refining ML tools for complex mixtures, improving inter-laboratory validation, and further integrating computational models into environmental risk assessment frameworks. By addressing these challenges, ML-assisted NTA can significantly enhance the detection, quantification, and evaluation of EECs, ultimately contributing to more effective environmental monitoring and public health protection.

Machine Learning-based Classification for the Prioritization of Potentially Hazardous Chemicals with Structural Alerts in Nontarget Screening

Nontarget screening (NTS) with liquid chromatography high-resolution mass spectrometry (LC-HRMS) is commonly used to detect unknown organic micropollutants in the environment. One of the main challenges in NTS is the prioritization of relevant LC-HRMS features. A novel prioritization strategy based on structural alerts to select NTS features that correspond to potentially hazardous chemicals is presented here. This strategy leverages raw tandem mass spectra (MS2) and machine learning models to predict the probability that NTS features correspond to chemicals with structural alerts. The models were trained on fragments and neutral losses from the experimental MS2 data. The feasibility of this approach is evaluated for two groups: aromatic amines and organophosphorus structural alerts. The neural network classification model for organophosphorus structural alerts achieved an Area Under the Curve of the Receiver Operating Characteristics (AUC-ROC) of 0.97 and a true positive rate of 0.65 on the test set. The random forest model for the classification of aromatic amines achieved an AUC-ROC value of 0.82 and a true positive rate of 0.58 on the test set. The models were successfully applied to prioritize LC-HRMS features in surface water samples, showcasing the high potential to develop and implement this approach further.

From High Resolution Tandem Mass Spectrometry to Pollutant Toxicity AI-Based Prediction: A Case Study of 7 Endocrine Disruptors Endpoints

Based on high-resolution mass spectrometry (HRMS), nontarget analysis (NTA) can rapidly identify and characterize numerous hazardous substances in complex environmental samples. However, the intricate identification process often results in the underutilization of many mass spectrometry features. Even when chemical structures are identified, their toxicological effects and health outcomes may remain unknown. To address these challenges, this study introduces MSFragTox, a novel approach that leverages the rich fragmentation spectra inherent in high resolution tandem mass spectrometry (MS/MS) to directly predict toxicity. This method integrates MS/MS data with high-throughput screening (HTS) assays, focusing on seven endocrine disruption-related endpoints from Tox21, and uses MS-derived fingerprints: substructure fragmentation probability vectors to construct toxicity predictions using machine learning algorithms. The best model demonstrated robust performance with an average area under the receiver operating characteristic curve (AUROC) of 0.845 on the test set, outperforming models based on traditional molecular fingerprints and descriptors. Additionally, a web client (http://ms.envwind.site:8500) is provided for users to screen toxicity based on chemical MS/MS data. Furthermore, in-depth analyses of commonalities and differences in substructures reveal the mechanisms underlying across toxicity endpoints. Using MSFragTox, we validated the potential endocrine-disrupting effects of substances corresponding to MS/MS from real samples, highlighting the feasibility of directly studying toxicity through MS/MS and its potential applications in risk prediction and early warning for environmental samples.

Communicating with Stakeholders to Identify High-Impact Research Directions for Non-Targeted Analysis

Non-targeted analysis (NTA) using high-resolution mass spectrometry without defined chemical targets has the potential to expand and improve chemical monitoring in many fields. Despite rapid advancements within the research community, NTA methods and data remain underutilized by many potential beneficiaries. To better understand barriers toward widespread adoption, the Best Practices for Non-Targeted Analysis (BP4NTA) working group conducted focus group meetings and follow-up surveys with scientists (n = 61) from various sectors (e.g., drinking water utilities, epidemiologists, n = 9) where NTA is expected to provide future value. Meeting participants included producers and end-users of NTA data with a wide range of familiarity with NTA methods and outputs. Discussions focused on identifying specific barriers that limit adoption and on setting NTA product development priorities. Stated priorities fell into four major categories: 1) education and training materials; 2) QA/QC frameworks and study design guidance; 3) accessible compound databases and libraries; and 4) NTA data linkages with chemical fate and toxicity information. Based on participant feedback, this manuscript proposes research directions, such as standardization of training materials, that BP4NTA and other institutions can pursue to expand NTA use in various application scenarios and decision contexts.

MLinvitroTox reloaded for high-throughput hazard-based prioritization of high-resolution mass spectrometry data

MLinvitroTox is an automated Python pipeline developed for high-throughput hazard-driven prioritization of toxicologically relevant signals detected in complex environmental samples through high-resolution tandem mass spectrometry (HRMS/MS). MLinvitroTox is a machine learning (ML) framework comprising 490 independent XGBoost classifiers trained on molecular fingerprints from chemical structures and target-specific endpoints from the ToxCast/Tox21 invitroDBv4.1 database. For each analyzed HRMS feature, MLinvitroTox generates a 490-bit bioactivity fingerprint used as a basis for prioritization, focusing the time-consuming molecular identification efforts on features most likely to cause adverse effects. The practical advantages of MLinvitroTox are demonstrated for groundwater HRMS data. Among the 874 features for which molecular fingerprints were derived from spectra, including 630 nontargets, 185 spectral matches, and 59 targets, around 4% of the feature/endpoint relationship pairs were predicted to be active. Cross-checking the predictions for targets and spectral matches with invitroDB data confirmed the bioactivity of 120 active and 6791 nonactive pairs while mislabeling 88 active and 56 non-active relationships. By filtering according to bioactivity probability, endpoint scores, and similarity to the training data, the number of potentially toxic features was reduced by at least one order of magnitude. This refinement makes the analytical confirmation of the toxicologically most relevant features feasible, offering significant benefits for cost-efficient chemical risk assessment.Scientific Contribution:In contrast to the classical ML-based approaches for toxicity prediction, MLinvitroTox predicts bioactivity for HRMS features (i.e., distinct m/z signals) based on MS2 fragmentation spectra rather than the chemical structures from the identified features. While the original proof of concept study was accompanied by the release of a MLinvitroTox v1 KNIME workflow, in this study, we release a Python MLinvitroTox v2 package, which, in addition to automation, expands functionality to include predicting toxicity from structures, cleaning up and generating chemical fingerprints, customizing models, and retraining on custom data. Furthermore, as a result of improvements in bioactivity data processing, realized in the concurrently released pytcpl Python package for the custom processing of invitroDBv4.1 input data used for training MLinvitroTox, the current release introduces enhancements in model accuracy, coverage of biological mechanistic targets, and overall interpretability.

Progress, applications, and challenges in high-throughput effect-directed analysis for toxicity driver identification - is it time for HT-EDA?

The rapid increase in the production and global use of chemicals and their mixtures has raised concerns about their potential impact on human and environmental health. With advances in analytical techniques, in particular, high-resolution mass spectrometry (HRMS), thousands of compounds and transformation products with potential adverse effects can now be detected in environmental samples. However, identifying and prioritizing the toxicity drivers among these compounds remain a significant challenge. Effect-directed analysis (EDA) emerged as an important tool to address this challenge, combining biotesting, sample fractionation, and chemical analysis to unravel toxicity drivers in complex mixtures. Traditional EDA workflows are labor-intensive and time-consuming, hindering large-scale applications. The concept of high-throughput (HT) EDA has recently gained traction as a means of accelerating these workflows. Key features of HT-EDA include the combination of microfractionation and downscaled bioassays, automation of sample preparation and biotesting, and efficient data processing workflows supported by novel computational tools. In addition to microplate-based fractionation, high-performance thin-layer chromatography (HPTLC) offers an interesting alternative to HPLC in HT-EDA. This review provides an updated perspective on the state-of-the-art in HT-EDA, and novel methods/tools that can be incorporated into HT-EDA workflows. It also discusses recent studies on HT-EDA, HT bioassays, and computational prioritization tools, along with considerations regarding HPTLC. By identifying current gaps in HT-EDA and proposing new approaches to overcome them, this review aims to bring HT-EDA a step closer to monitoring applications.

Advances in environmental analysis of high molecular weight disinfection byproducts

The disinfection of drinking water, while critical for public health, leads to the formation of disinfection byproducts (DBPs). Toxicological and epidemiological studies have demonstrated that exposure to disinfected water samples may pose adverse effects on human health. Recent research highlights the potential greater toxicity contribution of DBP fractions with high molecular weight (MW) (with more than two carbon atoms) compared to regulated low MW DBPs, emphasizing the need for advanced analytical techniques to identify and characterize these fractions. In this review, we summarize different analytical techniques for indirectly assessing DBP precursors and directly analyzing DBPs, discussing their advantages and limitations. Additionally, since identifying DBP toxicity agents in complex water mixtures is crucial for further optimizing water disinfection and controlling DBP formation, key DBP identification methods based on both chemical and bioassay metrics are also included and discussed. Finally, we highlight three important aspects for the future development of analytical methods to enhance the understanding of high MW DBP formation.

Determining high priority disinfection byproducts based on experimental aquatic toxicity data and predictive models: Virtual screening and in vivo study

Only about 100 disinfection byproducts (DBPs) have been tested for their potential aquatic toxicity. It is not known which specific DBPs, DBP main groups, and DBP subgroups are more toxic due to the lack of experimental toxicity data. Herein, high priority specific DBPs, DBP main groups, DBP subgroups, most sensitive model aquatic species, potential PBT and PMT (persistent, bioaccumulative/mobile, and toxic) DBPs were virtually screened for 1187 updated DBPs inventory. Priority setting based on experimental and predicted acute and chronic aquatic toxicity data found that the aromatic and alicyclic DBPs in four DBPs main groups showed high priority because larger proportions of aromatic and alicyclic DBPs are in high hazard categories (i.e. Acute and/or Chronic Toxic-1 or Toxic-2) according to the criteria in GHS system compared to the aliphatic and heterocyclic DBPs. The halophenols, estrogen-DBPs, nonhalogenated esters, and nonhalogenated aldehydes were recognized as high priority DBPs subgroups. For specific DBPs, 19 and 31 DBPs should be highly concerned in the future study because both acute and chronic toxicity of those DBPs to all of the three aquatic life (algae, Daphnia magna, fish) were classified as Toxic-1 and Toxic-2, respectively. The Daphnia magna and algae were sensitive to the acute toxicity of DBPs, while the fish and Daphnia magna were sensitive to the chronic toxicity of DBPs. One potential PBT (Tetrachlorobisphenol A) and four potential PMT DBPs were identified. For verification, the acute toxicity of four DBPs on three aquatic organism were performed, and their tested acute toxicity data to three aquatic organisms were consistent with the predictions. Our results could be beneficial to government regulators to adopt effective measures to limit the discharge of high priority DBPs and help the scientific community to develop or improve disinfection processes to reduce the production of high priority DBPs.

Quantification Approaches in Non-Target LC/ESI/HRMS Analysis: An Interlaboratory Comparison

Nontargeted screening (NTS) utilizing liquid chromatography electrospray ionization high-resolution mass spectrometry (LC/ESI/HRMS) is increasingly used to identify environmental contaminants. Major differences in the ionization efficiency of compounds in ESI/HRMS result in widely varying responses and complicate quantitative analysis. Despite an increasing number of methods for quantification without authentic standards in NTS, the approaches are evaluated on limited and diverse data sets with varying chemical coverage collected on different instruments, complicating an unbiased comparison. In this interlaboratory comparison, organized by the NORMAN Network, we evaluated the accuracy and performance variability of five quantification approaches across 41 NTS methods from 37 laboratories. Three approaches are based on surrogate standard quantification (parent-transformation product, structurally similar or close eluting) and two on predicted ionization efficiencies (RandFor-IE and MLR-IE). Shortly, HPLC grade water, tap water, and surface water spiked with 45 compounds at 2 concentration levels were analyzed together with 41 calibrants at 6 known concentrations by the laboratories using in-house NTS workflows. The accuracy of the approaches was evaluated by comparing the estimated and spiked concentrations across quantification approaches, instrumentation, and laboratories. The RandFor-IE approach performed best with a reported mean prediction error of 15× and over 83% of compounds quantified within 10× error. Despite different instrumentation and workflows, the performance was stable across laboratories and did not depend on the complexity of water matrices.

Advancing the Prediction of MS/MS Spectra Using Machine Learning

Tandem mass spectrometry (MS/MS) is an important tool for the identification of small molecules and metabolites where resultant spectra are most commonly identified by matching them with spectra in MS/MS reference libraries. While popular, this strategy is limited by the contents of existing reference libraries. In response to this limitation, various methods are being developed for the in silico generation of spectra to augment existing libraries. Recently, machine learning and deep learning techniques have been applied to predict spectra with greater speed and accuracy. Here, we investigate the challenges these algorithms face in achieving fast and accurate predictions on a wide range of small molecules. The challenges are often amplified by the use of generic machine learning benchmarking tactics, which lead to misleading accuracy scores. Curating data sets, only predicting spectra for sufficiently high collision energies, and working more closely with experimental mass spectrometrists are recommended strategies to improve overall prediction accuracy in this nuanced field.

Evaluation of Nontargeted Mass Spectral Data Acquisition Strategies for Water Analysis and Toxicity-Based Feature Prioritization by MS2Tox

The machine-learning tool MS2Tox can prioritize hazardous nontargeted molecular features in environmental waters, by predicting acute fish lethality of unknown molecules based on their MS2 spectra, prior to structural annotation. It has yet to be investigated how the extent of molecular coverage, MS2 spectra quality, and toxicity prediction confidence depend on sample complexity and MS2 data acquisition strategies. We compared two common nontargeted MS2 acquisition strategies with liquid chromatography high-resolution mass spectrometry for structural annotation accuracy by SIRIUS+CSI:FingerID and MS2Tox toxicity prediction of 191 reference chemicals spiked to LC-MS water, groundwater, surface water, and wastewater. Data-dependent acquisition (DDA) resulted in higher rates (19-62%) of correct structural annotations among reference chemicals in all matrices except wastewaters, compared to data-independent acquisition (DIA, 19-50%). However, DIA resulted in higher MS2 detection rates (59-84% DIA, 37-82% DDA), leading to higher true positive rates for spectral library matching, 40-73% compared to 34-72%. DDA resulted in higher MS2Tox toxicity prediction accuracy than DIA, with root-mean-square errors of 0.62 and 0.71 log-mM, respectively. Given the importance of MS2 spectral quality, we introduce a "CombinedConfidence" score to convey relative confidence in MS2Tox predictions and apply this approach to prioritize potentially ecotoxic nontargeted features in environmental waters.

Non-target analysis of Danish wastewater treatment plant effluent: Statistical analysis of chemical fingerprinting as a step toward a future monitoring tool

In an attempt to discover and characterize the plethora of xenobiotic substances, this study investigates chemical compounds released into the environment with wastewater effluents. A novel non-targeted screening methodology based on ultra-high resolution Orbitrap mass spectrometry and nanoflow ultra-high performance liquid chromatography together with a newly optimized data-processing pipeline were applied to effluent samples from two state-of-the-art and one small wastewater treatment facility. In total, 785 molecular structures were obtained, of which 38 were identified as single compounds, while 480 structures were identified at a putative level. Most of these substances were therapeutics and drugs, present as parent compounds and metabolites. Using R packages Phyloseq and MetacodeR, originally developed for bioinformatics, significant differences in xenobiotic presence in the wastewater effluents between the three sites were demonstrated.

Prioritization, Identification, and Quantification of Emerging Contaminants in Recycled Textiles Using Non-Targeted and Suspect Screening Workflows by LC-ESI-HRMS

Recycled textiles are becoming widely available to consumers as manufacturers adopt circular economy principles to reduce the negative impact of garment production. Still, the quality of the source material directly impacts the final product, where the presence of harmful chemicals is of utmost concern. Here, we develop a risk-based suspect and non-targeted screening workflow for the detection, identification, and prioritization of the chemicals present in consumer-based recycled textile products after manufacture and transport. We apply the workflow to characterize 13 recycled textile products from major retail outlets in Sweden. Samples were extracted and analyzed by liquid chromatography coupled with high-resolution mass spectrometry (LC-HRMS). In positive and negative ionization mode, 20,119 LC-HRMS features were detected and screened against persistent, mobile, and toxic (PMT) as well as other textile-related chemicals. Six substances were matched with PMT substances that are regulated in the European Union (EU) with a Level 2/3 confidence. Forty-three substances were confidently matched with textile-related chemicals reported for use in Sweden. For estimating the relative priority score, aquatic toxicity and concentrations were predicted for 7416 features with tandem mass spectra (MS2) and used to rank the non-targeted features. The top 10 substances were evaluated due to elevated environmental risk linked to the recycling process and potential release at end-of-life.

Deep Learning Bridged Bioactivity, Structure, and GC-HRMS-Readable Evidence to Decipher Nontarget Toxicants in Sediments

Identifying causative toxicants in mixtures is critical, but this task is challenging when mixtures contain multiple chemical classes. Effect-based methods are used to complement chemical analyses to identify toxicants, yet conventional bioassays typically rely on an apical and/or single endpoint, providing limited diagnostic potential to guide chemical prioritization. We proposed an event-driven taxonomy framework for mixture risk assessment that relied on high-throughput screening bioassays and toxicant identification integrated by deep learning. In this work, the framework was evaluated using chemical mixtures in sediments eliciting aryl-hydrocarbon receptor activation and oxidative stress response. Mixture prediction using target analysis explained <10% of observed sediment bioactivity. To identify additional contaminants, two deep learning models were developed to predict fingerprints of a pool of bioactive substances (event driver fingerprint, EDFP) and convert these candidates to MS-readable information (event driver ion, EDION) for nontarget analysis. Two libraries with 121 and 118 fingerprints were established, and 247 bioactive compounds were identified at confidence level 2 or 3 in sediment extract using GC-qToF-MS. Among them, 12 toxicants were analytically confirmed using reference standards. Collectively, we present a "bioactivity-signature-toxicant" strategy to deconvolute mixtures and to connect patchy data sets and guide nontarget analysis for diverse chemicals that elicit the same bioactivity.

Recent Computer-Aided Studies on Herbicides: A Short Review

Current industrial herbicides have a negative impact on the environment and have widespread resistance, so computational studies on their properties, elimination, and overcoming resistance can be helpful. On the other hand, developing new herbicides, especially bioherbicides, is slow and costly. Therefore, computational studies that guide the design and search for new herbicides that exist in various plant sources, can alleviate the pain associated with the many obstacles. This review summarizes for the first time the most recent studies on both aspects of herbicides over 10 years.

Informing the public about chemical mixtures in the local environment: Currently applied indicators in the Netherlands and ways forward

The current use of chemicals puts pressure on human and ecological health. Based on the Aarhus Convention, citizens have the right to have access to information on substances in their local environment. Providing this information is a major challenge, especially considering complex mixtures, as the current substance-by-substance risk assessment may not adequately address the risk of co-exposure to multiple substances. Here, we provide an overview of the currently available indicators in the Netherlands to explore current scientific possibilities to indicate the impacts of complex chemical mixtures in the environment on human health and ecology at the local scale. This is limited to impact estimates on freshwater species for 701 substances, impact estimates of four metals on soil organisms, and impacts on human health for particulate matter (PM10) and nitrogen dioxide (NO2) in air. The main limiting factors in developing and expanding these indicators to cover more compartments and substances are the availability of emission and concentration data of substances and dose-response relationships at the population (human health) or community (ecology) level. As ways forward, we propose; 1) developing cumulative assessment groups (CAGs) for substances on the European Pollutant Transfer and Release Register and Water Framework Directive substance lists, to enable the development of mixture indicators based on mixture risk assessment and concentration addition principles; 2) to gain insight into local mixtures by also applying these CAGs to emission data, which is available for soil and air for more substances than concentrations data; 3) the application of analytical non-target screening methods as well as effect-based methods for whole-mixture assessment.

Machine learning methods to predict cadmium (Cd) concentration in rice grain and support soil management at a regional scale

Rice is a major dietary source of the toxic metal cadmium (Cd). Concentration of Cd in rice grain varies widely at the regional scale, and it is challenging to predict grain Cd concentration using soil properties. The lack of reliable predictive models hampers management of contaminated soils. Here, we conducted a three-year survey of 601 pairs of soil and rice samples at a regional scale. Approximately 78.3% of the soil samples exceeded the soil screening values for Cd in China, and 53.9% of rice grain samples exceeded the Chinese maximum permissible limit for Cd. Predictive models were developed using multiple linear regression and machine learning methods. The correlations between rice grain Cd and soil total Cd concentrations were poor (R 2 < 0.17). Both linear regression and machine learning methods identified four key factors that significantly affect grain Cd concentrations, including Fe-Mn oxide bound Cd, soil pH, field soil moisture content, and the concentration of soil reducible Mn. The machine learning-based support vector machine model showed the best performance (R 2 = 0.87) in predicting grain Cd concentrations at a regional scale, followed by machine learning-based random forest model (R 2 = 0.67), and back propagation neural network model (R 2 = 0.64). Scenario simulations revealed that liming soil to a target pH of 6.5 could be one of the most cost-effective approaches to reduce the exceedance of Cd in rice grain. Taken together, these results show that machine learning methods can be used to predict Cd concentration in rice grain reliably at a regional scale and to support soil management and safe rice production.

Predicting Odor Sensory Attributes of Unidentified Chemicals in Water Using Fragmentation Mass Spectra with Machine Learning Models

Knowing odor sensory attributes of odorants lies at the core of odor tracking when addressing waterborne odor issues. However, experimental determination covering tens of thousands of odorants in authentic water is not pragmatic due to the complexity of odorant identification and odor evaluation. In this study, we propose the first machine learning (ML) model to predict odor perception/threshold aiming at odorants in water, which can use either molecular structure or MS2 spectra as input features. We demonstrate that model performance using MS2 spectra is nearly as good as that using unequivocal structures, both with outstanding accuracy. We particularly show the model's robustness in predicting odor sensory attributes of unidentified chemicals by using the experimentally obtained MS2 spectra from nontarget analysis on authentic water samples. Interpreting the developed models, we identify the intricate interaction of functional groups as the predominant influence factor on odor sensory attributes. We also highlight the important roles of carbon chain length, molecular weight, etc., in the inherent olfactory mechanisms. These findings streamline the odor sensory attribute prediction and are crucial advancements toward credible tracking and efficient control of off-odors in water.

Grouping strategies for assessing and managing persistent and mobile substances

Background

Persistent, mobile and toxic (PMT), or very persistent and very mobile (vPvM) substances are a wide class of chemicals that are recalcitrant to degradation, easily transported, and potentially harmful to humans and the environment. Due to their persistence and mobility, these substances are often widespread in the environment once emitted, particularly in water resources, causing increased challenges during water treatment processes. Some PMT/vPvM substances such as GenX and perfluorobutane sulfonic acid have been identified as substances of very high concern (SVHCs) under the European Registration, Evaluation, Authorisation and Restriction of Chemicals (REACH) regulation. With hundreds to thousands of potential PMT/vPvM substances yet to be assessed and managed, effective and efficient approaches that avoid a case-by-case assessment and prevent regrettable substitution are necessary to achieve the European Union's zero-pollution goal for a non-toxic environment by 2050.

Main

Substance grouping has helped global regulation of some highly hazardous chemicals, e.g., through the Montreal Protocol and the Stockholm Convention. This article explores the potential of grouping strategies for identifying, assessing and managing PMT/vPvM substances. The aim is to facilitate early identification of lesser-known or new substances that potentially meet PMT/vPvM criteria, prompt additional testing, avoid regrettable use or substitution, and integrate into existing risk management strategies. Thus, this article provides an overview of PMT/vPvM substances and reviews the definition of PMT/vPvM criteria and various lists of PMT/vPvM substances available. It covers the current definition of groups, compares the use of substance grouping for hazard assessment and regulation, and discusses the advantages and disadvantages of grouping substances for regulation. The article then explores strategies for grouping PMT/vPvM substances, including read-across, structural similarity and commonly retained moieties, as well as the potential application of these strategies using cheminformatics to predict P, M and T properties for selected examples.

Conclusions

Effective substance grouping can accelerate the assessment and management of PMT/vPvM substances, especially for substances that lack information. Advances to read-across methods and cheminformatics tools are needed to support efficient and effective chemical management, preventing broad entry of hazardous chemicals into the global market and favouring safer and more sustainable alternatives.

Predicting the Activity of Unidentified Chemicals in Complementary Bioassays from the HRMS Data to Pinpoint Potential Endocrine Disruptors

The majority of chemicals detected via nontarget liquid chromatography high-resolution mass spectrometry (HRMS) in environmental samples remain unidentified, challenging the capability of existing machine learning models to pinpoint potential endocrine disruptors (EDs). Here, we predict the activity of unidentified chemicals across 12 bioassays related to EDs within the Tox21 10K dataset. Single- and multi-output models, utilizing various machine learning algorithms and molecular fingerprint features as an input, were trained for this purpose. To evaluate the models under near real-world conditions, Monte Carlo sampling was implemented for the first time. This technique enables the use of probabilistic fingerprint features derived from the experimental HRMS data with SIRIUS+CSI:FingerID as an input for models trained on true binary fingerprint features. Depending on the bioassay, the lowest false-positive rate at 90% recall ranged from 0.251 (sr.mmp, mitochondrial membrane potential) to 0.824 (nr.ar, androgen receptor), which is consistent with the trends observed in the models' performances submitted for the Tox21 Data Challenge. These findings underscore the informativeness of fingerprint features that can be compiled from HRMS in predicting the endocrine-disrupting activity. Moreover, an in-depth SHapley Additive exPlanations analysis unveiled the models' ability to pinpoint structural patterns linked to the modes of action of active chemicals. Despite the superior performance of the single-output models compared to that of the multi-output models, the latter's potential cannot be disregarded for similar tasks in the field of in silico toxicology. This study presents a significant advancement in identifying potentially toxic chemicals within complex mixtures without unambiguous identification and effectively reducing the workload for postprocessing by up to 75% in nontarget HRMS.

Innovative analytical methodologies for characterizing chemical exposure with a view to next-generation risk assessment

The chemical burden on the environment and human population is increasing. Consequently, regulatory risk assessment must keep pace to manage, reduce, and prevent adverse impacts on human and environmental health associated with hazardous chemicals. Surveillance of chemicals of known, emerging, or potential future concern, entering the environment-food-human continuum is needed to document the reality of risks posed by chemicals on ecosystem and human health from a one health perspective, feed into early warning systems and support public policies for exposure mitigation provisions and safe and sustainable by design strategies. The use of less-conventional sampling strategies and integration of full-scan, high-resolution mass spectrometry and effect-directed analysis in environmental and human monitoring programmes have the potential to enhance the screening and identification of a wider range of chemicals of known, emerging or potential future concern. Here, we outline the key needs and recommendations identified within the European Partnership for Assessment of Risks from Chemicals (PARC) project for leveraging these innovative methodologies to support the development of next-generation chemical risk assessment.

Online and Offline Prioritization of Chemicals of Interest in Suspect Screening and Non-targeted Screening with High-Resolution Mass Spectrometry

Recent advances in high-resolution mass spectrometry (HRMS) have enabled the detection of thousands of chemicals from a single sample, while computational methods have improved the identification and quantification of these chemicals in the absence of reference standards typically required in targeted analysis. However, to determine the presence of chemicals of interest that may pose an overall impact on ecological and human health, prioritization strategies must be used to effectively and efficiently highlight chemicals for further investigation. Prioritization can be based on a chemical's physicochemical properties, structure, exposure, and toxicity, in addition to its regulatory status. This Perspective aims to provide a framework for the strategies used for chemical prioritization that can be implemented to facilitate high-quality research and communication of results. These strategies are categorized as either "online" or "offline" prioritization techniques. Online prioritization techniques trigger the isolation and fragmentation of ions from the low-energy mass spectra in real time, with user-defined parameters. Offline prioritization techniques, in contrast, highlight chemicals of interest after the data has been acquired; detected features can be filtered and ranked based on the relative abundance or the predicted structure, toxicity, and concentration imputed from the tandem mass spectrum (MS2). Here we provide an overview of these prioritization techniques and how they have been successfully implemented and reported in the literature to find chemicals of elevated risk to human and ecological environments. A complete list of software and tools is available from https://nontargetedanalysis.org/.

Profiling emerging micropollutants in urban stormwater runoff using suspect and non-target screening via high-resolution mass spectrometry

Urban surface runoff contains chemicals that can negatively affect water quality. Urban runoff studies have determined the transport dynamics of many legacy pollutants. However, less attention has been paid to determining the first-flush effects (FFE) of emerging micropollutants using suspect and non-target screening (SNTS). Therefore, this study employed suspect and non-target analyses using liquid chromatography-high resolution mass spectrometry to detect emerging pollutants in urban receiving waters during stormwater events. Time-interval sampling was used to determine occurrence trends during stormwater events. Suspect screening tentatively identified 65 substances, then, their occurrence trend was grouped using correlation analysis. Non-target peaks were prioritized through hierarchical cluster analysis, focusing on the first flush-concentrated peaks. This approach revealed 38 substances using in silico identification. Simultaneously, substances identified through homologous series observation were evaluated for their observed trends in individual events using network analysis. The results of SNTS were normalized through internal standards to assess the FFE, and the most of tentatively identified substances showed observed FFE. Our findings suggested that diverse pollutants that could not be covered by target screening alone entered urban water through stormwater runoff during the first flush. This study showcases the applicability of the SNTS in evaluating the FFE of urban pollutants, offering insights for first-flush stormwater monitoring and management.

Machine Learning-Based Hazard-Driven Prioritization of Features in Nontarget Screening of Environmental High-Resolution Mass Spectrometry Data

Nontarget high-resolution mass spectrometry screening (NTS HRMS/MS) can detect thousands of organic substances in environmental samples. However, new strategies are needed to focus time-intensive identification efforts on features with the highest potential to cause adverse effects instead of the most abundant ones. To address this challenge, we developed MLinvitroTox, a machine learning framework that uses molecular fingerprints derived from fragmentation spectra (MS2) for a rapid classification of thousands of unidentified HRMS/MS features as toxic/nontoxic based on nearly 400 target-specific and over 100 cytotoxic endpoints from ToxCast/Tox21. Model development results demonstrated that using customized molecular fingerprints and models, over a quarter of toxic endpoints and the majority of the associated mechanistic targets could be accurately predicted with sensitivities exceeding 0.95. Notably, SIRIUS molecular fingerprints and xboost (Extreme Gradient Boosting) models with SMOTE (Synthetic Minority Oversampling Technique) for handling data imbalance were a universally successful and robust modeling configuration. Validation of MLinvitroTox on MassBank spectra showed that toxicity could be predicted from molecular fingerprints derived from MS2 with an average balanced accuracy of 0.75. By applying MLinvitroTox to environmental HRMS/MS data, we confirmed the experimental results obtained with target analysis and narrowed the analytical focus from tens of thousands of detected signals to 783 features linked to potential toxicity, including 109 spectral matches and 30 compounds with confirmed toxic activity.

Prioritizing molecular formulae identified by non-target analysis through high-throughput modelling: application to identify compounds with high human accumulation potential from house dust

Because it is typically not possible to pursue compound identification efforts for all chemical features detected during non-target analysis (NTA), the need for prioritization arises. Here we propose a strategy that ranks chemical features detected in environmental samples based on a model-derived metric that quantifies a feature's attribute that makes it desirable to elucidate its structure, e.g., a high potential for bioaccumulation in humans or wildlife. The procedure involves the identification of isomers that could plausibly represent the molecular formulae assigned to NTA-detected chemical features. For each isomer, the prioritization metric is calculated using properties predicted with high-throughput methods. After the molecular formulae are ranked based on the average values of the prioritization metric calculated for all isomers assigned to a formula, the highest ranked molecular formulae are prioritized for structure elucidation. We applied this workflow to features identified in house dust, using the ratio of chemical intake through dust ingestion to chemical concentration in blood (dose-to-concentration ratio, DCR) as the prioritization metric. Collections of isomers for the molecular formulae were assembled from the PubChem database and DCR was estimated using partitioning and biotransformation properties predicted for each isomer using quantitative structure property relationships. The ten top-ranked molecular formulae with notably lower average DCR-values represented mostly compounds already known to be indoor pollutants of concern, such as two polybrominated diphenyl ethers, bis(2-ethylhexyl) tetrabromophthalate, tetrabromobisphenol A, tris(1,3-dichloroisopropyl)phosphate and the azo dye disperse blue 373.

Predicting the Diagnostic Information of Tandem Mass Spectra of Environmentally Relevant Compounds Using Machine Learning

Acquisition and processing of informative tandem mass spectra (MS2) is crucial for numerous applications, including library-based (tentative) identification, feature prioritization, and prediction of chemical and toxicological characteristics. However, for environmentally relevant compounds, approaches to automatically assess the quality of the MS2 spectra are missing. This work focused on developing a machine learning-based approach to automatically evaluate the diagnostic information of MS2 spectra (e.g., number, distribution, and intensity of diagnostic fragments) of environmentally relevant compounds analyzed with electrospray ionization. For this, approximately 1400 MS2 spectra of 204 environmental contaminants, acquired with different collision energies using liquid chromatography coupled to high-resolution mass spectrometry, were used to train a random forest classifier to distinguish between spectra providing good or poor diagnostic information. Prior to training, validation, and testing, spectra were manually labeled based on criteria such as number, intensity, range of fragments present, molecular ion intensity, and noise levels. Subsequently, feature engineering and selection were applied to retrieve relevant variables from raw MS2 spectra as inputs for the classifier. The optimal set of features based on model performances was selected and used to train a final model, which showed an accuracy of 84%, a precision of 88%, and a recall of 75%. Results show that the combination of selected features and the machine learning model used here can effectively distinguish between MS2 spectra providing good or poor diagnostic information according to the defined criteria. The developed model has the potential to improve a broad range of applications that rely on MS2 data.

Complementary methods for structural assignment of isomeric candidate structures in non-target liquid chromatography ion mobility high-resolution mass spectrometric analysis

Non-target screening with LC/IMS/HRMS is increasingly employed for detecting and identifying the structure of potentially hazardous chemicals in the environment and food. Structural assignment relies on a combination of multidimensional instrumental methods and computational methods. The candidate structures are often isomeric, and unfortunately, assigning the correct structure among a number of isomeric candidate structures still is a key challenge both instrumentally and computationally. While practicing non-target screening, it is usually impossible to evaluate separately the limitations arising from (1) the inability of LC/IMS/HRMS to resolve the isomeric candidate structures and (2) the uncertainty of in silico methods in predicting the analytical information of isomeric candidate structures due to the lack of analytical standards for all candidate structures. Here we evaluate the feasibility of structural assignment of isomeric candidate structures based on in silico-predicted retention time and database collision cross-section (CCS) values as well as based on matching the empirical analytical properties of the detected feature with those of the analytical standards. For this, we investigated 14 candidate structures corresponding to five features detected with LC/HRMS in a spiked surface water sample. Considering the predicted retention times and database CCS values with the accompanying uncertainty, only one of the isomeric candidate structures could be deemed as unlikely; therefore, the annotation of the LC/IMS/HRMS features remained ambiguous. To further investigate if unequivocal annotation is possible via analytical standards, the reversed-phase LC retention times and low- and high-resolution ion mobility spectrometry separation, as well as high-resolution MS2 spectra of analytical standards were studied. Reversed-phase LC separated the highest number of candidate structures while low-resolution ion mobility and high-resolution MS2 spectra provided little means for pinpointing the correct structure among the isomeric candidate structures even if analytical standards were available for comparison. Furthermore, the question arises which prediction accuracy is required from the in silico methods to par the analytical separation. Based on the experimental data of the isomeric candidate structures studied here and previously published in the literature (516 retention time and 569 CCS values), we estimate that to reduce the candidate list by 95% of the structures, the confidence interval of the predicted retention times would need to decrease to below 0.05 min for a 15-min gradient while that of CCS values would need to decrease to 0.15%. Hereby, we set a clear goal to the in silico methods for retention time and CCS prediction.

Bypassing the Identification: MS2Quant for Concentration Estimations of Chemicals Detected with Nontarget LC-HRMS from MS2 Data

Nontarget analysis by liquid chromatography-high-resolution mass spectrometry (LC-HRMS) is now widely used to detect pollutants in the environment. Shifting away from targeted methods has led to detection of previously unseen chemicals, and assessing the risk posed by these newly detected chemicals is an important challenge. Assessing exposure and toxicity of chemicals detected with nontarget HRMS is highly dependent on the knowledge of the structure of the chemical. However, the majority of features detected in nontarget screening remain unidentified and therefore the risk assessment with conventional tools is hampered. Here, we developed MS2Quant, a machine learning model that enables prediction of concentration from fragmentation (MS2) spectra of detected, but unidentified chemicals. MS2Quant is an xgbTree algorithm-based regression model developed using ionization efficiency data for 1191 unique chemicals that spans 8 orders of magnitude. The ionization efficiency values are predicted from structural fingerprints that can be computed from the SMILES notation of the identified chemicals or from MS2 spectra of unidentified chemicals using SIRIUS+CSI:FingerID software. The root mean square errors of the training and test sets were 0.55 (3.5×) and 0.80 (6.3×) log-units, respectively. In comparison, ionization efficiency prediction approaches that depend on assigning an unequivocal structure typically yield errors from 2× to 6×. The MS2Quant quantification model was validated on a set of 39 environmental pollutants and resulted in a mean prediction error of 7.4×, a geometric mean of 4.5×, and a median of 4.0×. For comparison, a model based on PaDEL descriptors that depends on unequivocal structural assignment was developed using the same dataset. The latter approach yielded a comparable mean prediction error of 9.5×, a geometric mean of 5.6×, and a median of 5.2× on the validation set chemicals when the top structural assignment was used as input. This confirms that MS2Quant enables to extract exposure information for unidentified chemicals which, although detected, have thus far been disregarded due to lack of accurate tools for quantification. The MS2Quant model is available as an R-package in GitHub for improving discovery and monitoring of potentially hazardous environmental pollutants with nontarget screening.

Identifying potential toxic organic substances in leachates from tire wear particles and their mechanisms of toxicity to Scenedesmus obliquus

Tire wear particles (TWPs) are increasingly being found in the aquatic environment. However, there is limited information available on the environmental consequences of TWP constituents that may be release into water. In this study, TWP leachate samples were obtained by immersing TWPs in ultrapure water. Using high-resolution mass spectrometry and toxicity identification, we identified potentially toxic organic substances in the TWP leachates. Additionally, we investigated their toxicity and underlying mechanisms. Through our established workflow, we structurally identified 13 substances using reference standards. The median effective concentration (EC₅₀) of TWP leachates on Scenedesmus obliquus growth was comparable to that of simulated TWP leachates prepared with consistent concentrations of the 13 identified substances, indicating their dominance in the toxicity of TWP leachates. Among these substances, cyclic amines (EC₅₀: 1.04-3.65 mg/L) were found to be toxic to S. obliquus. We observed significant differential metabolites in TWP leachate-exposed S. obliquus, primarily associated with linoleic acid metabolism and purine metabolism. Oxidative stress was identified as a crucial factor in algal growth inhibition. Our findings shed light on the risk posed by TWP leachable substances to aquatic organisms.

Machine Learning and Nontargeted Liquid Chromatography–Mass Spectrometry to Assess Ecotoxicity

Anneli Kruve and Pilleriin Peets from Stockholm University in Sweden, discuss their latest research in machine learning and nontargeted liquid chromatography–mass spectrometry (LC–MS) to assess ecotoxicity.