Quantitative structure-activity relationships for primary aerobic biodegradation of organic chemicals in pristine surface waters: starting points for predicting biodegradation under acclimatization

Improving predictions and understanding of primary and ultimate biodegradation rates with machine learning models

This study aimed to develop machine learning based quantitative structure biodegradability relationship (QSBR) models for predicting primary and ultimate biodegradation rates of organic chemicals, which are essential parameters for environmental risk assessment. For this purpose, experimental primary and ultimate biodegradation rates of high consistency were compiled for 173 organic compounds. A significant number of descriptors were calculated with a collection of quantum/computational chemistry software and tools to achieve comprehensive representation and interpretability. Following a pre-screening process, multiple QSBR models were developed for both primary and ultimate endpoints using three algorithms: extreme gradient boosting (XGBoost), support vector machine (SVM), and multiple linear regression (MLR). Furthermore, a unified QSBR model was constructed using the knowledge transfer technique and XGBoost. Results demonstrated that all QSBR models developed in this study had good performance. Particularly, SVM models exhibited high level of goodness of fit (coefficient of determination on the training set of 0.973 for primary and 0.980 for ultimate), robustness (leave-one-out cross-validated coefficient of 0.953 for primary and 0.967 for ultimate), and external predictive ability (external explained variance of 0.947 for primary and 0.958 for ultimate). The knowledge transfer technique enhanced model performance by learning from properties of two biodegradation endpoints. Williams plots were used to visualize the application domains of the models. Through SHapley Additive exPlanations (SHAP) analysis, this study identified key features affecting biodegradation rates. Notably, MDEO-12, APC2D1_C_O, and other features contributed to primary biodegradation, while AATS0v, AATS2v, and others inhibited it. For ultimate biodegradation, features like No. of Rotatable Bonds, APC2D1_C_O, and minHBa were contributors, while C1SP3, Halogen Ratio, GGI4, and others hindered the process. Also, the study quantified the contributions of each feature in predictions for individual chemicals. This research provides valuable tools for predicting both primary and ultimate biodegradation rates while offering insights into the mechanisms.

Scientific concepts and methods for moving persistence assessments into the 21st century

The evaluation of a chemical substance's persistence is key to understanding its environmental fate, exposure concentration, and, ultimately, environmental risk. Traditional biodegradation test methods were developed many years ago for soluble, nonvolatile, single-constituent test substances, which do not represent the wide range of manufactured chemical substances. In addition, the Organisation for Economic Co-operation and Development (OECD) screening and simulation test methods do not fully reflect the environmental conditions into which substances are released and, therefore, estimates of chemical degradation half-lives can be very uncertain and may misrepresent real environmental processes. In this paper, we address the challenges and limitations facing current test methods and the scientific advances that are helping to both understand and provide solutions to them. Some of these advancements include the following: (1) robust methods that provide a deeper understanding of microbial composition, diversity, and abundance to ensure consistency and/or interpret variability between tests; (2) benchmarking tools and reference substances that aid in persistence evaluations through comparison against substances with well-quantified degradation profiles; (3) analytical methods that allow quantification for parent and metabolites at environmentally relevant concentrations, and inform on test substance bioavailability, biochemical pathways, rates of primary versus overall degradation, and rates of metabolite formation and decay; (4) modeling tools that predict the likelihood of microbial biotransformation, as well as biochemical pathways; and (5) modeling approaches that allow for derivation of more generally applicable biotransformation rate constants, by accounting for physical and/or chemical processes and test system design when evaluating test data. We also identify that, while such advancements could improve the certainty and accuracy of persistence assessments, the mechanisms and processes by which they are translated into regulatory practice and development of new OECD test guidelines need improving and accelerating. Where uncertainty remains, holistic weight of evidence approaches may be required to accurately assess the persistence of chemicals. Integr Environ Assess Manag 2022;18:1454-1487. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Classification and Regression Machine Learning Models for Predicting Aerobic Ready and Inherent Biodegradation of Organic Chemicals in Water

Machine learning (ML) is viewed as a promising tool for the prediction of aerobic biodegradation, one of the most important elimination pathways of organic chemicals from the environment. However, available models only have small datasets (<3200 records), make binary classification predictions, evaluate ready biodegradability, and do not incorporate experimental conditions (e.g., system setup and reaction time). This study addressed all these limitations by first compiling a large database of 12,750 records, considering both ready and inherent biodegradation under different conditions, and then developing regression and classification models using different chemical representations and ML algorithms. The best regression model (R² = 0.54 and root mean square error of 0.25) and classification model (the prediction accuracy from 85.1%) achieved very good performance. The model interpretation indicated that the models correctly captured the effects of chemical substructures, following the order of C═O > O═C-O > OH > CH₃ > halogen > branching > N > 6-member ring. The consideration of chemical speciation based on pK_a and α notations did not affect the regression model performance but significantly improved the classification model performance (the accuracy increased to 87.6%). The models also showed large applicability domains and provided reasonable predictions for more than 98% of over 850,000 environmentally relevant chemicals in the Distributed Structure-Searchable Toxicity database. These robust, trustable models were finally made widely accessible through two free online predictors with graphical user interface.

