Prediction of chemical biodegradability using computational methods

Prediction of Drug-Induced Liver Injury: From Molecular Physicochemical Properties and Scaffold Architectures to Machine Learning Approaches

The process of developing new drugs is widely acknowledged as being time-intensive and requiring substantial financial investment. Despite ongoing efforts to reduce time and expenses in drug development, ensuring medication safety remains an urgent problem. One of the major problems involved in drug development is hepatotoxicity, specifically known as drug-induced liver injury (DILI). The popularity of new drugs often poses a significant barrier during development and frequently leads to their recall after launch. In silico methods have many advantages compared with traditional in vivo and in vitro assays. To establish a more precise and reliable prediction model, it is necessary to utilize an extensive and high-quality database consisting of information on drug molecule properties and structural patterns. In addition, we should also carefully select appropriate molecular descriptors that can be used to accurately depict compound characteristics. The aim of this study was to conduct a comprehensive investigation into the prediction of DILI. First, we conducted a comparative analysis of the physicochemical properties of extensively well-prepared DILI-positive and DILI-negative compounds. Then, we used classic substructure dissection methods to identify structural pattern differences between these two different types of chemical molecules. These findings indicate that it is not feasible to establish property or substructure-based rules for distinguishing between DILI-positive and DILI-negative compounds. Finally, we developed quantitative classification models for predicting DILI using the naïve Bayes classifier (NBC) and recursive partitioning (RP) machine learning techniques. The optimal DILI prediction model was obtained using NBC, which combines 21 physicochemical properties, the VolSurf descriptors and the LCFP_10 fingerprint set. This model achieved a global accuracy (GA) of 0.855 and an area under the curve (AUC) of 0.704 for the training set, while the corresponding values were 0.619 and 0.674 for the test set, respectively. Moreover, indicative substructural fragments favorable or unfavorable for DILI were identified from the best naïve Bayesian classification model. These findings may help prioritize lead compounds in the early stage of drug development pipelines.

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.

An improved opposition-based crow search algorithm for biodegradable material classification

The development of a reliable quantitative structure-activity relationship (QSAR) classification model with a small number of molecular descriptors is a crucial step in chemometrics. In this study, an improvement of crow search algorithm (CSA) is proposed by adapting the opposite-based learning (OBL) approach, which is named as OBL-CSA, to improve the exploration and exploitation capability of the CSA in quantitative structure-biodegradation relationship (QSBR) modelling of classifying the biodegradable materials. The results reveal that the performance of OBL-CSA not only manifest in improving the classification performance, but also in reduced computational time required to complete the process when compared to the standard CSA and other four optimization algorithms tested, which are the particle swarm algorithm (PSO), black hole algorithm (BHA), grey wolf algorithm (GWA), and whale optimization algorithm (WOA). In conclusion, the OBL-CSA could be a valuable resource in the classification of biodegradable materials.

In silico analytical toolset for predictive degradation and toxicity of hazardous pollutants in water sources

Different phenolic compounds, including multimeric lignin derivatives in the β-O-4 form, are among the most prevalent compounds in wastewater, often generated from paper industries. Relatively small concentrations of lignin are hazardous to aquatic organisms and can trigger severe environmental hazards. Herein, we present a predictive toolset to insight the induced toxic hazards prediction, and their Lignin peroxidase (LiP)-assisted degradation mechanism of selected multimeric lignin model compounds. T.E.ST and Toxtree toolset were deployed for toxic hazards estimation in different endpoints. To minimize the concerning hazards, we screened multimeric compounds for binding affinity with LiP. The binding affinity was found to be significantly lower than the reference compound. An Extra precision (XP) Glide score of -6.796 kcal/mol was found for dimer (guaiacyl 4-O-5 guaiacyl) complex as lowest compared to reference compound (-4.007 kcal/mol). The active site residues ASP-153, HIP-226, VAL-227, ARG-244, GLU-215, 239, PHE-261 were identified as site-specific key binding AA residues actively involved with corresponding ligands, forming Hydrophobic, H-Bond, π-Stacking, π-π type interactions. The DESMOND-assisted molecular dynamics simulation's (MDS) trajectories of protein-ligand revealed the considerable binding behavior and attained stability and system equilibrium state. Such theoretical and predictive conclusions indicted the feasibility of LiP assisted sustainable mitigation of lignin-based compounds, and such could be used to protect the environment from the potential hazards posed by recognized similar pollutants.

