In silico approaches for the prediction of the breakthrough of organic contaminants in wastewater treatment plants.

The removal efficiency (RE) of organic contaminants in wastewater treatment plants (WWTPs) is a major determinant of the environmental impact of chemicals which are discharged to wastewater. In a recent study, non-target screening analysis was applied to quantify the percentage removal efficiency (RE%) of more than 300 polar contaminants, by analyzing influent and effluent samples from a Swedish WWTP with direct injection UHPLC-Orbitrap-MS/MS. Based on subsets extracted from these data, we developed quantitative structure-property relationships (QSPRs) for the prediction of WWTP breakthrough (BT) to the effluent water. QSPRs were developed by means of multiple linear regression (MLR) and were selected after checking for overfitting and chance relationships by means of bootstrap and randomization procedures. A first model provided good fitting performance, showing that the proposed approach for the development of QSPRs for the prediction of BT is reasonable. By further populating the dataset with similar chemicals using a Tanimoto index approach based on substructure count fingerprints, a second QSPR indicated that the prediction of BT is also applicable to new chemicals sufficiently similar to the training set. Finally, a class-specific QSPR for PEGs and PPGs showed BT prediction trends consistent with known degradation pathways.

QSAR Models for the Prediction of Dietary Biomagnification Factor in Fish.

Xenobiotics released in the environment can be taken up by aquatic and terrestrial organisms and can accumulate at higher concentrations through the trophic chain. Bioaccumulation is therefore one of the PBT properties that authorities require to assess for the evaluation of the risks that chemicals may pose to humans and the environment. The use of an integrated testing strategy (ITS) and the use of multiple sources of information are strongly encouraged by authorities in order to maximize the information available and reduce testing costs. Moreover, considering the increasing demand for development and the application of new approaches and alternatives to animal testing, the development of in silico cost-effective tools such as QSAR models becomes increasingly important. In this study, a large and curated literature database of fish laboratory-based values of dietary biomagnification factor (BMF) was used to create externally validated QSARs. The quality categories (high, medium, low) available in the database were used to extract reliable data to train and validate the models, and to further address the uncertainty in low-quality data. This procedure was useful for highlighting problematic compounds for which additional experimental effort would be required, such as siloxanes, highly brominated and chlorinated compounds. Two models were suggested as final outputs in this study, one based on good-quality data and the other developed on a larger dataset of consistent Log BMFL values, which included lower-quality data. The models had similar predictive ability; however, the second model had a larger applicability domain. These QSARs were based on simple MLR equations that could easily be applied for the predictions of dietary BMFL in fish, and support bioaccumulation assessment procedures at the regulatory level. To ease the application and dissemination of these QSARs, they were included with technical documentation (as QMRF Reports) in the QSAR-ME Profiler software for QSAR predictions available online.

Aquatic toxicity prediction of diverse pesticides on two algal species using QSTR modeling approach.

With the aim of identification of toxic nature of the diverse pesticides on the aquatic compartment, a large dataset of pesticides (n = 325) with experimental toxicity data on two algal test species (Pseudokirchneriella subcapitata (PS) (synonym: Raphidocelis subcapitata, Selenastrum capricornutum) and Scenedemus subspicatus (SS)) was gathered and subjected to quantitative structure toxicity relationship (QSTR) analysis to predict aquatic toxicity of pesticides. The QSTR models were developed by multiple linear regressions (MLRs), and the genetic algorithm (GA) was used for the variable selection. The developed GA-MLR models were statistically robust enough internally (Q2LOO = 0.620-0.663) and externally (Q2Fn = 0.693-0.868, CCCext = 0.843-0.877). The leverage approach of applicability domain (AD) and prediction reliability indicator assured the reliability of the developed models. The mechanistic interpretation highlighted that the presence of SO2, F and aromatic rings influenced the toxicity of pesticides towards PS species while the presence of alkyl, alkyl halide, aromatic rings and carbonyl was responsible for the toxicity of pesticides towards SS species. Additionally, we have reported the application of developed models to pesticides without experimental value and the cumulative toxicity of pesticides on the aquatic environment by using principal component analysis (PCA). The reliable prediction and prioritization of toxic compounds from the developed models will be useful in the aquatic toxicity assessment of pesticides.