A review on effective soil health bio-indicators for ecosystem restoration and sustainability

Preventing degradation, facilitating restoration, and maintaining soil health is fundamental for achieving ecosystem stability and resilience. A healthy soil ecosystem is supported by favorable components in the soil that promote biological productivity and provide ecosystem services. Bio-indicators of soil health are measurable properties that define the biotic components in soil and could potentially be used as a metric in determining soil functionality over a wide range of ecological conditions. However, it has been a challenge to determine effective bio-indicators of soil health due to its temporal and spatial resolutions at ecosystem levels. The objective of this review is to compile a set of effective bio-indicators for developing a better understanding of ecosystem restoration capabilities. It addresses a set of potential bio-indicators including microbial biomass, respiration, enzymatic activity, molecular gene markers, microbial metabolic substances, and microbial community analysis that have been responsive to a wide range of ecosystem functions in agricultural soils, mine deposited soil, heavy metal contaminated soil, desert soil, radioactive polluted soil, pesticide polluted soil, and wetland soils. The importance of ecosystem restoration in the United Nations Sustainable Development Goals was also discussed. This review identifies key management strategies that can help in ecosystem restoration and maintain ecosystem stability.

Identifying Functional Groups that Determine Rates of Micropollutant Biotransformations Performed by Wastewater Microbial Communities

The goal of this research was to identify functional groups that determine rates of micropollutant (MP) biotransformations performed by wastewater microbial communities. To meet this goal, we performed a series of incubation experiments seeded with four independent wastewater microbial communities and spiked them with a mixture of 40 structurally diverse MPs. We collected samples over time and used high-resolution mass spectrometry to estimate biotransformation rate constants for each MP in each experiment and to propose structures of 46 biotransformation products. We then developed random forest models to classify the biotransformation rate constants based on the presence of specific functional groups or observed biotransformations. We extracted classification importance metrics from each random forest model and compared them across wastewater microbial communities. Our analysis revealed 30 functional groups that we define as either biotransformation promoters, biotransformation inhibitors, structural features that can be biotransformed based on uncharacterized features of the wastewater microbial community, or structural features that are not rate-determining. Our experimental data and analysis provide novel insights into MP biotransformations that can be used to more accurately predict MP biotransformations or to inform the design of new chemical products that may be more readily biodegradable during wastewater treatment.

Microbial Technologies Employed for Biodegradation of Neonicotinoids in the Agroecosystem

Neonicotinoids are synthetic pesticides widely used for the control of various pests in agriculture throughout the world. They mainly attack the nicotinic acetylcholine receptors, generate nervous stimulation, receptor clot, paralysis and finally cause death. They are low volatile, highly soluble and have a long half-life in soil and water. Due to their extensive use, the environmental residues have immensely increased in the last two decades and caused many hazardous effects on non-target organisms, including humans. Hence, for the protection of the environment and diversity of living organism's the degradation of neonicotinoids has received widespread attention. Compared to the other methods, biological methods are considered cost-effective, eco-friendly and most efficient. In particular, the use of microbial species makes the degradation of xenobiotics more accessible fast and active due to their smaller size. Since this degradation also converts xenobiotics into less toxic substances, the various metabolic pathways for the microbial degradation of neonicotinoids have been systematically discussed. Additionally, different enzymes, genes, plasmids and proteins are also investigated here. At last, this review highlights the implementation of innovative tools, databases, multi-omics strategies and immobilization techniques of microbial cells to detect and degrade neonicotinoids in the environment.

Methodological Advances to Study Contaminant Biotransformation: New Prospects for Understanding and Reducing Environmental Persistence?

Complex microbial communities in environmental systems play a key role in the detoxification of chemical contaminants by transforming them into less active metabolites or by complete mineralization. Biotransformation, i.e., transformation by microbes, is well understood for a number of priority pollutants, but a similar level of understanding is lacking for many emerging contaminants encountered at low concentrations and in complex mixtures across natural and engineered systems. Any advanced approaches aiming to reduce environmental exposure to such contaminants (e.g., novel engineered biological water treatment systems, design of readily degradable chemicals, or improved regulatory assessment strategies to determine contaminant persistence a priori) will depend on understanding the causal links among contaminant removal, the key driving agents of biotransformation at low concentrations (i.e., relevant microbes and their metabolic activities), and how their presence and activity depend on environmental conditions. In this Perspective, we present the current understanding and recent methodological advances that can help to identify such links, even in complex environmental microbiomes and for contaminants present at low concentrations in complex chemical mixtures. We discuss the ensuing insights into contaminant biotransformation across varying environments and conditions and ask how much closer we have come to designing improved approaches to reducing environmental exposure to contaminants.

