TyPol - a new methodology for organic compounds clustering based on their molecular characteristics and environmental behavior

Solid-liquid distribution of SARS-CoV-2 in primary effluent of a wastewater treatment plant

Distributions of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) and fecal viral biomarkers between solid and liquid phases of wastewater are largely unknown. Herein, distributions of SARS-CoV-2, Pepper Mild Mottle Virus (PMMoV), and F-RNA bacteriophage group II (FRNAPH-II) were determined by viral RNA RT-qPCR. Comparison of viral recovery using three conventional fractionation methods included membrane filtration, a combination of mid-speed centrifugation and membrane filtration, and high-speed centrifugation. SARS-CoV-2 partitioned to the solids fraction in greater abundance compared to liquid fractions in a combination of mid-speed centrifugation and membrane filtration and high-speed centrifugation, but not in membrane filtration method in a particular assay, while fecal biomarkers (PMMoV and FRNAPH-II) exhibited the reciprocal relationship. The wastewater fractionation method had minimal effects on the solids-liquids distribution for all viral and phage markers tested; however, viral RNA load was significantly greater in solid-liquid fractions viral RNA loads compared with the than whole-wastewater PEG precipitation. A RNeasy PowerWater Kit with PCR inhibitor removal resulted in greater viral RNA loads and lesser PCR inhibition compared to a QIAamp Viral RNA Mini Kit without PCR inhibitor removal. These results support the development of improved methods and interpretation of WBE of SARS-CoV-2. •Distribution of SARS-CoV-2 to liquid and solid portions was addressed.•Addressing PCR inhibition is important in wastewater-based epidemiology.•Fraction methods have minimal effect.

Consideration of unmeasured micropollutants released from WWTP for potential impact estimations

During wastewater treatment, micropollutants are only partly eliminated and may present a risk for human health and aquatic ecosystems. The potential impacts these substances may have are currently underestimated due to the lack in available concentrations that lie below the limit of quantification (LOQ) for an important set of micropollutants. Here, the potential impacts due to 261 organic micropollutants on human health and aquatic environments were investigated at the scale of France. Even with concentrations below the LOQ, certain micropollutants were found to have a significant potential impact. For unmeasured concentrations, a global concentration distribution built from several datasets with different LOQ was used. By disregarding the unmeasured micropollutants, the potential impacts have been underestimated by >300 % on both human health and aquatic environments. Certain substances, such as hydrazine, endrin, or 2,3,7,8-TetraCDD, could lead to very strong potential impacts, even with unmeasured concentration levels. Moreover, the usual convention of LOQ/2 to replace unmeasured concentrations also appeared to overestimate the potential impact. The present work can be adapted to any other compartment or geographical context.

Assessment of pesticides volatilization potential based on their molecular properties using the TyPol tool

Following treatment, amounts of pesticides can reach the atmosphere because of spray drift, volatilization from soil or plants, and/or wind erosion. Monitoring and risk assessment of air contamination by pesticides is a recent issue and more insights on pesticide transfer to atmosphere are needed. Thus, the objective of this work was to better understand and assess pesticides emission potential to air through volatilization. The TyPol tool was used to explore the relationships between the global, soil and plant volatilization potentials of 178 pesticides, and their molecular properties. The outputs of TyPol were then compared to atmospheric pesticide concentrations monitored in various French regions. TyPol was able to discriminate pesticides that were observed in air from those that were not. Clustering considering parameters driving the emission potential from soil (sorption characteristics) or plant (lipophilic properties), in addition to vapor pressure, allowed better discrimination of the pesticides than clustering considering all parameters for the global emission potential. Pesticides with high volatilization potential have high total energy, and low molecular weight, molecular connectivity indices and polarizability. TyPol helped better understand the volatilization potential of pesticides. It can be used as a first step to assess the risk of air contamination by pesticides.

Suspect screening of environmental contaminants by UHPLC-HRMS and transposable Quantitative Structure-Retention Relationship modelling