Predicting the Bioconcentration Factor in Fish from Molecular Structures.

The bioconcentration factor (BCF) is one of the metrics used to evaluate the potential of a substance to bioaccumulate into aquatic organisms. In this work, linear and non-linear regression QSARs were developed for the prediction of log BCF using different computational approaches, and starting from a large and structurally heterogeneous dataset. The new MLR-OLS and ANN regression models have good fitting with R2 values of 0.62 and 0.70, respectively, and comparable external predictivity with R2ext 0.64 and 0.65 (RMSEext of 0.78 and 0.76), respectively. Furthermore, linear and non-linear classification models were developed using the regulatory threshold BCF >2000. A class balanced subset was used to develop classification models which were applied to chemicals not used to create the QSARs. These classification models are characterized by external and internal accuracy up to 84% and 90%, respectively, and sensitivity and specificity up to 90% and 80%, respectively. QSARs presented in this work are validated according to regulatory requirements and their quality is in line with other tools available for the same endpoint and dataset, with the advantage of low complexity and easy application through the software QSAR-ME Profiler. These QSARs can be used as alternatives for, or in combination with, existing models to support bioaccumulation assessment procedures.

QSARINS-Chem standalone version: A new platform-independent software to profile chemicals for physico-chemical properties, fate, and toxicity.

The new software QSARINS-Chem standalone version is a multiplatform tool, freely downloadable, for the in silico profiling of multiple properties and activities of organic chemicals. This software, which is based on the concept of the QSARINS-chem module embedded in the QSARINS software, has been fully redesigned and redeveloped in the Java™ language. In addition to a selection of models included in the old module, the new software predicts biotransformation rates and aquatic toxicities of pharmaceuticals and personal care products in multiple organisms, and offers a suite of tools for the analysis of predictions. Furthermore, a comprehensive and transparent database of molecular structures is provided. The new QSARINS-Chem standalone version is an informative and solid tool, which is useful to support the assessment of the potential hazard and risks related to organic chemicals and is dedicated to users which are interested in the application of QSARs to generate reliable predictions.

New BODIPYs for photodynamic therapy (PDT): Synthesis and activity on human cancer cell lines.

A new class of compounds based on the 4,4-difluoro-4-bora-3a,4a-diaza-s-indacene core, known as BODIPYs, has attracted significant attention as photosensitizers suitable for application in photodynamic therapy (PDT), which is a minimally invasive procedure to treat cancer. In PDT the combination of a photosensitizer (PS), light, and oxygen leads to a series of photochemical reactions generating reactive oxygen species (ROS) exerting cytotoxic action on tumor cells. Here we present the synthesis and the study of the in vitro photodynamic effects of two BODIPYs which differ in the structure of the substituent placed on the meso (or 8) position of the dipyrrolylmethenic nucleus. The two compounds were tested on three human cancer cell lines of different origin and degree of malignancy. Our results indicate that the BODIPYs are very effective in reducing the growth/viability of HCT116, SKOV3 and MCF7 cells when irradiated with a green LED source, whereas they are practically devoid of activity in the dark. Phototoxicity occurs mainly through apoptotic cell death, however necrotic cell death also seems to play a role. Furthermore, singlet oxygen generation and induction of the increase of reactive oxygen species also appear to be involved in the photodynamic effect of the BODIPYs. Finally, it is worth noting that the two BODIPYs are also able to exert anti-migratory activity.

Are In Silico Approaches Applicable As a First Step for the Prediction of e-Liquid Toxicity in e-Cigarettes?

Recent studies have raised concerns about e-cigarette liquid inhalation toxicity by reporting the presence of chemicals with European Union CLP toxicity classification. In this scenario, the regulatory context is still developing and is not yet up to date with vaping current reality. Due to the paucity of toxicological studies, robust data regarding which components in e-liquids exhibit potential toxicities, are still inconsistent. In this study we applied computational methods for estimating the toxicity of poorly studied chemicals as a useful tool for predicting the acute toxicity of chemicals contained in e-liquids. The purpose of this study was 3-fold: (a) to provide a lower tier assessment of the potential health concerns associated with e-liquid ingredients, (b) to prioritize e-liquid ingredients by calculating the e-tox index, and (c) to estimate acute toxicity of e-liquid mixtures. QSAR models were generated using QSARINS software to fill the acute toxicity data gap of 264 e-liquid ingredients. As a second step, the potential acute toxicity of e-liquids mixtures was evaluated. Our preliminary data suggest that a computational approach may serve as a roadmap to enable regulatory bodies to better regulate e-liquid composition and to contribute to consumer health protection.

