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.

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.

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.

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.

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.

Experimental assessment of the environmental fate and effects of triazoles and benzotriazole.

The environmental fate and effects of triazoles and benzotriazoles are of concern within the context of chemical regulation. As part of an intelligent testing strategy, experimental tests were performed on endpoints that are relevant for risk assessment. The experimental tests included the assessment of ecotoxicity to an alga, a daphnid and zebrafish embryos, and the assessment of ready biodegradability. Triazole and benzotriazole compounds were selected for testing, based on existing toxicity data for vertebrate and invertebrate species, as well as on the principal component analysis of molecular descriptors aimed at selecting the minimum number of test compounds in order to maximise the chemical domain spanned for both compound classes. The experimental results show that variation in the toxicities of triazoles and benzotriazole across species was relatively minor; in general, the largest factor was approximately 20. The study conducted indicated that triazoles are not readily biodegradable.

Evaluation of CADASTER QSAR models for the aquatic toxicity of (benzo)triazoles and prioritisation by consensus prediction.

QSAR regression models of the toxicity of triazoles and benzotriazoles ([B]TAZs) to an alga (Pseudokirchneriella subcapitata), Daphnia magna and a fish (Onchorhynchus mykiss), were developed by five partners in the FP7-EU Project, CADASTER. The models were developed by different methods - Ordinary Least Squares (OLS), Partial Least Squares (PLS), Bayesian regularised regression and Associative Neural Network (ASNN) - by using various molecular descriptors (DRAGON, PaDEL-Descriptor and QSPR-THESAURUS web). In addition, different procedures were used for variable selection, validation and applicability domain inspection. The predictions of the models developed, as well as those obtained in a consensus approach by averaging the data predicted from each model, were compared with the results of experimental tests that were performed by two CADASTER partners. The individual and consensus models were able to correctly predict the toxicity classes of the chemicals tested in the CADASTER project, confirming the utility of the QSAR approach. The models were also used for the prediction of aquatic toxicity of over 300 (B)TAZs, many of which are included in the REACH pre-registration list, and were without experimental data. This highlights the importance of QSAR models for the screening and prioritisation of untested chemicals, in order to reduce and focus experimental testing.

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.

Antiproliferative Pt(IV) complexes: synthesis, biological activity, and quantitative structure-activity relationship modeling.