A Quantitative Structure-Retention Relationship (QSRR) model is proposed and aims at increasing the confidence level associated to the identification of organic contaminants by Ultra-High Performance Liquid Chromatography hyphenated to High Resolution Mass Spectrometry (UHPLC-HRMS) in environmental samples under a suspect screening approach. The model was built from a selection of 8 easily accessible physicochemical descriptors, and was validated from a set of 274 organic compounds commonly found in environmental samples. The proposed predictive figure approach is based on the mobile phase composition at solute elution (expressed as % acetonitrile), that has the major advantage of making the model reusable by other laboratories, since the elution composition is independent of both the column geometry and the UHPLC-system. The model quality was assessed and was altered neither by the columns from different lots, nor by the complex matrices of environmental water samples. Then, the solute retention of any organic compound present in water samples is expected to be predicted within ± 14.3% acetonitrile by our model. Solute retention can therefore be used as a supplementary tool for the identification of environmental contaminants by UHPLC-HRMS, in addition to mass spectrometry data already used in the suspect screening approach.

Influence of Glyphosate Formulations on the Behavior of Sulfentrazone in Soil in Mixed Applications

The selection of weed biotypes that are resistant to glyphosate has increased the demand for its use mixed with other herbicides, such as sulfentrazone. However, when chemical molecules are mixed, interactions may occur, modifying the behavior of these molecules in the environment, such as the sorption and desorption in soil. In this study, we hypothesized that the presence of glyphosate-formulated products might increase the sorption or decrease the desorption of sulfentrazone, thereby increasing the risk of the contamination of water resources. Therefore, our work aimed to evaluate the sorption, desorption, and leaching of sulfentrazone in the soil in an isolated and mixed application with different glyphosate formulations. The sorption coefficients (Kfs) for the sulfentrazone, sulfentrazone + Roundup Ready, sulfentrazone + Roundup Ultra, and sulfentrazone + Zapp Qi foram were 1.3, 2.1, 2.3, and 1.9, respectively. The desorption coefficients (Kfd) for the sulfentrazone, sulfentrazone + Roundup Ready, sulfentrazone + Roundup Ultra, and sulfentrazone + Zapp Qi foram were 65.7, 125.2, 733.3 and 239.8, respectively. The experiments demonstrated that the sorption and desorption of sulfentrazone in combination with the other formulated glyphosate products are altered, supporting the hypothesis suggested by this work, i.e., that the presence of other molecules is a factor that affects the behavior of herbicides in the soil. This phenomenon altered the vertical mobility of sulfentrazone. Situations involving mixtures of pesticides should be evaluated in order to improve our understanding of the dynamics of these molecules and thus avoid environmental contamination.

Natural Chlordecone Degradation Revealed by Numerous Transformation Products Characterized in Key French West Indies Environmental Compartments

Production and use of the insecticide chlordecone has caused long-term environmental pollution in the James River area and the French West Indies (FWI) that has resulted in acute human-health problems and a social crisis. High levels of chlordecone in FWI soils, even after its ban decades ago, and the absence of detection of transformation products (TPs), have suggested that chlordecone is virtually nonbiodegradable in the environment. Here, we investigated laboratory biodegradation, consisting of bacterial liquid cultures and microcosms inoculated with FWI soils, using a dual nontargeted GC-MS and LC-HRMS approach. In addition to previously reported, partly characterized hydrochlordecones and polychloroindenes (families A and B), we discovered 14 new chlordecone TPs, assigned to four families (B, C, D, and E). Organic synthesis and NMR analyses allowed us to achieve the complete structural elucidation of 19 TPs. Members of TP families A, B, C, and E were detected in soil, sediment, and water samples from Martinique and include 17 TPs not initially found in commercial chlordecone formulations. 2,4,5,6,7-Pentachloroindene was the most prominent TP, with levels similar to those of chlordecone. Overall, our results clearly show that chlordecone pollution extends beyond the parent chlordecone molecule and includes a considerable number of previously undetected TPs. Structural diversity of the identified TPs illustrates the complexity of chlordecone degradation in the environment and raises the possibility of extensive worldwide pollution of soil and aquatic ecosystems by chlordecone TPs.

Interaction of Metamitron and Fenhexamid with Ca2+ -Montmorillonite Clay Surfaces: A Density Functional Theory Molecular Dynamics Study