Effect of concentration and hydraulic reaction time on the removal of pharmaceutical compounds in a membrane bioreactor inoculated with activated sludge

Pharmaceuticals are often not fully removed in wastewater treatment plants (WWTPs) and are thus being detected at trace levels in water bodies all over the world posing a risk to numerous organisms. These organic micropollutants (OMPs) reach WWTPs at concentrations sometimes too low to serve as growth substrate for microorganisms; thus, co-metabolism is thought to be the main conversion mechanism. In this study, the microbial removal of six pharmaceuticals was investigated in a membrane bioreactor at increasing concentrations (4-800 nM) of the compounds and using three different hydraulic retention times (HRT; 1, 3.5 and 5 days). The bioreactor was inoculated with activated sludge from a municipal WWTP and fed with ammonium, acetate and methanol as main growth substrates to mimic co-metabolism. Each pharmaceutical had a different average removal efficiency: acetaminophen (100%) > fluoxetine (50%) > metoprolol (25%) > diclofenac (20%) > metformin (15%) > carbamazepine (10%). Higher pharmaceutical influent concentrations proportionally increased the removal rate of each compound, but surprisingly not the removal percentage. Furthermore, only metformin removal improved to 80-100% when HRT or biomass concentration was increased. Microbial community changes were followed with 16S rRNA gene amplicon sequencing in response to the increment of pharmaceutical concentration: Nitrospirae and Planctomycetes 16S rRNA relative gene abundance decreased, whereas Acidobacteria and Bacteroidetes increased. Remarkably, the Dokdonella genus, previously implicated in acetaminophen metabolism, showed a 30-fold increase in abundance at the highest concentration of pharmaceuticals applied. Taken together, these results suggest that the incomplete removal of most pharmaceutical compounds in WWTPs is dependent on neither concentration nor reaction time. Accordingly, we propose a chemical equilibrium or a growth substrate limitation as the responsible mechanisms of the incomplete removal. Finally, Dokdonella could be the main acetaminophen degrader under activated sludge conditions, and non-antibiotic pharmaceuticals might still be toxic to relevant WWTP bacteria.

Trends in predictive biodegradation for sustainable mitigation of environmental pollutants: Recent progress and future outlook

The feasibility of in-silico techniques, together with the computational framework, has been applied to predictive bioremediation aiming to clean-up contaminants, toxicity evaluation, and possibilities for the degradation of complex recalcitrant compounds. Emerging contaminants from different industries have posed a significant hazard to the environment and public health. Given current bioremediation strategies, it is often a failure or inadequate for sustainable mitigation of hazardous pollutants. However, clear-cut vital information about biodegradation is quite incomplete from a conventional remediation techniques perspective. Lacking complete information on bio-transformed compounds leads to seeking alternative methods. Only scarce information about the transformed products and toxicity profile is available in the published literature. To fulfill this literature gap, various computational or in-silico technologies have emerged as alternating techniques, which are being recognized as in-silico approaches for bioremediation. Molecular docking, molecular dynamics simulation, and biodegradation pathways predictions are the vital part of predictive biodegradation, including the Quantitative Structure-Activity Relationship (QSAR), Quantitative structure-biodegradation relationship (QSBR) model system. Furthermore, machine learning (ML), artificial neural network (ANN), genetic algorithm (GA) based programs offer simultaneous biodegradation prediction along with toxicity and environmental fate prediction. Herein, we spotlight the feasibility of in-silico remediation approaches for various persistent, recalcitrant contaminants while traditional bioremediation fails to mitigate such pollutants. Such could be addressed by exploiting described model systems and algorithm-based programs. Furthermore, recent advances in QSAR modeling, algorithm, and dedicated biodegradation prediction system have been summarized with unique attributes.

Development of models predicting biodegradation rate rating with multiple linear regression and support vector machine algorithms

Biodegradation is a significant process for removing organic chemicals from water, soil and sediment environments, and therefore biodegradability is critical to evaluate the environmental persistence of organic chemicals. In this study, based on a dataset with 171 compounds, four quantitative structure-activity relationship (QSAR) models were developed for predicting primary and ultimate biodegradation rate rating with multiple linear regression (MLR) and support vector machine (SVM) algorithms. Two MLR models were built with a dataset with carbon atom number ≤9, and two SVM models were built with a dataset with carbon atom number >9. In the MLR models, n_ArX (number of X on aromatic ring) is the most important descriptor governing primary and ultimate biodegradation of organic chemicals. For the SVM models, determination coefficient (R²) values, cross-validated coefficients (Q²_LOO) and external validation coefficient (Q²_ext) values are over 0.9, indicating the SVM models have satisfactory goodness-of-fit, robustness and external predictive abilities. The applicability domains of these models were visualized by the Williams plot. The developed models can be used as effective tools to predict biodegradability of organic chemicals.