Application of chemometric methods and QSAR models to support pesticide risk assessment starting from ecotoxicological datasets.

The EFSA 'Guidance on tiered risk assessment for edge-of-field surface waters' underscores the importance of in silico models to support the pesticide risk assessment. The aim of this work was to use in silico models starting from an available, structured and harmonized pesticide dataset that was developed for different purposes, in order to stimulate the use of QSAR models for risk assessment. The present work focuses on the development of a set of in silico models, developed to predict the aquatic toxicity of heterogeneous pesticides with incomplete/unknown toxic behavior in the water compartment. The generated models have good fitting performances (R2: 0.75-0.99), they are internally robust (Q2loo: 0.66-0.98) and can handle up to 30% of perturbation of the training set (Q2 lmo: 0.64-0.98). The absence of chance correlation was guaranteed by low values of R2 calculated on scrambled responses (R2 Yscr: 0.11-0.38). Different statistical parameters were used to quantify the external predictivity of the models (CCCext: 0.73-0.91, Q2 ext-Fn: 0.53-0.96). The results indicate that all the best models are predictive when applied to chemicals not involved in the models development. In addition, all models have similar accuracy both in fitting and in prediction and this represents a good degree of generalization. These models may be useful to support the risk assessment procedure when experimental data for key species are missing or to create prioritization lists for the general a priori assessment of the potential toxicity of existing and new pesticides which fall in the applicability domain.

CoMPARA: Collaborative Modeling Project for Androgen Receptor Activity.

BACKGROUND: Endocrine disrupting chemicals (EDCs) are xenobiotics that mimic the interaction of natural hormones and alter synthesis, transport, or metabolic pathways. The prospect of EDCs causing adverse health effects in humans and wildlife has led to the development of scientific and regulatory approaches for evaluating bioactivity. This need is being addressed using high-throughput screening (HTS) in vitro approaches and computational modeling. OBJECTIVES: In support of the Endocrine Disruptor Screening Program, the U.S. Environmental Protection Agency (EPA) led two worldwide consortiums to virtually screen chemicals for their potential estrogenic and androgenic activities. Here, we describe the Collaborative Modeling Project for Androgen Receptor Activity (CoMPARA) efforts, which follows the steps of the Collaborative Estrogen Receptor Activity Prediction Project (CERAPP). METHODS: The CoMPARA list of screened chemicals built on CERAPP's list of 32,464 chemicals to include additional chemicals of interest, as well as simulated ToxCast™ metabolites, totaling 55,450 chemical structures. Computational toxicology scientists from 25 international groups contributed 91 predictive models for binding, agonist, and antagonist activity predictions. Models were underpinned by a common training set of 1,746 chemicals compiled from a combined data set of 11 ToxCast™/Tox21 HTS in vitro assays. RESULTS: The resulting models were evaluated using curated literature data extracted from different sources. To overcome the limitations of single-model approaches, CoMPARA predictions were combined into consensus models that provided averaged predictive accuracy of approximately 80% for the evaluation set. DISCUSSION: The strengths and limitations of the consensus predictions were discussed with example chemicals; then, the models were implemented into the free and open-source OPERA application to enable screening of new chemicals with a defined applicability domain and accuracy assessment. This implementation was used to screen the entire EPA DSSTox database of ∼875,000 chemicals, and their predicted AR activities have been made available on the EPA CompTox Chemicals dashboard and National Toxicology Program's Integrated Chemical Environment. https://doi.org/10.1289/EHP5580.

High-throughput evaluation of organic contaminant removal efficiency in a wastewater treatment plant using direct injection UHPLC-Orbitrap-MS/MS.