Several Pt(IV) complexes of the general formula [Pt(L)2(L')2(L'')2] [axial ligands L are Cl-, RCOO-, or OH-; equatorial ligands L' are two am(m)ine or one diamine; and equatorial ligands L'' are Cl- or glycolato] were rationally designed and synthesized in the attempt to develop a predictive quantitative structure-activity relationship (QSAR) model. Numerous theoretical molecular descriptors were used alongside physicochemical data (i.e., reduction peak potential, Ep, and partition coefficient, log Po/w) to obtain a validated QSAR between in vitro cytotoxicity (half maximal inhibitory concentrations, IC50, on A2780 ovarian and HCT116 colon carcinoma cell lines) and some features of Pt(IV) complexes. In the resulting best models, a lipophilic descriptor (log Po/w or the number of secondary sp3 carbon atoms) plus an electronic descriptor (Ep, the number of oxygen atoms, or the topological polar surface area expressed as the N,O polar contribution) is necessary for modeling, supporting the general finding that the biological behavior of Pt(IV) complexes can be rationalized on the basis of their cellular uptake, the Pt(IV)-->Pt(II) reduction, and the structure of the corresponding Pt(II) metabolites. Novel compounds were synthesized on the basis of their predicted cytotoxicity in the preliminary QSAR model, and were experimentally tested. A final QSAR model, based solely on theoretical molecular descriptors to ensure its general applicability, is proposed.

QSAR modeling and prediction of the endocrine-disrupting potencies of brominated flame retardants.

In the European Union REACH regulation, the chemicals with particularly harmful behaviors, such as endocrine disruptors (EDs), are subject to authorization, and the identification of safer alternatives to these chemicals is required. In this context, the use of quantitative structure-activity relationships (QSAR) becomes particularly useful to fill the data gap due to the very small number of experimental data available to characterize the environmental and toxicological profiles of new and emerging pollutants with ED behavior such as brominated flame retardants (BFRs). In this study, different QSAR models were developed on different responses of endocrine disruption measured for several BFRs. The multiple linear regression approach was applied to a variety of theoretical molecular descriptors, and the best models, which were identified from all of the possible combinations of the structural variables, were internally validated for their performance using the leave-one-out (Q(LOO)(2) = 73-91%) procedure and scrambling of the responses. External validation was provided, when possible, by splitting the data sets in training and test sets (range of Q(EXT)(2) = 76-90%), which confirmed the predictive ability of the proposed equations. These models, which were developed according to the principles defined by the Organization for Economic Co-operation and Development to improve the regulatory acceptance of QSARs, represent a simple tool for the screening and characterization of BFRs.

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.

Evaluation and QSAR modeling on multiple endpoints of estrogen activity based on different bioassays.

There is a great need for an effective means of rapidly assessing endocrine-disrupting activity, especially estrogen-simulating activity, due to the large number of chemicals that have serious adverse effects on the environment. Many approaches using a variety of biological screening assays are used to identify endocrine disrupting chemicals. The present investigation analyzes the consistency and peculiarity of information from different experimental assays collected from a literature survey, by studying the correlation of the different endpoints. In addition, the activity values of more widely used selected bioassays have been combined by principle components analysis (PCA) to build one cumulative endpoint, the estrogen activity index (EAI), for priority setting to identify chemicals most likely possessing estrogen activity for early entry into screening. This index was then modeled using only a few theoretical molecular descriptors. The constructed MLR-QSAR model has been statistically validated for its predictive power, and can be proposed as a preliminary evaluative method to screen/prioritize estrogens according to their integrated estrogen activity, just starting from molecular structure.

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.

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.

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.

Statistical external validation and consensus modeling: a QSPR case study for Koc prediction.

The soil sorption partition coefficient (log K(oc)) of a heterogeneous set of 643 organic non-ionic compounds, with a range of more than 6 log units, is predicted by a statistically validated QSAR modeling approach. The applied multiple linear regression (ordinary least squares, OLS) is based on a variety of theoretical molecular descriptors selected by the genetic algorithms-variable subset selection (GA-VSS) procedure. The models were validated for predictivity by different internal and external validation approaches. For external validation we applied self organizing maps (SOM) to split the original data set: the best four-dimensional model, developed on a reduced training set of 93 chemicals, has a predictivity of 78% when applied on 550 validation chemicals (prediction set). The selected molecular descriptors, which could be interpreted through their mechanistic meaning, were compared with the more common physico-chemical descriptors log K(ow) and log S(w). The chemical applicability domain of each model was verified by the leverage approach in order to propose only reliable data. The best predicted data were obtained by consensus modeling from 10 different models in the genetic algorithm model population.

In silico screening of estrogen-like chemicals based on different nonlinear classification models.

Increasing concern is being shown by the scientific community, government regulators, and the public about endocrine-disrupting chemicals that are adversely affecting human and wildlife health through a variety of mechanisms. There is a great need for an effective means of rapidly assessing endocrine-disrupting activity, especially estrogen-simulating activity, because of the large number of such chemicals in the environment. In this study, quantitative structure activity relationship (QSAR) models were developed to quickly and effectively identify possible estrogen-like chemicals based on 232 structurally-diverse chemicals (training set) by using several nonlinear classification methodologies (least-square support vector machine (LS-SVM), counter-propagation artificial neural network (CP-ANN), and k nearest neighbour (kNN)) based on molecular structural descriptors. The models were externally validated by 87 chemicals (prediction set) not included in the training set. All three methods can give satisfactory prediction results both for training and prediction sets, and the most accurate model was obtained by the LS-SVM approach through the comparison of performance. In addition, our model was also applied to about 58,000 discrete organic chemicals; about 76% were predicted not to bind to Estrogen Receptor. The obtained results indicate that the proposed QSAR models are robust, widely applicable and could provide a feasible and practical tool for the rapid screening of potential estrogens.