Metamitron (Meta), an herbicide, and fenhexamid (Fen), a fungicide, are authorized by the European Union to be used in agriculture. This article reports theoretical calculations about Meta and Fen in interaction with a clay surface: a Ca-montmorillonite (Mont). Conformational searches have been performed thanks to Car-Parrinello molecular dynamics simulations from which geometries have been extracted. Interaction and adsorption energies have been calculated for isomers of Meta or Fen in interaction with Mont to understand the relative stability of various kinds of complexation. Substantial adsorption energies are comparable for Meta and Fen: around -40 kcal/mol. For Fen-Mont, the CO monodentate family is surprisingly the lowest in energy. Moreover, the 10 lowest-energy isomers involve complexation on Fen carbonyl oxygens. The Meta-Mont lowest-energy family, N-N, does not involve π delocalization breaking within Meta. At the same time, the stronger the interaction energy is, the larger the structural modifications within Mont are, particularly concerning the interacting cation distance to the surface. The non-negligible charge transfer and the magnitude of the adsorption energy speak in favor of the chemisorption of the pesticide on the surface. © 2019 Wiley Periodicals, Inc.

Clustering pesticides according to their molecular properties, fate, and effects by considering additional ecotoxicological parameters in the TyPol method

Understanding the fate and ecotoxicological effects of pesticides largely depends on their molecular properties. We recently developed "TyPol" (Typology of Pollutants), a classification method of organic compounds based on statistical analyses. It combines several environmental (sorption coefficient, degradation half-life) and one ecotoxicological (bioconcentration factor) parameters, to structural molecular descriptors (number of atoms in the molecule, molecular surface, dipole moment, energy of orbitals, etc.). The present study attempts to extend TyPol to the ecotoxicological effects of pesticides on non-target organisms, based on data analysis from available literature and databases. It revealed that relevant ecotoxicological endpoints for terrestrial organisms (e.g., soil microorganisms, invertebrates) that support a range of ecosystemic services are lacking as compared to aquatic organisms. The availability of ecotoxicological parameters was also lower for chronic than for acute ecotoxicity endpoints. Consequently, seven parameters were included for acute (EC50, LC50) and chronic (NOEC) ecotoxicological effects for one terrestrial (Eisenia sp.) and three aquatic (Daphnia sp., algae, Lemna sp.) organisms. In this new configuration, we used TyPol to classify 50 pesticides into different clusters that gather molecules with similar environmental behaviors and ecotoxicological effects. The classification results evidenced relationships between molecular descriptors, environmental parameters, and the added ecotoxicological endpoints. This proof-of-concept study also showed that TyPol in silico classification can successfully address new scientific questions and be expanded with other parameters of interest.

Categorizing chlordecone potential degradation products to explore their environmental fate

Chlordecone (C₁₀Cl₁₀O; CAS number 143-50-0) has been used extensively as an organochlorine insecticide but is nowadays banned and listed on annex A in The Stockholm Convention on Persistent Organic Pollutants (POPs). Although experimental evidences of biodegradation of this compound are scarce, several dechlorination products have been proposed by Dolfing et al. (2012) using Gibbs free energy calculations to explore different potential transformation routes. We here present the results of an in silico classification (TyPol - Typology of Pollutants) of chlordecone transformation products (TPs) based on statistical analyses combining several environmental endpoints and structural molecular descriptors. Starting from the list of putative chlordecone TPs and considering available data on degradation routes of other organochlorine compounds, we used different clustering strategies to explore the potential environmental behaviour of putative chlordecone TPs from the knowledge on their molecular descriptors. The method offers the possibility to focus on TPs present in different classes and to infer their environmental fate. Thus, we have deduced some hypothetical trends for the environmental behaviour of TPs of chlordecone assuming that TPs, which were clustered away from chlordecone, would have different environmental fate and ecotoxicological impact compared to chlordecone. Our findings suggest that mono- and di-hydrochlordecone, which are TPs of chlordecone often found in contaminated soils, may have similar environmental behaviour in terms of persistence.

Towards a better pesticide policy for the European Union

This opinion article aims to foster the debate about pesticide legislation in the European Union (EU). Numerous formerly authorized and widely used pesticides are now banned in the EU because unexpected and unacceptable risks emerged after their initial introduction to the market. Throughout this time lapse, environmental quality and human health have been threatened by the use of these compounds. These hazards could have been prevented by a more responsive pesticide regulatory framework. This article provides detailed insights into the pros and cons of pesticides, and points out weaknesses of the current pesticide environmental risk assessment procedures. Possibilities for improving the robustness and reliability of the pesticide regulatory framework are discussed.

QSAR Study for Carcinogenic Potency of Aromatic Amines Based on GEP and MLPs