The removal efficiency (RE) of organic contaminants in wastewater treatment plants (WWTPs) is a major determinant of the environmental impact of these contaminants. However, RE data are available for only a few chemicals due to the time and cost required for conventional target analysis. In the present study, we applied non-target screening analysis to evaluate the RE of polar contaminants, by analyzing influent and effluent samples from a Swedish WWTP with direct injection UHPLC-Orbitrap-MS/MS. Matrix effects were evaluated by spiking the samples with isotope-labeled standards of 40 polar contaminants. For 85% of the compounds, the matrix effects in the influent and effluent were not significantly different. Approximately 10 000 compounds were detected in the wastewater, of which 319 were identified by using the online database mzCloud. Level 1 identification confidence was achieved for 31 compounds for which we had reference standards, and level 2 was achieved for the remainder. RE was calculated from the ratio of the peak areas in the influent and the effluent from the non-target analysis. Good agreement was found with RE determined from the target analysis of the target compounds. The method generated reliable estimates of RE for large numbers of contaminants with comparatively low effort and is foreseen to be particularly useful in applications where information on a large number of chemicals is needed.

Development of human biotransformation QSARs and application for PBT assessment refinement.

Toxicokinetics heavily influence chemical toxicity as the result of Absorption, Distribution, Metabolism (Biotransformation) and Elimination (ADME) processes. Biotransformation (metabolism) reactions can lead to detoxification or, in some cases, bioactivation of parent compounds to more toxic chemicals. Moreover, biotransformation has been recognized as a key process determining chemical half-life in an organism and is thus a key determinant for bioaccumulation assessment for many chemicals. This study addresses the development of QSAR models for the prediction of in vivo whole body human biotransformation (metabolism) half-lives measured or empirically-derived for over 1000 chemicals, mainly represented by pharmaceuticals. Models presented in this study meet regulatory standards for fitting, validation and applicability domain. These QSARs were used, in combination with literature models for the prediction of biotransformation half-lives in fish, to refine the screening of the potential PBT behaviour of over 1300 Pharmaceuticals and Personal Care Products (PPCPs). The refinement of the PBT screening allowed, among others, for the identification of PPCPs, which were predicted as PBTs on the basis of their chemical structure, but may be easily biotransformed. These compounds are of lower concern in comparison to potential PBTs characterized by large predicted biotransformation half-lives.

QSAR modeling of cumulative environmental end-points for the prioritization of hazardous chemicals.

The hazard of chemicals in the environment is inherently related to the molecular structure and derives simultaneously from various chemical properties/activities/reactivities. Models based on Quantitative Structure Activity Relationships (QSARs) are useful to screen, rank and prioritize chemicals that may have an adverse impact on humans and the environment. This paper reviews a selection of QSAR models (based on theoretical molecular descriptors) developed for cumulative multivariate endpoints, which were derived by mathematical combination of multiple effects and properties. The cumulative end-points provide an integrated holistic point of view to address environmentally relevant properties of chemicals.

Modeling ready biodegradability of fragrance materials.

In the present study, quantitative structure activity relationships were developed for predicting ready biodegradability of approximately 200 heterogeneous fragrance materials. Two classification methods, classification and regression tree (CART) and k-nearest neighbors (kNN), were applied to perform the modeling. The models were validated with multiple external prediction sets, and the structural applicability domain was verified by the leverage approach. The best models had good sensitivity (internal ≥80%; external ≥68%), specificity (internal ≥80%; external 73%), and overall accuracy (≥75%). Results from the comparison with BIOWIN global models, based on group contribution method, show that specific models developed in the present study perform better in prediction than BIOWIN6, in particular for the correct classification of not readily biodegradable fragrance materials.

The QSPR-THESAURUS: the online platform of the CADASTER project.

The aim of the CADASTER project (CAse Studies on the Development and Application of in Silico Techniques for Environmental Hazard and Risk Assessment) was to exemplify REACH-related hazard assessments for four classes of chemical compound, namely, polybrominated diphenylethers, per and polyfluorinated compounds, (benzo)triazoles, and musks and fragrances. The QSPR-THESAURUS website (http: / /qspr-thesaurus.eu) was established as the project's online platform to upload, store, apply, and also create, models within the project. We overview the main features of the website, such as model upload, experimental design and hazard assessment to support risk assessment, and integration with other web tools, all of which are essential parts of the QSPR-THESAURUS.