A Comparative Study of the Performance for Predicting Biodegradability Classification: The Quantitative Structure-Activity Relationship Model vs the Graph Convolutional Network

The prediction and evaluation of the biodegradability of molecules with computational methods are becoming increasingly important. Among the various methods, quantitative structure-activity relationship (QSAR) models have been demonstrated to predict the ready biodegradation of chemicals but have limited functionality owing to their complex implementation. In this study, we employ the graph convolutional network (GCN) method to overcome these issues. A biodegradability dataset from previous studies was trained to generate prediction models by (i) the QSAR models using the Mordred molecular descriptor calculator and MACCS molecular fingerprint and (ii) the GCN model using molecular graphs. The performance comparison of the methods confirms that the GCN model is more straightforward to implement and more stable; the specificity and sensitivity values are almost identical without specific descriptors or fingerprints. In addition, the performance of the models was further verified by randomly dividing the dataset into 100 different cases of training and test sets and by varying the test set ratio from 20 to 80%. The results of the current study clearly suggest the promise of the GCN model, which can be implemented straightforwardly and can replace conventional QSAR prediction models for various types and properties of molecules.

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.

Classification of Biodegradable Substances Using Balanced Random Trees and Boosted C5.0 Decision Trees

Substances that do not degrade over time have proven to be harmful to the environment and are dangerous to living organisms. Being able to predict the biodegradability of substances without costly experiments is useful. Recently, the quantitative structure-activity relationship (QSAR) models have proposed effective solutions to this problem. However, the molecular descriptor datasets usually suffer from the problems of unbalanced class distribution, which adversely affects the efficiency and generalization of the derived models. Accordingly, this study aims at validating the performances of balanced random trees (RTs) and boosted C5.0 decision trees (DTs) to construct QSAR models to classify the ready biodegradation of substances and their abilities to deal with unbalanced data. The balanced RTs model algorithm builds individual trees using balanced bootstrap samples, while the boosted C5.0 DT is modeled using cost-sensitive learning. We employed the two-dimensional molecular descriptor dataset, which is publicly available through the University of California, Irvine (UCI) machine learning repository. The molecular descriptors were ranked according to their contributions to the balanced RTs classification process. The performance of the proposed models was compared with previously reported results. Based on the statistical measures, the experimental results showed that the proposed models outperform the classification results of the support vector machine (SVM), K-nearest neighbors (KNN), and discrimination analysis (DA). Classification measures were analyzed in terms of accuracy, sensitivity, specificity, precision, false positive rate, false negative rate, F1 score, receiver operating characteristic (ROC) curve, and area under the ROC curve (AUROC).

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.

Evaluation of phenylethylamine type entactogens and their metabolites relevant to ecotoxicology - a QSAR study

The impact of the selected entactogens and their o-quinone metabolites on the environment was explored in QSAR studies by the use of predicted molecular descriptors, ADMET properties and environmental toxicity parameters, i.e., acute toxicity in Tetrahymena pyriformis (TOX_ATTP) expressed as Th_pyr_pIGC50/mmol L-1, acute toxicity in Pimephales promelas, the fathead minnow (TOX_FHM) expressed as Minnow LC50/mg L-1, the acute toxicity in Daphnia magna (TOX_DM) expressed as Daphnia LC50/mg L-1 and bioconcentration factor (BCF). The formation of corresponding o-quinones via benzo-dioxo-lone ring, O-demethylenation was predicted as the main metabolic pathway for all entactogens except for 1-(2,2-difluorobenzo[d][1,3]dioxol-5-yl)propan-2-amine (DiFMDA). The least favourable ADMET profile was revealed for N-(1-(benzo[d][1,3]dioxol-5-yl)propan-2-yl)-O-methylhydroxylamine (MDMEO). QSAR studies revealed significant linear correlations between MlogP of entactogens and MlogP of o-quinone metabolites (R = 0.99), and Th_pyr_pIGC50/mmol L-1 (R = 0.94), also their MlogPs with Minnow_LC50/mg L-1 (R = 0.80 and R = 0.78), BCF (R = 0.86 and R = 0.82) and percentage of o-quinones' yields (R = 0.73 and R = 0.80). Entactogens were predicted as non-biodegradable molecules, whereas the majority of their o-quinones were biodegradable.