A new analysis strategy was used to classify the carcinogenicity of aromatic amines. The physical-chemical parameters are closely related to the carcinogenicity of compounds. Quantitative structure activity relationship (QSAR) is a method of predicting the carcinogenicity of aromatic amine, which can reveal the relationship between carcinogenicity and physical-chemical parameters. This study accessed gene expression programming by APS software, the multilayer perceptrons by Weka software to predict the carcinogenicity of aromatic amines, respectively. All these methods relied on molecular descriptors calculated by CODESSA software and eight molecular descriptors were selected to build function equations. As a remarkable result, the accuracy of gene expression programming in training and test sets are 0.92 and 0.82, the accuracy of multilayer perceptrons in training and test sets are 0.84 and 0.74 respectively. The precision of the gene expression programming is obviously superior to multilayer perceptrons both in training set and test set. The QSAR application in the identification of carcinogenic compounds is a high efficiency method.

Identification and characterization of tebuconazole transformation products in soil by combining suspect screening and molecular typology

Pesticides generate transformation products (TPs) when they are released into the environment. These TPs may be of ecotoxicological importance. Past studies have demonstrated how difficult it is to predict the occurrence of pesticide TPs and their environmental risk. The monitoring approaches mostly used in current regulatory frameworks target only known ecotoxicologically relevant TPs. Here, we present a novel combined approach which identifies and categorizes known and unknown pesticide TPs in soil by combining suspect screening time-of-flight mass spectrometry with in silico molecular typology. We used an empirical and theoretical pesticide TP library for compound identification by both non-target and target time-of-flight (tandem) mass spectrometry, followed by structural proposition through a molecular structure correlation program. In silico molecular typology was then used to group TPs according to common molecular descriptors and to indirectly elucidate their environmental parameters by analogy to known pesticide compounds with similar molecular descriptors. This approach was evaluated via the identification of TPs of the triazole fungicide tebuconazole occurring in soil during a field dissipation study. Overall, 22 empirical and 12 yet unknown TPs were detected, and categorized into three groups with defined environmental properties. This approach combining suspect screening time-of-flight mass spectrometry with molecular typology could be extended to other organic pollutants and used to rationalize the choice of TPs to be investigated towards a more comprehensive environmental risk assessment scheme.

Consensus toxicity factors for polychlorinated dibenzo-p-dioxins, dibenzofurans, and biphenyls combining in silico models and extensive in vitro screening of AhR-mediated effects in human and rodent cells

Consensus toxicity factors (CTFs) were developed as a novel approach to establish toxicity factors for risk assessment of dioxin-like compounds (DLCs). Eighteen polychlorinated dibenzo-p-dioxins, dibenzofurans (PCDD/Fs), and biphenyls (PCBs) with assigned World Health Organization toxic equivalency factors (WHO-TEFs) and two additional PCBs were screened in 17 human and rodent bioassays to assess their induction of aryl hydrocarbon receptor-related responses. For each bioassay and compound, relative effect potency values (REPs) compared to 2,3,7,8-tetrachlorodibenzo-p-dioxin were calculated and analyzed. The responses in the human and rodent cell bioassays generally differed. Most notably, the human cell models responded only weakly to PCBs, with 3,3',4,4',5-pentachlorobiphenyl (PCB126) being the only PCB that frequently evoked sufficiently strong responses in human cells to permit us to calculate REP values. Calculated REPs for PCB126 were more than 30 times lower than the WHO-TEF value for PCB126. CTFs were calculated using score and loading vectors from a principal component analysis to establish the ranking of the compounds and, by rescaling, also to provide numerical differences between the different congeners corresponding to the TEF scheme. The CTFs were based on rat and human bioassay data and indicated a significant deviation for PCBs but also for certain PCDD/Fs from the WHO-TEF values. The human CTFs for 2,3,4,7,8-pentachlorodibenzofuran, 1,2,3,4,7,8-hexachlorodibenzofuran, 1,2,3,4,6,7,8-heptachlorodibenzo-p-dioxin, and 1,2,3,4,7,8,9-heptachlorodibenzofuran were up to 10 times greater than their WHO-TEF values. Quantitative structure-activity relationship models were used to predict CTFs for untested WHO-TEF compounds, suggesting that the WHO-TEF value for 1,2,3,7,8-pentachlorodibenzofuran could be underestimated by an order of magnitude for both human and rodent models. Our results indicate that the CTF approach provides a powerful tool for condensing data from batteries of screening tests using compounds with similar mechanisms of action, which can be used to improve risk assessment of DLCs.