Comparison of Small Molecule Biotransformation Half-Lives between Activated Sludge and Soil: Opportunities for Read-Across?

Compartment-specific degradation half-lives are essential pieces of information in the regulatory risk assessment of synthetic chemicals. However, their measurement according to regulatory testing guidelines is laborious and costly. Despite the obvious ecological and economic benefits of knowing environmental degradability as early as possible, its consideration in the early phases of rational chemical design is therefore challenging. Here, we explore the possibility to use half-lives determined in highly time- and work-efficient biotransformation experiments with activated sludge and mixtures of chemicals to predict soil half-lives from regulatory simulation studies. We experimentally determined half-lives for 52 structurally diverse agrochemical active ingredients in batch reactors with three concentrations of the same activated sludge. We then developed bi- and multivariate models for predicting half-lives in soil by regressing the experimentally determined half-lives in activated sludge against average soil half-lives of the same chemicals extracted from regulatory data. The models differed in how we accounted for sorption-related bioavailability differences in soil and activated sludge. The best-performing models exhibited good coefficients of determination (R² of around 0.8) and low average errors (<factor of 3 in half-life predictions) and were robust in cross-validation. From a practical perspective, these results suggest that it may indeed be possible to read across from half-lives determined in highly efficient biotransformation experiments in activated sludge to soil half-lives, which are obtained from much more work- and resource-intense regulatory studies, and that these predictions are clearly superior to predictions based on the output of BIOWIN, a publicly available quantitative structure-biodegradation relationship (QSBR) model. From a theoretical perspective, these results suggest that soil and activated sludge microbial communities, although certainly different in terms of taxonomic composition, may be functionally similar with respect to the enzymatic transformation of environmentally relevant concentrations of a diverse range of chemical compounds.

Disentanglement of the chemical, physical, and biological processes aids the development of quantitative structure-biodegradation relationships for aerobic wastewater treatment

Attenuation of organic compounds in sewage treatment plants (STPs) is affected by a complex interplay between chemical (e.g. ionization, hydrolysis), physical (e.g. sorption, volatilization), and biological (e.g. biodegradation, microbial acclimation) processes. These effects should be accounted for individually, in order to develop predictive cheminformatics tools for STPs. Using measured data from 70 STPs in the Netherlands for 69 chemicals (pharmaceuticals, herbicides, etc.), we highlighted the influences of 1) chemical ionization, 2) sorption to sludge, and 3) acclimation of the microbial consortia on the primary removal of chemicals. We used semi-empirical corrections for each of these influences to deduce biodegradation rate constants upon which quantitative structure-biodegradation relationships (QSBRs) were developed. As shown by a global QSBR, biodegradation in STPs generally relates to structural complexity, size, energetics, and charge distribution. Statistics of the global QSBR were reasonable, being R²_training=0.69 (training set of 51 compounds) and R²_validation=0.50 (validation set of 18 compounds). Class-specific QSBRs utilized electronic properties potentially relating to rate-limiting enzymatic steps. For class-specific QSBRs, values of R² of in between 0.7 and 0.8 were obtained. With caution, environmental risk assessment methodologies may apply these models to estimate biodegradation rates for 'data-poor' compounds. The approach also highlights 'meta data' on STP operational parameters needed to develop QSBRs of better predictability in the future.

Gene Editing and Systems Biology Tools for Pesticide Bioremediation: A Review

Bioremediation is the degradation potential of microorganisms to dissimilate the complex chemical compounds from the surrounding environment. The genetics and biochemistry of biodegradation processes in datasets opened the way of systems biology. Systemic biology aid the study of interacting parts involved in the system. The significant keys of system biology are biodegradation network, computational biology, and omics approaches. Biodegradation network consists of all the databases and datasets which aid in assisting the degradation and deterioration potential of microorganisms for bioremediation processes. This review deciphers the bio-degradation network, i.e., the databases and datasets (UM-BBD, PAN, PTID, etc.) aiding in assisting the degradation and deterioration potential of microorganisms for bioremediation processes, computational biology and multi omics approaches like metagenomics, genomics, transcriptomics, proteomics, and metabolomics for the efficient functional gene mining and their validation for bioremediation experiments. Besides, the present review also describes the gene editing tools like CRISPR Cas, TALEN, and ZFNs which can possibly make design microbe with functional gene of interest for degradation of particular recalcitrant for improved bioremediation.