Assessing predictive uncertainty in comparative toxicity potentials of triazoles.

Comparative toxicity potentials (CTPs) quantify the potential ecotoxicological impacts of chemicals per unit of emission. They are the product of a substance's environmental fate, exposure, and hazardous concentration. When empirical data are lacking, substance properties can be predicted. The goal of the present study was to assess the influence of predictive uncertainty in substance property predictions on the CTPs of triazoles. Physicochemical and toxic properties were predicted with quantitative structure-activity relationships (QSARs), and uncertainty in the predictions was quantified with use of the data underlying the QSARs. Degradation half-lives were based on a probability distribution representing experimental half-lives of triazoles. Uncertainty related to the species' sample size that was present in the prediction of the hazardous aquatic concentration was also included. All parameter uncertainties were treated as probability distributions, and propagated by Monte Carlo simulations. The 90% confidence interval of the CTPs typically spanned nearly 4 orders of magnitude. The CTP uncertainty was mainly determined by uncertainty in soil sorption and soil degradation rates, together with the small number of species sampled. In contrast, uncertainty in species-specific toxicity predictions contributed relatively little. The findings imply that the reliability of CTP predictions for the chemicals studied can be improved particularly by including experimental data for soil sorption and soil degradation, and by developing toxicity QSARs for more species.

Metabolic biotransformation half-lives in fish: QSAR modeling and consensus analysis.

Bioaccumulation in fish is a function of competing rates of chemical uptake and elimination. For hydrophobic organic chemicals bioconcentration, bioaccumulation and biomagnification potential are high and the biotransformation rate constant is a key parameter. Few measured biotransformation rate constant data are available compared to the number of chemicals that are being evaluated for bioaccumulation hazard and for exposure and risk assessment. Three new Quantitative Structure-Activity Relationships (QSARs) for predicting whole body biotransformation half-lives (HLN) in fish were developed and validated using theoretical molecular descriptors that seek to capture structural characteristics of the whole molecule and three data set splitting schemes. The new QSARs were developed using a minimal number of theoretical descriptors (n=9) and compared to existing QSARs developed using fragment contribution methods that include up to 59 descriptors. The predictive statistics of the models are similar thus further corroborating the predictive performance of the different QSARs; Q(2)ext ranges from 0.75 to 0.77, CCCext ranges from 0.86 to 0.87, RMSE in prediction ranges from 0.56 to 0.58. The new QSARs provide additional mechanistic insights into the biotransformation capacity of organic chemicals in fish by including whole molecule descriptors and they also include information on the domain of applicability for the chemical of interest. Advantages of consensus modeling for improving overall prediction and minimizing false negative errors in chemical screening assessments, for identifying potential sources of residual error in the empirical HLN database, and for identifying structural features that are not well represented in the HLN dataset to prioritize future testing needs are illustrated.

Daphnia and fish toxicity of (benzo)triazoles: validated QSAR models, and interspecies quantitative activity-activity modelling.

Due to their chemical properties synthetic triazoles and benzo-triazoles ((B)TAZs) are mainly distributed to the water compartments in the environment, and because of their wide use the potential effects on aquatic organisms are cause of concern. Non testing approaches like those based on quantitative structure-activity relationships (QSARs) are valuable tools to maximize the information contained in existing experimental data and predict missing information while minimizing animal testing. In the present study, externally validated QSAR models for the prediction of acute (B)TAZs toxicity in Daphnia magna and Oncorhynchus mykiss have been developed according to the principles for the validation of QSARs and their acceptability for regulatory purposes, proposed by the Organization for Economic Co-operation and Development (OECD). These models are based on theoretical molecular descriptors, and are statistically robust, externally predictive and characterized by a verifiable structural applicability domain. They have been applied to predict acute toxicity for over 300 (B)TAZs without experimental data, many of which are in the pre-registration list of the REACH regulation. Additionally, a model based on quantitative activity-activity relationships (QAAR) has been developed, which allows for interspecies extrapolation from daphnids to fish. The importance of QSAR/QAAR, especially when dealing with specific chemical classes like (B)TAZs, for screening and prioritization of pollutants under REACH, has been highlighted.